Peristaltic pumping of a second-order fluid in a planar channel

The peristaltic motion of a non-Newtonian fluid represented by the constitutive equation for a second-order fluid was studied for the case of a planar channel with harmonically undulating extensible walls. A perturbation series for the parameter ( half-width of channel/wave length) obtained explicit terms of 0(²), 0(²Re²) and 0(₁Re²) respectively representing curvature, inertia and the non-Newtonian character of the fluid. Numerical computations were performed and compared to the perturbation analysis in order to determine the range of validity of the terms.Presented at the second conference Recent Developments in Structured Continua, May 23–25, 1990, in Sherbrooke, Québec, Canada     

On periodic motions in magnetohydrodynamics

Summary A uniqueness theorem and an existence theorem for motions asymptotically stable in the mean and periodically depending on time, are demonstrated in magnetohydrodynamics, on the hypothesis that a motion sufficiently regular and asymptotically stable in the mean exists. In particular, similar theorems for steady motions are deduced. The procedure used is also a simple method for constructing periodic and steady motions.

---

Sommario Si dimostra un teorema di unicità ed un teorema di esistenza per i moti magnetoidrodinamici periodicî nel tempo ed asintoticamente stabili in media, subordinatamente all'ipotesi che esista un moto asintoticamente stabile in media sufficientemente regolare, e si deducono, in particolare, analoghi teoremi per i moti stazionari. Il procedimento adoperato costituisce anche un semplice metodo per la effettiva costruzione dei moti periodici e dei moti stazionari.

---

---

This work was done in the sphere of activity of the C.N.R. groups for mathematical research and, in Italian, has already been presented for publication in Rend. Acc. Sc. Fis. Mat., Napoli.     

Torsion of a hemisphere embedded in an elastic half-space

This paper deals with the stress distribution in a homogeneous isotropic elastic hemisphere embedded in a semi-infinite homogeneous isotropic elastic medium when a rigid circular disc is clamped to the plane face of the hemisphere and the stresses are caused by the rotation of the disc through an angle . The problem is reduced to the solution of a Fredholm integral equation of the second kind in the auxiliary function (t). An analytical expression for the torque T required to rotate the die through an angle is obtained in terms of (t). The Fredholm integral equation is solved numerically, and the numerical values of T are graphed.This work has been supported by the National Research Council of Canada through NRC-Research Grant No. A4177.     

Various reciprocal theorems and variational principles in the theories of nonlocal micropolar linear elastic mediums

In the first part of our paper, we have extended the concepts of the classical convolution and the convolution scalar product given by I. Hlaváck and presented the concepts of the convolution vector and the convolution vector scalar product, which enable us to extend the initial value as well as the initial-boundary value problems for the equation with the operator coefficients to those for the system of equations with the operator coefficients.In the second part of this paper, based on the concepts of the convolution vector and the convolution vector scalar product, two fundamental types of reciprocal theorems of the non-local micropolar linear elastodynamics for inhomogeneous and anisotropic solids are derived.In the third part of this paper, based on the concepts and results in the first and second parts as well as the Lagrange multiplies method which is presented by W. Z. Chien, four main types of variational principles are given for the nonlocal micropolar linear elastodynamics for inhomogeneous and anisotropic solids. These are the counterparts of the variational principles of Hu-Washizu type, Hellinger-Reissner type and Gurtin type in classical elasticity as well as Hlaváck type and Iesan type in local micropolar and nonlocal elasticity. Finally, we have proved the equivalence of the last two main variational principles which are given in this paper.     

A particle transport problem with a quadratic nonlinearity and general boundary conditions: The stationary solution

Summary We study a stationary, nonlinear, particle transport problem in slab geometry with general boundary conditions. The existence and uniqueness of the solution is proved by means of fixed point techniques, provided that the source term is sufficiently small.

---

Sommario Si studia un problema stazionario nonlineare di particelle in geometria piana con condizioni al contorno generali. L'esistenza e unicitá della soluzione è dimostrata con tecniche di punto fisso purchè il termine di sorgente sia sufficientemente piccolo.

     

Theoretical investigation of the effect of an internal heat source on turbulent liquid metal heat transfer in a channel with uniform heat flux

Summary The effect of an internal heat source on the heat transfer characteristics for turbulent liquid metal flow between parallel plates is studied analytically. The analysis is carried out for the conditions of uniform internal heat generation, uniform wall heat flux, and fully established temperature and velocity profiles. Consideration is given both to the uniform or slug flow approximation and the power law approximation for the turbulent velocity profile. Allowance is made for turbulent eddying within the liquid metal through the use of an idealized eddy diffusivity function. It is found that the Nusselt number is unaffected by the heat source strength when the velocity profile is assumed to be uniform over the channel cross section. In the case of a 1/7-power velocity expression, the Nusselt numbers are lower than those in the absence of internal heat generation, and decrease with diminishing eddy conduction. Nusselt numbers, in the absence of an internal heat source, are compared with existing calculations, and indications are that the present results are adequate for preliminary design purposes.Nomenclature A hydrodynamic parameter - a half height of channel - a ₁ a constant, 1+0.01 Pr Re ^0.9 - a ₂ a constant, 0.01 Pr Re ^0.9 - C _p specific heat at constant pressure - D _h hydraulic diameter of channel, 4a - h heat transfer coefficient, q _w/(t _w–t _b) - I ₁ integral defined by (17) - I ₂ integral defined by (18) - k diffusivity parameter, (1+0.01 Pr Re ^0.9)^1/2 - m exponent in power velocity expression - Nu Nusselt number, hD _h/ - Nu ₀ Nusselt number in absence of internal heat generation - Pr Prandtl number, / - Q heat generation rate per volume - q _w wall heat flux - Re Reynolds number for channel, 2/ - s ratio of heat generation rate to wall heat flux, Qa/q _w - T dimensionless temperature, (t _w–t)/(t _w–t _b) - t fluid temperature, t _w wall temperature, t _b fluid bulk temperature - u fluid velocity in x direction, , fluid mean velocity - x longitudinal coordinate measured from channel entrance - x ⁺ dimensionless longitudinal coordinate, 2(x/a)/Pr Re - y transverse coordinate measured from channel centerline - z transverse coordinate measured from channel wall, a–y - molecular diffusivity of heat, /C _p - dummy variable of integration - dummy variable of integration - _H eddy diffusivity of heat - _M eddy diffusivity of momentum - dummy variable of integration - fluid thermal conductivity - _T dimensionless diffusivity, Pr ( _H/) - fluid kinematic viscosity - dummy variable of integration - fluid density - dummy variable of integration - ratio of eddy diffusivity for heat transfer to that for momentum transfer, _H/ _M - average value of - dimensionless velocity distribution, u/     

Natural convection heat transfer through an enclosed horizontal layer of supercritical carbon dioxide

In natural convection heat transfer through a thin horizontal layer of carbon dioxide, maxima in the equivalent thermal conductivities are obtained in the vicinity of the respective pseudocritical temperatures at pressures of 75.8, 89.6 and 103.4 bar. The maxima are the more pronounced, the closer the critical point is approached.Comparison of experimental results with Nusselt equations shows good agreement except for the immediate vicinity of the pseudocritical temperature.In visual observations a distinct change in flow structure appears in the immediate vicinity of the pseudocritical temperature. A steady state polygon pattern and a boiling-like action could not be observed in this geometry.

---

Zusammenfassung Beim Wärmetransport durch freie Konvektion in einer dünnen waagerechten Schicht von Kohlendioxid ergaben sich Maxima der scheinbaren Wärmeleitfähigkeit in der Nähe der pseudokritischen Temperaturen bei Drükken von 75,8, 89,6 und 103,4 bar. Die Maxima sind um so ausgeprägter, je mehr man sich dem kritischen Punkt nähert.Ein Vergleich der Versuchsergebnisse mit Nusseltbeziehungen ergibt gute Übereinstimmung außer in unmittelbarer Umgebung der pseudokritischen Temperatur. Direkte Beobachtungen der Konvektionsmuster zeigen in unmittelbarer Umgebung der pseudokritischen Temperatur eine deutliche Strukturänderung. Ein stationäres Zellmuster und siedeähnliche Vorgänge konnten in dieser Anordnung nicht beobachtet werden.

---

Nomenclature A area of the heating or cooling plate - C constant in the correlation - g acceleration of gravity - h heat transfer coefficient - k thermal conductivity of fluid in the gap - k _e equivalent thermal conductivity - m, n exponents of dimensionless numbers - q heat flux - T _C,PC absolute temperature; critical C, pseudocritical PC - Gr Grashof numberg ( _h– _c) ³/ ² - Nu Nusselt numberh/k - Pr Prandtl number/ - thermal diffusivity - coefficient of volume expansion - width of gap - _c,h temperature of cooling (c)-, heating (h)-plate - _m arithmetic mean temperature ( _c+ _h)/2 - kinematic viscosity - _c,h fluid density at the temperature of the cooling (c)- or heating (h)-plate - heat flow rate through the gap     

Nonstationary filtration in partially saturated porous media

From the mathematical formulation of a one-dimensional flow through a partially saturated porous medium, we arrive at a nonlinear free boundary problem, the boundary being between the saturated and the unsaturated regions in the medium. In particular we obtain an equation which is parabolic in the unsaturated part of the domain and elliptic in the saturated part.Existence, uniqueness, a maximum principle and regularity properties are proved for weak solutions of a Cauchy-Dirichlet problem in the cylinder {(x,t): 0x1, t0} and the nature, in particular the regularity, of the free boundary is discussed.Finally, it is shown that solutions of a large class of Cauchy-Dirichlet problems converge towards a stationary solution as t and estimates are given for the rate of convergence.     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

An approximate theoretical analysis and experimental verification of turbulent entrance region flow of drag reducing fluids

Summary Entry lengths for pipe flows of moderately drag reducing fluids are determined using momentum integral technique. It is shown theoretically that the entry lengths for drag reducing fluids could be significantly larger than the Newtonian fluids flowing turbulently under otherwise identical conditions. The experimental data from the literature bear out the theoretical calculations.

---

Zusammenfassung Mit Hilfe der Impuls-Methode wird die Einlauflänge in einer Rohrströmung für Flüssigkeiten mit mäßig starker Widerstandsverminderung berechnet. Es wird vorausgesagt, daß die Einlauflänge für derartige Flüssigkeiten erheblich größer sein kann als für newtonsche Flüssigkeiten unter sonst identischen Bedingungen. Aus der Literatur entnommene experimentelle Daten bestätigen diese theoretischen Berechnungen.

---

Nomenclature A ₁ Coefficient in eq. [7] - A Slope of logarithmic velocity profile - a Exponent in eq. [10] - B Intercept function for logarithmic velocity profile - De Deborah number, - f Friction factor - F Function, eq. [30] - G Function given in eq. [11] - Static pressure, dynes/cm² - q Index of power law velocity profile - R Pipe radius, cm - r Radial distance, cm - R Core radius, cm - Re Reynolds number - Axial velocity, cm/s - u _c Core velocity, cm/s - u ⁺ Dimensionless velocity, eq. [5] - u ^* Friction velocity, , cm/s - Radial velocity, cm/s - V Average velocity, cm/s - x Axial distance, cm - x _e Entry length, cm - y Distance from the wall, cm - y ⁺ Dimensionless distance, eq. [5] - y _I ⁺ Dimensionless viscous sublayer thickness - coefficient in eq. [17] - exponent of Reynolds number in eq. [17] - shear rate, s^–1 - boundary layer thickness, cm - _fl fluid relaxation time, s - µ fluid viscosity, gm/cm s - v kinematic viscosity, cm²/s - _l laminar sublayer thickness, dimensionless - fluid density, gm/cm³ - shear stress, dynes/cm² - _w shear stress at the wall, dynes/cm² - ₁, ₂, ₃, ₄ functions in eq. [27] - ~ time averaged quantities - — dimensionless quantity With 3 figures and 1 table     

A possible mechanism for stabilizing elongational flow in certain polymeric liquids at constant temperature and composition

Summary A new calculation of the stability of elongational flow of a liquid filament is based on an idealized treatment of the effects of non-uniformity of filament diameter. Two uniform cylindrical filaments of different initial cross-sectional areas are connected by a fictitious deformable clamp which transmits tensile force without interfering with the assumed homogeneous elongational flow of each filament. A uniform tensile force is suddenly applied and kept constant. For a Newtonian liquid, the area ratio rapidly tends to infinity. For a rubberlike liquid, the area ratio at first increases but then becomes constant; this suggests that finite deformation of a Gaussian molecular network may provide a possible mechanism for spinnability at constant temperature and composition.

---

Zusammenfassung Eine neuartige Berechnung der Stabilität des Dehnungs-Fließverhaltens wird auf einer idealisierten Behandlung des Einflusses der Ungleichförmigkeit des Fadendurchmessers aufgebaut. Das Modell besteht aus zwei gleichförmigen zylindrischen Fäden von verschiedenem Ausgangsquerschnitt mit einer fiktiven Verbindung, die Zugkräfte überträgt ohne die angenommenen homogenen Dehnungs-Fließzustände der beiden Fäden zu beeinträchtigen. Eine konstante Zugkraft wird gleichförmig und sprungartig aufgebracht. Für eine Newtonsche-Flüssigkeit geht das Verhältnis der Fadenquerschnitte sehr schnell gegen unendlich. Für eine gummiartige Flüssigkeit nimmt das Verhältnis zuerst zu, bleibt dann aber konstant. Das erlaubt den Schluß, daß endliche Deformation eines Gauss schen molekularen Netzwerks einen möglichen Mechanismus für die Spinnbarkeit bei konstanter Temperatur und Zusammensetzung zu liefern vermag.

     

Laminar natural convection over a slender vertical frustum of a cone

The problem of laminar, natural convection flow over a slender frustum of a cone is treated in this paper. The governing differential equations are solved by a combination of quasi-linearization and finite-difference methods. Numerical solutions are obtained for Pr=0.7 and for a range of values of the transverse curvature parameter. It is shown that the effect of transverse curvature is of great significance in such flows.

---

Laminare natürliche Konvektion an einem dünnen, senkrechten Kegelstumpf

Zusammenfassung In diesem Bericht ist das Problem der laminaren natürlichen Konvektionsströmung an einem dünnen Kegelstumpf behandelt. Die maßgebliche Differentialgleichung ist durch eine Verbindung von Quasilinearisation und Differenzenverfahren gelöst. Eine numerische Lösung für Pr=0.7 wird für verschiedene Werte eines Krümmungsparameters angegeben. Es ist gezeigt, daß in solchen Strömungen dieser Krümmungsparameter eine große Bedeutung besitzt.

---

Nomenclature f dependent variable, defined in Eq. (7) - g dependent variable, defined in Eq. (7) - g_e gravitational acceleration - h heat transfer coefficient, or -grid - k heat conductivity, or -grid - L characteristic length - Nu Nusselt number - Pr Prandtl number - r radial distance from the axis of the cone - R transverse curvature effect ratio, defined in Eq. (23) - Re Reynold number - T temperature - u, v velocity components in the x- and y-directions, respectively - x, y rectangular coordinates Greek letters dimensionless temperature, definedinEq. (4) - bulk modulus - cone angle - dynamic viscosity - stream function - , independent variable, defined in (7) - transverse curvature parameter     

Spectral and correlation characteristics of the turbulent boundary layer on an extended flexible cylinder

In marine geophysical seismological prospecting extensive use is made of towed receiving systems consisting of extended flexible cylinders containing acoustic sensors over which the water flows in the longitudinal direction. The boundary layer pressure fluctuations on the cylinder surface are picked up by the sensors as hydrodynamic noise. This paper is concerned with the study of the energy and spacetime characteristics of the pressure fluctuations in the turbulent boundary layer on an extended flexible cylinder in a longitudinal flow. The pressure fluctuations on the cylinder surface have been investigated experimentally for Re_X=(2–4)·10⁷, a dimensionless diameter of the pressure fluctuation sensors d _p ⁺ =d_pu^*/=500, where d_p is the sensor diameter, u^* the dynamic viscosity, and the viscosity coefficient, and frequencies 0.02311.259 (=^*/U). The spectral and correlation characteristics of the pressure fluctuations on the surface of the flexible cylinder are found to differ from the corresponding characteristics for a rigid cylinder [1–4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i aza, No, 5, pp. 49–54, September–October, 1989.     

An experimental study on the two-dimensional opposed wall jet in a turbulent boundary layer: Change in the flow pattern with velocity ratio

Results of the measurement of flow properties in a two-dimensional turbulent wall jet which is injected into the turbulent boundary layer in the direction opposite to that of the boundary layer flow are presented by varying the ratio of the jet issuing velocity to the mainstream velocity . This flow includes the flow separation and the recirculating flow, and so it requires the magnitude and direction of instantaneous velocity be measured. In the present work, a tandem hot-wire probe is manufactured and its characteristics are examined experimentally. With the use of this probe the change in the penetration length of the jet with the velocity ratio is investigated. It is clarified that two regimes of flow patterns exist: in the low velocity ratio the penetration length of the jet changes little with , and in the high velocity ratio it goes far from the nozzle with increasing . Streamlines, turbulence intensity contours and static pressure contours are presented in the two typical velocity ratios corresponding to each of two flow patterns, and they are compared.List of symbols b ₀ nozzle width (Fig. 1) - , e mean and fluctuating output voltages of hot-wire anemometer - p, p mean static pressure, p = p – p _o - p ₀ static pressure in undisturbed mainstream - p _w , p _w mean wall pressure, p _w = p _w – p _o - U ₀ mainstream velocity - U _j jet velocity at the nozzle exit - , u mean and fluctuating velocity components in x-direction - u _e effective cooling velocity - x distance along the wall from nozzle exit - x _pmax, x _pmin positions where the wall pressure has maximum and minimum values respectively - x _s penetration length of jet - y distance from the wall - forward flow fraction - ₁, ₂ yaw and pitch angles of flow direction (Fig. 4) - velocity ratio, = U _j /U _o - density of the fluid - nondimensional stream function The authors wish to express their appreciation to Prof. Toshio Tanaka of Gifu University for his advice in the course of the experiment. They also would like to thank the Research Foundation for the Electrotechnology of Chubu which partly supported this work.     

Base-pressure pulsations behind a cylinder and disk in a subsonic stream flow

Base-pressure fluctuations behind a long cylinder (l/d 5–10) and the disk (l/ d 0.0) is investigated experimentally in this paper. The spectral and correlation characteristics of the base-pressure fluctuations behind axisymmetric bodies at a Mach number M 1.0 are generalized on the basis of the data obtained and the results of other authors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 181–183, January–February, 1977.     

Impulsive motion of a semi-infinite flat plate in a viscous fluid

The unsteady laminar boundary layer flow is investigated for a semi-infinite flat plate subjected to impulsive motion. An approximate solution is obtained by utilizing Meksyn's method. These results vary smoothly from Rayleigh's unsteady solution to the steady state solution of Blasius. Results are compared to those of Lam and Crocco.Nomenclature A expansion coefficient, see eq. (13) - a expansion coefficient, see eq. (10) - B expansion coefficient, see eq. (14) - b expansion coefficient, see eq. (12) - G function defined by eq. (6) - U free stream velocity - u velocity in x direction - v velocity in y direction - x coordinate along plate - y coordinate normal to plate Greek symbols (l, ) incomplete gamma function - function defined by eq. (15) - y(U/x) ^1/2 - kinematic viscosity - x/Ut - (Uvx)^1/2 f(, )     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

Some observations of the permanent set of cross-linked natural rubber samples after heating in a state of pure shear

Experimental evidence is presented to show that when cross-linked natural rubber samples are heated in air whilst subject to a pure shear strain and subsequently freed of all forces the new force-free or set state is not related to the original force-free state by a pure shear. During a pure shear linear elements along one of the principal directions of strain remain unchanged in length; it is shown that these elements undergo an increase in length as a result of the recovery after heating. The results are compared with the predictions of the two-network theory of permanent set propounded by previous workers.

---

Zusammenfassung Vernetzte Naturkautschuk-Proben, in Luft erhitzt, zeigen bei gleicher Beanspruchung mit reiner Scherung und anschließender Entlastung einen neuen spannungsfreien (oder gesetzten) Zustand, der mit dem ursprünglichen spannungsfreien Zustand nicht durch eine reine Scherung in Beziehung steht. Durch eine solche bleiben lineare Elemente entlang einer der Hauptachsen der Spannung in ihrer Länge unverändert; es. wird gezeigt, daß diese Elemente an Länge bei der Erholung nach dem Aufheizen zunehmen. Die Resultate werden mit den Voraussagen der Zweinetzwerk-Theorie des permanent set der Autoren verglichen.

     

Effect of suspended particles on the stability of plane poiseuille flow

The behavior of the neutral stability curves is investigated for various values of the particle relaxation time and mass concentration 0 100 and 0 f 0.1. It is shown that as increases from zero the flow is at first destabilized and then at >6 becomes stable, while at >40 the stabilizing effect of the dispersed phase grows weaker. It is found that there is a certain interval 10< <40 on which the flow is most stable.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 46–53, January–February, 1986.