A general model for the effective viscosity of pseudoplastic and dilatant fluids

Summary A new and very general expression is proposed for correlation of data for the effective viscosity of pseudoplastic and dilatant fluids as a function of the shear stress. Most of the models which have been proposed previously are shown to be special cases of this expression. A straightforward procedure is outlined for evaluation of the arbitrary constants.

---

Zusammenfassung Eine neue und sehr allgemeine Formel wird für die Korrelation der Werte der effektiven Viskosität von strukturviskosen und dilatanten Flüssigkeiten in Abhängigkeit von der Schubspannung vorgeschlagen. Die meisten schon früher vorgeschlagenen Methoden werden hier als Spezialfälle dieser Gleichung gezeigt. Ein einfaches Verfahren für die Auswertung der willkürlichen Konstanten wird beschrieben.

---

Nomenclature b arbitrary constant inSisko model (eq. [5]) - n arbitrary exponent in eq. [1] - x independent variable - y(x) dependent variable - y ₀(x) limiting behavior of dependent variable asx 0 - y(x) limiting behavior of dependent variable asx - z original dependent variable - arbitrary constant inSisko model (eq. [5]) andBird-Sisko model (eq. [6]) - arbitrary exponent in eqs. [2] and [8] - effective viscosity = shear stress/rate of shear - _A effective viscosity at = _A - _B empirical constant in eqs. [2] and [8] - ₀ limiting value of effective viscosity as 0 - ₀() limiting behavior of effective viscosity as 0 - limiting value of effective viscosity as - () limiting behavior of effective viscosity as - rate of shear - arbitrary constant inBird-Sisko model (eq.[6]) - shear stress - _A arbitrary constant in eqs. [2] and [8] - ₀ shear stress at inBingham model - _1/2 shear stress at = ( ₀ + )/2 With 8 figures     

Radiative melting of a horizontal clear ice layer

In this paper the horizontal layer of clear ice sticking to the substrate is melted by comparatively short wave radiation similar to solar radiation for the purpose of removing ice from the surface of the material subject to atmospheric icing. The radiating source used for melting is 300 wattages halogen lamps whose color temperature is 3200K at 100 voltages. From the present investigation, a typical phenomenon of backmelting is observed clearly and it can be found that the predicted results including the melting rate of upper and lower layers which are melted by radiant energy impinged on or penetrated the ice layer are in good agreement with the experimental results.

---

Aufschmelzen einer waagerechten Klareisschicht durch Strahlung

Zusammenfassung Eine waagerechte Klareisschicht, die auf einer Unterlage aufgefroren war, wurde durch kurzwellige Strahlung, Ähnlich der Sonnenstrahlung, zum Schmelzen gebracht, um die Entfernung von Eis nach atmosphÄrischer Vereisung zu untersuchen. Die Strahlungsquelle war eine 300 Watt-Halogenlampe mit einer Farbtemperatur von 3200 Kelvin bei 100 Volt. Als typische Erscheinung wurde ein Rückseiten-Schmelzen gefunden, im übrigen sind die vorausberechneten Schmelzraten an der Ober- und der Unterseite durch aufgenommene oder durchgelassene Strahlungsenergie in guter übereinstimmung mit den Messungen.

---

Nomenclature a_v monochromatic absorption coefficient - A transmission (= q _r ⁺ {h_i}/q_ro) - c_p specific heat - E_bv monochromatic emissive power - h_D mass transfer coefficient - h_i initial thickness of ice layer - h_m thickness of substrate - L_i latent heat of melting - L_w latent heat of evaporation or condensation - heat flux absorbed at surface of substrate - q_r0 radiant heat flux impinged onto ice or free surface - q _r ⁺ {y} forward radiant heat flux - q _r ^– {y} backward radiant heat flux - S₁ thickness of upper melt layer - S₂ thickness of lower melt layer - S'₂ distance from free surface to bottom surface of ice layer - t time - T temperature - T₁ temperature of air-water or air-ice interface - T₂ temperature of substrate surface - T₃ temperature of back side surface of substrate - T_b temperature of radiating source - T_i temperature in ice layer - T_w1 temperature in upper melt layer - T_w2 temperature in lower melt layer - T environmental temperature - W_w saturated vapor concentration at free surface - W_t8 vapor concentration at environment - y distance from free or ice surface - y grid size of water or ice - y_m grid size of substrate Greek symbols heat transfer coefficient - spectral absorptivity - _t total absorptivity - _i thermal diffusivity of ice - _m thermal diffusivity of substrate - _w thermal diffusivity of water - _i thermal conductivity of ice - _m thermal conductivity of substrate - _w thermal conductivity of water - wavelength - _av densitiy of air-vapor mixture - _i density of ice - Stefan-Boltzman constant     

Analytical approximation for nonlinear adsorbing solute transport and first-order degradation

Under some constraints, solutes undergoing nonlinear adsorption migrate according to a traveling wave. Analytical traveling wave solutions were used to obtain an approximation for the solute front shape,c(z, t), for the situation of equilibrium nonlinear adsorption and first-order degradation. This approximation describes numerically obtained fronts and breakthrough curves well. It is shown to describe fronts more accurately than a solution based on linearized adsorption. The latter solution accounts neither for the relatively steep downstream solute front nor for the deceleration in time of the nonlinear front.Notation A parameter - c concentration [mol/m³] - c ₀ ^* depth-dependent local maximum concentration [mol/m³] - c; c ₀;c _i concentration difference, feed, and initial resident concentrations, respectively [mol/m³] - D pore scale diffusion/dispersion coefficient [m²/yr] - f adsorption isotherm - f derivative off toc - f second derivative off toc - G ^* parameter - K nonlinear adsorption coefficient [mol/m³)^1–n ] - l column length [m] - L _d dispersivity [m] - m parameter - n Freundlich sorption parameter - P function ofc ₀ ^* - _q change inq [mol/m³] - q adsorbed amount (volumetric basis) [mol/m³] - q derivative ofq toc - R nonlinear retardation factor - retardation factor for concentrationc - R _l linear retardation factor - R(z ^*) depth-dependent average retardation factor, for front at depthz ^* - s adsorbed amount (mass basis) [mol/kg] - t time [years] - u parameter - v flow velocity [m] - z ^* downstream front depth [m] - z depth [m] - transformed coordinate [m] - ^* reference point value of [m] - first-order decay parameter [y^–1] - dry bulk density [kg/m³] - volumetric water fraction - parameter     

Laminar source flow between parallel porous disks with equal suction and injection

An analysis is presented for laminar source flow between parallel stationary porous disks with suction at one of the disks and equal injection at the other. The solution is in the form of an infinite series expansion about the solution at infinite radius, and is valid for all suction and injection rates. Expressions for the velocity, pressure, and shear stress are presented and the effect of the cross flow is discussed.Nomenclature a distance between disks - A, B, ..., J functions of R _w only - F static pressure - p dimensionless static pressure, p(a ²/ ²) - Q volumetric flow rate of the source - r radial coordinate - r dimensionless radial coordinate, r/a - R radial coordinate of a point in the flow region - R dimensionless radial coordinate of a point in the flow region, R - Re source Reynolds number, Q/2a - R _w wall Reynolds number, Va/ - reduced Reynolds number, Re/r ² - critical Reynolds number - velocity component in radial direction - u dimensionless velocity component in radial direction, a/ - average radial velocity, Q/2a - u dimensionless average radial velocity, Re/r - ratio of radial velocity to average radial velocity, u/u - velocity component in axial direction - v dimensionless velocity component in axial direction, v - V magnitude of suction or injection velocity - z axial coordinate - z dimensionless axial coordinate, z a - viscosity - density - kinematic viscosity, / - shear stress at lower disk - shear stress at upper disk - ₀ dimensionless shear stress at lower disk, - ₁ dimensionless shear stress at upper disk, - dimensionless stream function     

Determination of the steady-state shear viscosity from measurements of the apparent viscosity for some common types of viscometers

Considering a number of model fluids, the relation between the (measurable) apparent viscosity _a and the (true) shear viscosity is studied for some commonly used viscometers, like capillary, slit, plate-plate and concentric cylinders (including the influence of the bottom of the cylinder), as well as for one laboratory type of viscometer. As long as is a purely monotonic function, a shift factor < 1 allows one to deduce from _a . Though in general variable, it frequently suffices for practical purposes to use a constant shift factor (the constant being characteristic of the type of viscometer used). This does not apply to dilute solutions or any fluids with two plateau values for . For plastic fluids, it is shown that Casson or Bingham behavior can — if valid at all — only describe the high shear stress limit of _a .     

The shear viscosity dependence on concentration,molecular weight,and shear rate of polystyrene solutions

The solution viscosity of narrow molecular weight distribution polystyrene samples dissolved in toluene and trans-decalin was investigated. The effect of polymer concentration, molecular weight and shear rate on viscosity was determined. The molecular weights lay between 5 10⁴ and 24 10⁶ and the concentrations covered a range of values below and above the critical valuec ^*, at which the macromolecular coils begin to overlap. Flow curves were generated for the solutions studied by plotting log versus log . Different molecular weights were found to have the same viscosity in the non-Newtonian region of the flow curves and follow a straight line with a slope of – 0.83. A plot of log ₀ versus logM _w for 3 wt-% polystyrene in toluene showed a slope of approximately 3.4 in the high molecular weight regime. Increasing the shear rate resulted in a viscosity that was independent of molecular weight. The sloped (log)/d (logM _w ) was found to be zero for molecular weights at which the corresponding viscosities lay on the straight line in the power-law region.On the basis of a relation between _sp and the dimensionless productc · [], simple three-term equations were developed for polystyrene in toluene andt-decalin to correlate the zero-shear viscosity with the concentration and molecular weight. These are valid over a wide concentration range, but they are restricted to molar masses greater than approximately 20000. In the limit of high molecular weights the exponent ofM _w in the dominant term in the equations for both solvents is close to the value 3.4. That is, the correlation between _sp andc · [] results in a sloped(log _sp)/d(logc · []) of approximately 3.4/a at high values ofc · [] wherea is the Mark-Houwink constant. This slope of 3.4/a is also the power ofc in the plot of ₀ versusc at high concentrations. a Mark-Houwink constant - B ₁,B ₂,B _n constants - c concentration (g · cm^–3) - c ^* critical concentration (g · cm^–3) - K, K constants - K _H Huggins constant - M molecular weight - M _c critical molecular weight - M _n number-average molecular weight - M _w weight-average molecular weight - n sloped(log _sp)/d (logc · []) at highc · [] - PS polystyrene - T temperature (K) - shear rate (s^–1) - critical shear rate (s^–1) - viscosity (Pa · s) - ₀ zero-shear viscosity (Pa · s) - _s solvent viscosity (Pa · s) - _sp specific viscosity - [] intrinsic viscosity (cm³ · g^–1) - dynamic viscosity (Pa · s) - | ^*| complex dynamic viscosity (Pa · s) - angular frequency (rad/s) - density of polymer solution (g · cm^–3) - ₁₂ shear stress (Pa) Dedicated to Prof. Dr. J. Schurz on the occasion of his 60^th birthday.Excerpt from the dissertation of Reinhard Kniewske: Bedeutung der molekularen Parameter von Polymeren auf die viskoelastischen Eigenschaften in wäßrigen und nichtwäßrigen Medien, Technische Universität Braunschweig 1983.     

The measurement of the material parameters of viscoelastic fluids using a rotating sphere and a rheogoniometer

Summary In this work, measurement of the flow field around a rotating sphere has been used to obtain the material parameters of a second-order Rivlin-Ericksen fluid. Experiments were carried out with a Laser-Doppler anemometer to obtain the velocity distribution and usingGiesekus' analysis, the material parameters for the second-order fluid were obtained.

---

Zusammenfassung In dieser Untersuchung wird die Ausmessung des Strömungsfeldes um eine rotierende Kugel dazu verwendet, um die Stoffparameter einer Rivlin-Ericksen-Flüssigkeit zweiter Ordnung zu erhalten. Die Experimente zur Bestimmung der Geschwindigkeitsverteilung werden mit einem Laser-Doppler-Anemometer durchgeführt, und zur Auswertung der Parameter der Flüssigkeit zweiter Ordnung wird eine Analyse vonGiesekus benutzt.

---

Notations A ₁,A₂ Rivlin-Ericksen tensor - A ² Parameter used in eq. [12] - a Radius of the sphere - B Parameter used in eq. [12] - I Unit tensor - m ₀(₁ – ₂)/a², parameter used by ref. (8) - N ₁,N₂ First and second normal stress difference - p Isotropic pressure - Radial distance from the centre of the rotating body - S ₁,S₂ Stress tensor - v _r,v,v Velocity components in a spherical coordinate system - ₀,₁,₂ Material parameters used in eq. [2] - Shear rate - _a Apparent voscosity - ₀ Zero-shear viscosity - Angle measured from the axis of rotation - Fluid density - Stream function - Shear stress - Angular velocity With 3 figures     

Maximum density effects on thermal instability of horizontal laminar boundary layers

Linear stability theory is used to investigate the onset of longitudinal vortices in laminar boundary layers along horizontal semi-infinite flat plates heated or cooled isothermally from below by considering the density inversion effect for water using a cubic temperature-density relationship. The analysis employs non-parallel flow model incorporating the variation of the basic flow and temperature fields with the streamwise coordinate as well as the transverse velocity component in the disturbance equations. Numerical results for the critical Grashof number Gr _L ^* =Gr _X ^* /Re _X< ^{Emphasis>/3/2} are presented for thermal conditions corresponding to –0.5 ₁–2.0 and –0.8 ₂1.2.Nomenclature a wavenumber, 2/ - D operator, d/d - F (f–f)/2 - f dimensionless stream function - g gravitational acceleration - G eigenvalue, Gr _L/Re_L - Gr _L Grashof number based on L - Gr _X Grashof number based on X - L characteristic length, (X/U)^1/2 - M number of divisions in y direction - P pressure - Pr Prandtl number, / - p dimensionless pressure, P/( ² /Re _L) - Re _L, Re_X Reynolds numbers, (U L/)=Re _X< ^1/2 and (U), respectively - T temperature - U, V, W velocity components in X, Y, Z directions - u, v, w dimensionless perturbation velocities, (U, V, W)/U - X, Y, Z rectangular coordinates - x, y, z dimensionless coordinates, (X, Y, Z)/L - thermal diffusivity - coefficient of thermal expansion - ₁, ₂ temperature coefficients for density-temperature relationship - similarity variable, Y/L=y - dimensionless temperature disturbance, /T - dimensionless wavelength of vortex rolls, 2/a - ₁, ₂ thermal parameters defined by equation (12) - kinematic viscosity - density - dimensionless basic temperature, (T _b T )/T - –1 - T temperature difference, (T _w–T ) - * critical value or dimensionless disturbance amplitude - prime, disturbance quantity or differentiation with respect to - b basic flow quantity - max value at a density maximum - w value at wall - free stream condition     

Helical flow of molten polymers in a cylindrical annulus

Summary The paper is concerned with an analytical investigation of helical flow of a non-Newtonian fluid through an annulus with a rotating inner cylinder. The shear dependence of viscosity is described by a power law and the temperature dependence by an exponential function.Velocity and temperature profiles, energy input and shear along the stream lines, pressure drop, and torque are presented for the range of input parameters encountered in polymer extrusion.The results of the study can be applied to a mixing element in a screw extruder and for a device to control extrudate temperature and output.Nomenclature a thermal diffusivity [m²/s] - b temperature coefficient [K^–1], see eq. [4] - c heat capacity [J/kg K] - h slot width [m] - I ₁,I ₂,I ₃ invariants of the rate of deformation tensor, see eq. [5] - k thermal conductivity [J/m s K] - l, L = 1/h length of the slot - l _T ,l _K thermal and kinematic entrance length - m power law exponent, see eq. [3] - M torque [m N] - p pressure [N/m²] - P dimensionless pressure gradient, see eq. [24] - P _R,P _RZ dimensionless components of the shear stress tensor, see eq. [25] and eq. [26] - r, R = r/r _wa radial coordinate - r _wa, r_wi outer and inner radius of annulus [m] - t time [s]; dwell time in the annulus - T temperature [K] - v , v_r, V_z velocity components [m/s] - v ₀ angular velocity at inner wall [m/s] - average velocity inz-direction [m/s] - V , V_R, V_Z dimensionless velocity components,v /v₀, v_r/v₀, v_z/v₀ - V _z velocity ratio, helical parameter - Y coordinate inr-direction, see eq. [20] - z, Z = z/h Pe axial coordinate - deformation - rate of deformation tensor [s^–1] - apparent viscosity [N s/m²], see eq. [3] - dimensionless temperature,b (T – T ₀) - azimuth coordinate - ratio of radii,r _wi/r_wa - density [kg/m³] - , kl shear stress tensor [N/m²] - fluidity [m²w/Nw s], see eq. [4] - Gf Griffith number, see eq. [12] - Pe Péclet number, see eq. [13] - Re Reynolds number, - 0 initial state, reference state - equilibrium state - e entrance - wi, wa at surface of inner or outer wall - r, R, z, Z, coordinates - i, j radial and axial position of nodal point in the grid - k, l tensor components Presented at Euromech 37, Napoli 6. 20–23. 1972.With 15 figuresDedicated to Prof. Dr.-Ing. G. Schenkel on his 60th birthday     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

A three-parameter model describing the experimental relation between shear stress and shear rate for laminar flow of aqueous polymer solutions

Summary A three-parameter model is introduced to describe the shear rate — shear stress relation for dilute aqueous solutions of polyacrylamide (Separan AP-30) or polyethylenoxide (Polyox WSR-301) in the concentration range 50 wppm – 10,000 wppm. Solutions of both polymers show for a similar rheological behaviour. This behaviour can be described by an equation having three parameters i.e. zero-shear viscosity ₀, infinite-shear viscosity , and yield stress ₀, each depending on the polymer concentration. A good agreement is found between the values calculated with this three-parameter model and the experimental results obtained with a cone-and-plate rheogoniometer and those determined with a capillary-tube rheometer.

---

Zusammenfassung Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit von strukturviskosen Flüssigkeiten wird durch ein Modell mit drei Parametern beschrieben. Mit verdünnten wäßrigen Polyacrylamid-(Separan AP-30) sowie Polyäthylenoxidlösungen (Polyox WSR-301) wird das Modell experimentell geprüft. Beide Polymerlösungen zeigen im untersuchten Schergeschwindigkeitsbereich von ein ähnliches rheologisches Verhalten. Dieses Verhalten kann mit drei konzentrationsabhängigen Größen, nämlich einer Null-Viskosität ₀, einer Grenz-Viskosität und einer Fließgrenze ₀ beschrieben werden. Die Ergebnisse von Experimenten mit einem Kegel-Platte-Rheogoniometer sowie einem Kapillarviskosimeter sind in guter Übereinstimmung mit den Werten, die mit dem Drei-Parameter-Modell berechnet worden sind.

---

a Pa^–1 physical quantity defined by:a = {1 – ( / ₀)}/ ₀ - c l concentration (wppm) - D m capillary diameter - L m length of capillary tube - P Pa pressure drop - R m radius of capillary tube - u m s^–1 average velocity - v _r m s^–1 local axial velocity at a distancer from the axis of the tube - shear rate (–dv _r /dr) - local shear rate in capillary flow - s^–1 wall shear rate in capillary flow - Pa s dynamic viscosity - _a Pa s apparent viscosity defined by eq. [2] - ( _a ) Pa s apparent viscosity in capillary tube at a distanceR from the axis - ₀ Pa s zero-shear viscosity defined by eq. [4] - Pa s infinite-shear viscosity defined by eq. [5] - l ratior/R - kg m^– density - Pa shear stress - ₀ Pa yield stress - _r Pa local shear stress in capillary flow - _R Pa wall shear stress in capillary flow _R = (PR/2L) - _v m³ s^–1 volume rate of flow With 8 figures and 1 table     

A note on three-dimensional boundary layer growth

Summary This note is an extension of the work of Görtler²) on two-dimensional boundary layer growth to the three-dimensional case. The solutions of three-dimensional boundary layer equations are obtained by considering the potential flow of the body to be governed by the functions At U ₀(, ) and At U ₀(, ) where is any positive number.     

The chaotic response of a fluttering panel: The influence of maneuvering

The influence of maneuvering on the chaotic response of a fluttering buckled plate on an aircraft has been studied. The governing equations, derived using Lagrangian mechanics, include geometric non-linearities associated with the occurrence of tensile stresses, as well as coupling between the angular velocity of the maneuver and the elastic degrees of freedom. Numerical simulation for periodic and chaotic responses are conducted in order to analyze the influence of the pull-up maneuver on the dynamic behavior of the panel. Long-time histories phase-plane plots, and power spectra of the responses are presented. As the maneuver (load factor) increases, the system exhibits complicated dynamic behavior including a direct and inverse cascade of subharmonic bifurcations, intermittency, and chaos. Beside these classical routes of transition from a periodic state to chaos, our calculations suggest amplitude modulation as a possible new mode of transition to chaos. Consequently this research contributes to the understanding of the mechanisms through which the transition between periodic and strange attractors occurs in, dissipative mechanical systems. In the case of a prescribed time dependent maneuver, a remarkable transition between the different types of limit cycles is presented.Nomenclature a plate length - a _r u _r /h - D plate bending stiffness - E modulus of elasticity - g acceleration due to gravity - h plate thickness - j₁,j₂,j₃ base vectors of the body frame of reference - K spring constant - M Mach number - n 1 + ₀/g - N ₁ applied in-plane force - p–p aerodynamic pressure - P pa ⁴/Dh - q ₀/2 - Q _r generalized Lagrangian forces - R rotation matrix - R ₄ N, a ²/D - t time - kinetic energy - u plate deflection - u displacement of the structure - u _r modal amplitude - v₀ velocity - x coordinates in the inertial frame of reference - z coordinates in the body frame of reference - Ka/(Ka+Eh) - - elastic energy - 2qa ³/D - a/_mh - Poisson's ratio - material coordinates - air density - _m plate density - - _r prescribed functions - _r sin(r z/a) - angular velocity - a/v₀ - skew-symmetric matrix form of the angular velocity     

Uniqueness and the maximum principle for quasilinear elliptic equations with quadratic growth conditions

In this paper, we show that the maximum principle holds for quasilinear elliptic equations with quadratic growth under general structure conditions.Two typical particular cases of our results are the following. On one hand, we prove that the equation (1) {ie77-01} where {ie77-02} and {ie77-03} satisfies the maximum principle for solutions in H ¹()L(), i.e., that two solutions u ₁, u ₂H¹() L() of (1) such that u ₁u₂ on , satisfy u ₁u₂ in . This implies in particular the uniqueness of the solution of (1) in H ₀ ¹ ()L().On the other hand, we prove that the equation (2) {ie77-04} where fH^–1() and g(u)>0, g(0)=0, satisfies the maximum principle for solutions uH¹() such that g(u)¦Du|{²L¹(). Again this implies the uniqueness of the solution of (2) in the class uH ₀ ¹ () with g(u)¦Du|{²L¹().In both cases, the method of proof consists in making a certain change of function u=(v) in equation (1) or (2), and in proving that the transformed equation, which is of the form (3) {ie77-05}satisfies a certain structure condition, which using ((v₁ -v ₂)⁺)ⁿ for some n>0 as a test function, allows us to prove the maximum principle.     

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     

An analytical solution for the velocity and shear rate distribution of non-ideal Bingham fluids in concentric cylinder viscometers

An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/_y - c = 2µ - d = ₀ ²–2c_y - f() = (–₀)²+2c(–_y) - r Radius - r _b Bob radius - r _c Cup radius - r _y Yield radius - r ₀ Stationary surface radius - r Rotating Stationary radius - Y ₀ Shear rate parameter = /µ Greek letters Shear rate - = (r _y /r _b )²– 1 - µ Bingham viscosity - µ₀ Initial differential viscosity - µ µ₀-µ - Rotation speed - Angular velocity - Shear stress - _b Bob shear stress - _B Bingham stress - _y (True) yield stress - ₀ Stress parameter = _B +µ Y ₀ - _B - _y      

The theory of a vibrating-rod densimeter

The paper presents a theory of a device for the accurate determination of the density of fluids over a wide range of thermodynamic states. The instrument is based upon the measurement of the characteristics of the resonance of a circular section tube, or rod, performing steady, transverse oscillations in the fluid. The theory developed accounts for the fluid motion external to the rod as well as the mechanical motion of the rod and is valid over a defined range of conditions. A complete set of working equations and corrections is obtained for the instrument which, together with the limits of the validity of the theory, prescribe the parameters of a practical design capable of high accuracy.Nomenclature A, B, C, D constants in equation (60) - A _j , B _j constants in equation (18) - a _j ⁺ , a _j ^– wavenumbers given by equation (19) - C _f drag coefficient defined in equation (64) - C _f ^/0 , C _f ^/1 components of C _f in series expansion in powers of - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - force per unit length - F _j ⁺ , F _j ^– constants in equation (24) - f, g functions of defined in equations (56) - G modulus of rigidity - I second moment of area - K constant in equation (90) - k, k constants defined in equations (9) - L half-length of oscillator - Ma Mach number - m _a mass per unit length of fluid a - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equation (17) - P power (energy per cycle) - P _a , P _b power in fluids a and b - p pressure - R radius of rod or outer radius of tube - R _c radius of container - R _i inner radius of tube - r radial coordinate - T tension - T _visc temperature rise due to heat generation by viscous dissipation - t time - v _r , v radial and angular velocity components - y lateral displacement - z axial coordinate - dimensionless tension - _a dimensionless mass of fluid a - _b dimensionless added mass of fluid b - _b dimensionless drag of fluid b - dimensionless parameter associated with - ₀ dimensionless coefficient of internal damping - dimensionless half-width of resonance curve - dimensionless frequency difference defined in equation (87) - spatial resolution of amplitude - R, , , _s , increments in R, , , _s , - dimensionless amplitude of oscillation - dimensionless axial coordinate - ratio of to - _a , _b ratios of to for fluids a and b - angular coordinate - parameter arising from distortion of initially plane cross-sections - _f thermal conductivity of fluid - dimensionless parameter associated with - viscosity of fluid - _a , _b viscosity of fluids a and b - dimensionless displacement - _j jth component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - density of fluid calculated on assumption that * - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - _{rr rr}, _r radial normal and shear stress components - spatial component of defined in equation (13) - _j jth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - phase angle - _r phase difference - _ra , _rb phase difference for fluids a and b - streamfunction - _j jth component defined in equation (22) - dimensionless frequency (based on ) - _a , _b dimensionless frequency in fluids a and b - _s dimensionless frequency (based on _s ) - angular frequency - ₀ resonant frequency in absence of fluid and internal damping - _r resonant frequency in absence of internal fluid - _ra , _rb resonant frequencies in fluids a and b - dimensionless frequency - dimensionless frequency when _a vanishes - dimensionless frequencies when _a vanishes in fluids a and b - dimensionless resonant frequency when _a , _b, _b and ₀ vanish - dimensionless resonant frequency when _a , _b and _b vanish - dimensionless resonant frequency when _b and _b vanish - dimensionless frequencies at which amplitude is half that at resonance     

The rheology of pigment dispersions

Summary A kinetic model is developed to relate the measured shear stress in a dispersion with the rate of deformation, and with the level of structure caused by the competing effects of flocculation and deflocculation.The model parameters are determined from experimental data obtained from equilibrium and transient oscillatory shear, using dispersions of a pigment in three different oil-based media. It is found that the model can successfully describe the flow behaviour of the dispersions under all three types of deformation, and account for different concentrations and temperatures.

---

Zusammenfassung Es wird ein kinetisches Modell entwickelt, das die in einer Dispersion gemessene Schubspannung mit der Deformationsgeschwindigkeit in Beziehung setzt unter Berücksichtigung der im Wettbewerb stehenden Flokkulations-und Deflokkulationseffekte.Die Modellparameter werden mit Hilfe experimenteller Daten bestimmt, die zu stationären und oszillatorischen Scherströmungen unter Einbeziehung wechselnder Beanspruchungsarten erhalten wurden. Dabei wurden Pigment-Dispersionen in drei verschiedenen Medien auf Öl-Basis verwendet. Man findet, daß das Modell das Fließverhalten der Dispersionen unter allen betrachteten Deformationstypen sowie bei den verschiedenen angewandten Konzentrationen und Temperaturen erfolgreich zu beschreiben vermag.

---

a, b experimental constants - c dispersion concentration by weight - d ₃ mean volume to surface diameter of a floc. - f (·,·) function defined by eq. [20] - g(·) a function of volume fraction - k Boltzmann's constant - n _s number of floccules containings flocs, per unit volume - ratio of the number of floccules containings flocs, per unit volume, to the total number of flocs per unit volume - t present time - B, D derived constants - E ₀,E ₁ activation energies for viscous flow of a dispersion at low and high shear rates - E _m activation energy for viscous flow of the medium - E _f ,E _d activation energies for flocculation and deflocculation of a dispersion - G shear rate - K ₁,K ₂ model constants - M parameter related to the total number of flow units - N number of flocs per unit volume - R, R ₀ the ratio of the rate function for deflocculation to the rate function for flocculation, and its value in the absence of shear - T absolute temperature - ₀, ₀ constants - (·,·),(·) rate functions for flocculation and deflocculation - , , v model parameters - , ₀, ₁ viscosity of a dispersion, and its value at low and high shear rates - _r the viscosity of a dispersion of floccules each containingr flocs - dynamic viscosity of a dispersion - µ, µ ₀ viscosity of a medium at temperatureT, and in the limit of high temperature - _p , _f , _F volume fraction of pigment, of flocs, of floccules - measured shear stress - non-dimensional time - a characteristic time for flocculation With 8 figures and 7 tables     

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      

Heat transfer and equilibrium temperature of a sphere in a supersonic rarefied gas flow

Some results are presented of experimental studies of the equilibrium temperature and heat transfer of a sphere in a supersonic rarefied air flow.The notations D sphere diameter - u, , T,,l, freestream parameters (u is velocity, density, T the thermodynamic temperature,l the molecular mean free path, the viscosity coefficient, the thermal conductivity) - T₀ temperature of the adiabatically stagnated stream - T_e mean equilibrium temperature of the sphere - T_w surface temperature of the cold sphere (T_we) - mean heat transfer coefficient - _e air thermal conductivity at the temperature T_e - P Prandtl number - M Mach number