The motion of a detaching elastic body

Conclusions The qualitative behavior of the displacement (t) and the radius R(t) during the different phases of the motion is illustrated in the diagram of Fig. 6.1.After the first impact at t = 0 the displacement (t) varies according to (5.2). If the first maximum of (t) is higher than ₁ then at time t ₁ the graph of (t) intersects the straight line = _cand detachment first occurs. In the second phase the dependance of on t is expressed by (5.6). The detachment will end at the instant t ₂ when vanishes.The radius R remains equal to R ₀ until (t) reaches the critical value ₁ = _c that is at t = t ₁. After t ₁, R(t) will decrease according to (4.4) up to its final value ₂.A rather unexpected property of the solution is that the greatest elongation is finite for every non-vanishing value of the ratio .To Jerry Ericksen for his 60^th birthday     

On the stability of periodic motions of an unbalanced rigid rotor on lubricated journal bearings

A theoretical investigation is carried out on the orbital motions of a symmetrical, unbalanced, rigid rotor subjected to a constant vertical load and supported on two lubricated journal bearings. In order to determine the fluid film forces, the short bearing theory is adopted.A method is illustrated that makes it possible to determine the analytical equation of the orbit as an approximated solution of the system of non-linear differential equations of motion of the journal axis. A procedure is also described for evaluating the stability of the solution found. Diagrams of the curves delimiting, in the working plane of the rotor -m, the areas of stability of the various periodic solutions determined are provided.Finally, the results obtained are compared and combined with those provided by a direct integration of the motion equation made using the Runge-Kutta method.Nomenclature C radial clearance - D = 2R bearing diameter - E mass unbalance cecentricity - Fx, Fy fluid film force components - fi = Fi/W dimensionless fluid film force components - L bearing length - M one half rotor mass - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbnL2yY9% 2CVzgDGmvyUnhitvMCPzgarmWu51MyVXgaruWqVvNCPvMCG4uz3bqe% fqvATv2CG4uz3bIuV1wyUbqee0evGueE0jxyaibaieYlf9irVeeu0d% Xdh9vqqj-hEeeu0xXdbba9frFf0-OqFfea0dXdd9vqaq-JfrVkFHe9% pgea0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaaca% qabeaadaabauaaaOqaaiaad2gacqGH9aqpcaWGnbGaam4qaiabeM8a% 3naaCaaaleqabaGaaGOmaaaakiaac+cacqaHdpWCcaWGxbaaaa!471F!\[m = MC\omega ^2 /\sigma W\] dimensionless one half rotor mass - R bearing radius - T = 2 synchronous orbit period - t time - W load per bearing - X, Y, Z coordinates - x = X/C; y = y/C; z = Z/L dimensionless coordinates - oil dynamic viscosity - = E/C dimensionless mass unbalance eccentricity - = (RL/W)/(R/C) ² (L/D) ² modified Sommerfeld number - = t dimensionless time = periodic orbit frequency - = 2/ frequency ratio - journal angular velocity - (·) dimensionless time derivative     

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 an inversion theorem for conical regions in elasticity theory

Summary Two-dimensional stress singularities in wedges have already drawn attention since a long time. An inverse square-root stress singularity (in a 360° wedge) plays an important role in fracture mechanics.Recently some similar three-dimensional singularities in conical regions have been investigated, from which one may be also important in fracture mechanics.Spherical coordinates are r, , . The conical region occupied by the elastic homogeneous body (and possible anisotropic) has its vertex at r=0. The mantle of the cone is described by an arbitrary function f(, )=0. The displacement components be u. For special values of (eigenvalues) there exist states of displacements (eigenstates) % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],which may satisfy rather arbitrary homogeneous boundary conditions along the generators.The paper brings a theorem which expresses that if is an eigenvalue, then also-1- is an eigenvalue. Though the theorem is related to a known theorem in Potential Theory (Kelvin's theorem), the proof has to be given along quite another line.

---

Zusammenfassung Zwei-dimensionale Spannungssingularitäten in keilförmigen Gebieten sind schon längere Zeit untersucht worden und neuerdings auch ähnliche drei-dimensionale Singularitäten in konischen Gebieten.Kugelkoordinaten sind r, , . Das konische Gebiet hat seine Spitze in r=0. Der Mantel des Kegels lässt sich beschreiben mittels einer willkürlichen Funktion f(, )=0. Die Verschiebungskomponenten seien u. Für spezielle Werte von (Eigenwerte) bestehen Verschiebunszustände % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],welche homogene Randwerte der Beschreibenden des Kegels entlang genügen.Das Bericht bringt ein Theorem, welches aussagt, das und =–1– beide Eigenwerte sind.

     

Asymptotics of the homogenized moduli for the elastic chess-board composite

We find the asymptotic behavior of the homogenized coefficients of elasticity for the chess-board structure. In the chess board white and black cells are isotropic and have Lamé constants (, ,) and (, ) respectively. We assume that the black cells are soft, so 0. It turns out that the Poisson ratio for this composite tends to zero with .     

Uniqueness theorems for displacement fields with locally finite energy in linear elastostatics

Uniquencess theorems are proved for the fundamental boundary value problems of linear elastostatics in bodies of arbitrary shape. The displacement fields are required to have finite strain energy in bounded portions of the bodies and satisfy the principle of virtual work. For bounded bodies, the total strain energy is finite and uniquencess is proved without additional hypotheses. In particular, no restrictions other than the energy condition are placed on the field singularities that may occur at sharp edges and corners. For unbounded bodies, uniqueness can be proved as in the bounded case if the total strain energy is finite. Sufficient conditions for this are shown to be the finiteness of the strain energy in bounded portions of the body together with the growth restriction % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakeaadaWdraqaaiaabwhadaWg% aaWcbaGaaeyAaaqabaGccaGGOaGaaeiEaiaacMcacaqG1bWaaSbaaS% qaaiaabMgaaeqaaOGaaiikaiaabIhacaGGPaGaaeizaiaabIhacaqG% 9aGaaGimaiaacIcacaqGYbGaaiykaiaacYcacaqGYbGaeyOKH4Qaey% OhIukaleaacqGHPoWvdaWgaaadbaGaaeOCaiaacYcacqaH0oazaeqa% aaWcbeqdcqGHRiI8aaaa!5E73!\[\int_{\Omega _{{\text{r}},\delta } } {{\text{u}}_{\text{i}} ({\text{x}}){\text{u}}_{\text{i}} ({\text{x}}){\text{dx = }}0({\text{r}}),{\text{r}} \to \infty } \] on the displacement fieldu _i , where _r, is the portion of the body that lies between concentric spheres with radiir andr+ and >0.This research was supported by the Air Force Office of Scientific Research. Reproduction in whole or part is permitted for any purpose of the United States Government.Prepared under Contract No. F 49620-77-C-0053 for Air Force Office of Scientific Research.     

Conditions for noncircular whirling of nonlinear isotropic rotors

Nonlinear rotors are often considered as potential sources of chaotic vibrations. The aim of the present paper is that of studying in detail the behaviour of a nonlinear isotropic Jeffcott rotor, representing the simplest nonlinear rotor. The restoring and damping forces have been expanded in Taylor series obtaining a Duffing-type equation. The isotropic nature of the system allows circular whirling to be a solution at all rotational speeds. However there are ranges of rotational speed in which this solution is unstable and other, more complicated, solutions exist.The conditions for stability of circular whirling are first studied from closed form solutions of the mathematical model and then the conditions for the existence of solutions of other type are studied by numerical experimentation. Although attractors of the limit cycle type are often found, chaotic attractors were identified only in few very particular cases. An attractor supposedly of the last type reported in the literature was found, after a detailed analysis, to be related to a nonchaotic polyharmonic solution.As the typical unbalance response of isotropic nonlinear rotors has been shown to be a synchronous circular whirling motion, the convergence characteristics of Newton-Raphson algorithm applied to the solution of the set of nonlinear algebraic equations obtained from the differential equations of motion are studied in some detail. c damping coefficient i imaginaty unit (i=% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbqfgBHr% xAU9gimLMBVrxEWvgarmWu51MyVXgaruWqVvNCPvMCG4uz3bqefqvA% Tv2CG4uz3bIuV1wyUbqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9% vqqj-hEeeu0xXdbba9frFf0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea% 0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabe% aadaabauaaaOqaamaakaaabaGaeyOeI0IaaGymaaWcbeaaaaa!3E66!\[\sqrt { - 1}\]) k stiffness m mass t time x _istate variables i=1, 4 z complex co-ordinate (z=x+iy)[J] Jacobian matrix Oxyz inertial co-ordinate frame Oz rotating co-ordinate frame perturbation term eccentricity complex co-ordinate (=+i) system eigenvalues nonlinearity parameter nondimensional time phase spin speed u nonrotating t rotating0 amplitude t nondimensional termsNomenclature     

Electrical parameters of a channel with finite electrodes with allowance for the electrode potential drop

Spatial problems involving the electric field in an MHD channel were formulated in [1] with allowance for the electrode potential drop. It was assumed that the electrode layer had a small thickness, so that relationships on the boundary of the layer could be applied to the surface of the electrode. It was assumed that the electrode potential drop ° could be represented as a function of the current density jn at the electrode in the form of a known function ° =f (j_n) determined experimentally or deduced from the appropriate electrode-layer theory. An approximate method was then put forward for solving such problems by reducing them to the determination of the electric field from a known distribution of the magnetic field and the gas-dynamic parameters. It was shown that when =°/ E is small (E is the characteristic induced or applied potential difference), the solution can be sought in the form of series in powers of . In the zero-order approximation, the electric field is determined without taking into account the electrode processes. The first approximation gives a correction of the order of . The quantity °, which is present in the boundary conditions on the electrode in the first-order approximation, is determined from the current density calculated in the zero-order approximation.One of the problems discussed in [1] was concerned with the electric current in a channel with one pair of symmetric electrodes. Its solution was found in the first approximation in the form of the integral Keldysh-Sedov formula. In this paper we report an analysis of the solution for ° taken in the form of a step function.     

Alloy separation of a binary mixture in a stressed elastic sphere

We consider the static state of a spherical isotropic binary elastic solid mixture whose boundary is given a uniform radial displacement. The elastic volumetric strain energy is given by the classical quadratic form from linear elasticity theory,% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakeaacaqGxbGaeyypa0ZaaSaa% aeaacaGGXaaabaGaaiOmaaaacaGG7bGaaCOUdiaacIcacaqGJbGaam% ykaiaadIcacaWG0bGaamOCaiaadwgacaWGPaWaaWbaaSqabeaacaGG% YaaaaOGaey4kaSIaaiOmaiabeY7aTjaacIcacaqGJbGaaiykaiaacY% hacaWGLbGaeyOeI0IaaiiiamaalaaabaGaamymaaqaaiaadodaaaGa% aiikaiaadshacaWGYbGaamyzaiaacMcacaWGXaGaaiiFamaaCaaale% qabaGaaiOmaaaakiaac2hacaGGUaaaaa!63E0!\[{\text{W}} = \frac{1}{2}\{ {\mathbf{\kappa }}({\text{c}})(tre)^2 + 2\mu ({\text{c}})|e - \frac{1}{3}(tre)1|^2 \} .\]Here, e is the infinitesimal strain tensor, c[0, 1] is the volumetric concentration of the mixture, and (·) and (·) are the (positive) bulk and shear material moduli, respectively, which are given functions of the concentration. As a function of c and e, the strain energy function is generally nonconvex. Thus, we consider the nonconvex problem of minimizing the potential energy of the body, among all spatial concentration and displacement fields, subject to a given boundary displacement and a fixed amount of component materials. Assuming spherical symmetry, we find that the two component materials must be separated in the optimal state of minimum potential energy. The harder material forms the central core of the sphere, and the softer material is segregated into a surrounding shell. This behavior is remindful of a general notion in metallurgy that in the casting of materials the harder material tends to migrate toward the center.Partial support of the NSF under grant MSS-9024637 and Alliant Techsystems Inc. is gratefully acknowledged.Professor R. Bartelletti of the Università di Pisa is gratefully acknowledged.     

Methods of Small and Large δ in the Nonlinear Dynamics – A Comparative Analysis

New asymptotic approaches for dynamical systems containing a power nonlinear term x ⁿ are proposed and analyzed. Two natural limiting cases are studied: n 1 + , 1 and n . In the firstcase, the 'small method' (SM)is used and its applicability for dynamical problems with the nonlinearterm sin as well as the usefulness of the SMfor the problem with small denominators are outlined. For n , a new asymptotic approach is proposed(conditionally we call it the 'large method' –LM). Error estimations lead to the followingconclusion: the LM may be used, even for smalln, whereas the SM has a narrow application area. Both of the discussed approaches overlap all values ofthe parameter n.     

The virtual origin of riblets in an open channel flow

In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f _r skin friction coefficient - f _ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H ⁺ u*H/ - H false water depth from top of riblets to water surface - H ⁺ u*H/ - I _e streamwise energy slope - I _b bed slope - k riblet height - k ⁺ u*k/ - l mixing length - l _s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y ₀ distance from top of riblet to virtual origin - y _v distance from top of riblet to edge of viscous layer - y ⁺ u*y/ - y ⁺ u*y/ - y ₀ ⁺ u*y ₀/ - u ⁺ u*y/ - shifting coefficient for standardization - thickness of viscous layer=y ₀+y - ⁺ u*/ - ⁺ u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress     

Effects of viscous dissipation on heat transfer in an oscillating flow past a flat plate

Theoretical investigation has been carried out of laminar thermal boundary layer response to harmonic oscillations in velocity associated with a progressive wave imposed on a steady free stream velocity and convected in the free stream direction. Series solutions are derived both to velocity and temperature field and the resulting equations are solved numerically. The functions affecting the temperature field are shown graphically for different values of Prandtl number. It is observed that there is more reduction in the rate of heat transfer for P _r<1 and a rise in the rate of heat transfer for P _r>1 due to the presence of oscillatory free-stream.Nomenclature u, v velocity components in the x and y direction - x, y Cartesian coordinate axes - t time - U, U ₀ instantaneous value of and mean free stream velocity - density of fluid - kinematic viscosity - T, T _w, T temperature of the fluid, wall and free stream fluid - c _p specific heat at constant pressure - thermal diffusivity - amplitude of free stream velocity - frequency - _p non-dimensional temperature (T–T /T _w–T ) - P _r Prandtl number (c _p/K) - E _c Eckert number (U ₀ ² /c _p(T _w–T )) - a parameter ( ) - ₀ boundary layer thickness of the oscillation of a harmonic oscillation of frequency ( ) - ordinary boundary layer thickness ( ) - time-averaged, time-independent external velocity - A, B, C, D, E, K, L, M, N, P functions used in expansion for u and - Nu Nusselt number (hx/k) - T _w–% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8V4rqqrFfpeea0Jc9yq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepGe9fr-xfr-x% frpeWZqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaacIcadaGcaa% qaaiaadAhacaWG4bGaai4laiqadwfagaqeaaWcbeaakiaacMcaaaa!3CA6!\[(\sqrt {vx/\bar U} )\] - k thermal conductivity     

An improved method for the measurement of strong flow birefringence

An improved procedure is proposed to solve the problem of measuring strong birefringence on flowing samples which have a retardation larger than /2 The procedure can be applied to devices which are based on the modulation of the polarization vector with a high-speed rotating half-wave plate. It uses two optical configurations which have different responses to the retardation. By combining the information obtained from the two configurations, one can easily and accurately determine the actual retardation of the sample, irrespective of the quadrant in which the measured retardation is located, provided the sample is isotropic. This technique can also overcome the problem of large errors occurring in the vicinity of = (2m+1)/2 due to the limited sensitivity around these values of . The use of the technique is illustrated with measurements on an isotropic solution of poly(benzyl glutamate) in m-cresol, which is strongly birefringent during shear flow.     

Heat transfer between gas-liquid spray stream flowing perpendicularly to the row of the cylinders

Experimental investigation and analysis of heat transfer process between a gas-liquid spray flow and the row of smooth cylinders placed in the surface perpendicular to the flow has been performed. Among others, there was taken into account in the analysis the phenomenon of droplets bouncing and omitting the cylinder as well as the phenomenon of the evaporation process from the liquid film surface.In the experiments test cylinders were used, which were placed between two other cylinders standing in the row.From the experiments and the analysis the conclusion can be drawn that the heat transfer coefficients values for a row of the cylinders are higher than for a single cylinder placed in the gasliquid spray flow.

---

Wärmeübergang an eine senkrecht anf eine Zylinderreihe auftreffende Gas-Flüssigkeits-Sprüh-Strömung

Zusammenfassung Es wurden Messungen und theoretische Analysen des Wärmeübergangs zwischen einer Gas-FlüssigkeitsSprüh-Strömung und den glatten Oberflächen einer Zylinderreihe durchgeführt, die senkrecht zum Sprühstrahl angeordnet waren. Dabei wurde in der Analyse unter anderem das Phänomen betrachtet, daß die Tropfen die Zylinderwand treffen und verfehlen können und daß sich ein Verdampfungsprozeß aus dem flüssigen Film an der Zylinderoberfläche einstellt.Gemessen wurde an einem zwischen zwei Randzylindern befindlichen Zylinder.Die Experimente und die Analyse gestatten die Schlußfolgerung, daß der Wärmeübergangskoeffizient für eine Zylinderreihe höher ist als für einen einzelnen Zylinder in der Sprühströmung.

---

Nomenclature a distance between axes of cylinders, m - c _l specific heat capacity of liquid, J/kg K - c _g specific heat capacity of gas, J/kg K - D cylinder diameter, m - g _l mass velocity of liquid, kg/m²s - ¯k average volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - k() local volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - L specific latent heat of vaporisation, J/kg - m mass fraction of water in gas-liquid spray flow, dimensionless - M constant in Eq. (9) - p pressure, Pa - p _g statical pressure of gas, Pa - p _w pressure of gas on the cylinder surface, Pa - p external pressure on the liquid film surface, Pa - r cylindrical coordinate, m - R radius of cylinder, m - T temperature, K, °C - T _l, t_l liquid temperature in the gas-liquid spray, K, °C - T _w,t_w temperature of cylinder surface, K, °C - T temperature of gas-liquid film interface, K - U liquid film velocity, m/s - w gas velocity on cylinder surface, m/s - w _g gas velocity in free stream, m/s - W _l liquid vapour mass ratio in free stream, dimensionless - W liquid vapour mass ratio at the edge of a liquid film, dimensionless - x coordinate, m - y coordinate, m - z complex variable, dimensionless - average heat transfer coefficient, W/m²K - local heat transfer coefficient, W/m² K - average heat transfer coefficient between cylinder surface and gas, W/m² K - _g, local heat transfer coefficient between cylinder surface and gas, W/m² K - mass transfer coefficient, kg/m²s - liquid film thickness, m - _lg dynamic diffusion coefficient of liquid vapour in gas, kg/m s - pressure distribution function on a cylinder surface - function defined by Eq. (3) - _l liquid dynamic viscosity, kg/m s - _g gas dynamic viscosity, kg/m s - cylindrical coordinate, rad, deg - _l thermal conductivity of liquid, W/m K - _g thermal conductivity of gas, W/m K - mass transfer driving force, dimensionless - _l density of liquid, kg/m³ - _g density of gas, kg/m³ - _w shear stress on the cylinder surface, N/m² - _w shear stress exerted by gas at the liquid film surface, N/m² - air relative humidity, dimensionless - T -T _w - _w =T _wT_l Dimensionless parameters I= enhancement factor of heat transfer - m ^*=M _l/M_g molar mass of liquid to the molar mass of gas ratio - Nu _g= D/ _g gas Nusselt number - Pr _g=c _{g g}/_g gas Prandtl number - Pr _l=c_lll liquid Prandtl number - ¯r=(r–R)/ dimensionless coordinate - Re _g=w_gD _g/_g gas Reynolds number - Re _g,max=w_g,max D _g/_g gas Reynolds number calculated for the maximal gas velocity between the cylinders - Sc=m ^* _g/_l–g Schmidt number =/R dimensionless film thickness     

Asymptotic Nusselt numbers for Newtonian flow through a parallel-plate channel with recycle

General expressions for evaluating the asymptotic Nusselt number for a Newtonian flow through a parallel-plate channel with recycle at the ends have been derived. Numerical results with the ratio of thicknesses as a parameter for various recycle ratios are obtained. A regression analysis shows that the results can be expressed by log Nur^0.83=0.3589 (log)² -0.2925 (log) + 0.3348 forR 3, 0.1 0.9; logNu=0.5982(log)² +0.3755 × 10^–2 (log) +0.8342 forR 10^–2, 0.1 0.9.

---

Asymptotische Nusselt-Zahlen für die Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung

Zusammenfassung In dieser Untersuchung wurden allgemeine Ausdrücke hergeleitet um die asymptotische Nusselt-Zahl für eine Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung an den Enden berechnen zu können. Es wurden numerische Ergebnisse mit den Dicken-Verhältnissen, als Parameter für verschiedene Rückführungs-verhältnisse, erhalten. Eine Regressionsanalyse zeigt, daß die Ergebnisse wie folgt ausgedrückt werden können: log Nur^0,83=0,3589 (log)² -0,2925 (log) + 0,3348 fürR 3, 0,1 0,9; logNu=0,5982(log)² +0,3755 × 10^–2 (log) + 0,8342 fürR 10^–2, 0,1 0,9.

---

Nomenclature A₁ shooting value,d(0)/d - A₂ shooting value,d(1)/d - B channel width - Gz Graetz number, U_bW²/L - h _m logarithmic average convective heat transfer coefficient - h _x average local convective heat transfer coefficient - k thermal conductivity - L channel length - Nu average local Nusselt number, 2 h_xW/k - Nu _m logarithmic average Nusselt number, 2h_mW/k - R recycle ratio, reverse volume flow rate divided by input volume flow rate - T temperature of fluid - T _m bulk temperature, Eq. (8) - T ₀ temperature of feed stream - T _s wall temperature - U velocity distribution - U _b reference velocity,V/BW - V input volume flow rate - v dimensionless velocity distribution, U/U_b - W channel thickness - x longitudinal coordinate - y transversal coordinate - Z₁-z₆ functions defined in Eq. (A1) - thermal diffusivity - least squares error, Eq. (A7) - weight, Eqs. (A8), (A9) - dimensionless coordinate,y/W - dimensionless coordinate,x/GzL - function, Eq. (7)     

A similarity law for acute circular cones at large angles of attack in a supersonic flow

For thin bodies placed in a hypersonic flow at a small angle of attack the similarity law is known. From this law it follows that for various numbers M, angles of attack , and relative thicknesses the similarity conditions will be observed if in the flows under consideration the parameters M and / are the same. This similarity law is obtained with the assumption M 1, 1. But even for M=3 and 1/3 the results of solving the complete system of gasdynamic equations for affino-similar bodies is in a good agreement with the similarity law [1], In [2] it is shown that this similarity law is generalized for the case of a flow around a thin pointed body at large angles of attack. According to the similarity law, at large angles of attack the flows near bodies with an identical distribution of cross-sectional shapes will be similar if the parameters K₁= cotan and K₂=m sin for all cases have one and the same value. As the angle of attack decreases, the requirements of constancy of K₁ and K₂ become analogous to the conditions M=const, /=const.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 78–83, May–June, 1976.The authors thank V. V. Lunev for the useful discussions and valuable observations.     

ON THE UNIFICATION OF THE HAMILTON PRINCIPLES IN NONHOLONOMIC SYSTEM AND IN HOLONOMIC SYSTEM

ＯＮＴＨＥＵＮＩＦＩＣＡＴＩＯＮＯＦＴＨＥＨＡＭＩＬＴＯＮＰＲＩＮＣＩＰＬＥＳＩＮＮＯＮＨＯＬＯＮＯＭＩＣＳＹＳＴＥＭＡＮＤＩＮＨＯＬＯＮＯＭＩＣＳＹＳＴＥＭ（梁立孚）（韦扬）ＯＮＴＨＥＵＮＩＦＩＣＡＴＩＯＮＯＦＴＨＥＨＡＭＩＬＴＯＮＰＲＩＮＣＩＰＬ...     

Response of an elastic Bingham fluid to oscillatory shear

The response of an elastic Bingham fluid to oscillatory strain has been modeled and compared with experiments on an oil-in-water emulsion. The newly developed model includes elastic solid deformation below the yield stress (or strain), and Newtonian flow above the yield stress. In sinusoidal oscillatory deformations at low strain amplitudes the stress response is sinusoidal and in phase with the strain. At large strain amplitudes, above the yield stress, the stress response is non-linear and is out of phase with strain because of the storage and release of elastic recoverable strain. In oscillatory deformation between parallel disks the non-uniform strain in the radial direction causes the location of the yield surface to move in-and-out during each oscillation. The radial location of the yield surface is calculated and the resulting torque on the stationary disk is determined. Torque waveforms are calculated for various strains and frequencies and compared to experiments on a model oil-in-water emulsion. Model parameters are evaluated independently: the elastic modulus of the emulsion is determined from data at low strains, the yield strain is determined from the phase shift between torque and strain, and the Bingham viscosity is determined from the frequency dependence of the torque at high strains. Using these parameters the torque waveforms are predicted quantitatively for all strains and frequencies. In accord with the model predictions the phase shift is found to depend on strain but to be independent of frequency.Notation A plate strain amplitude (parallel plates) - A _R plate strain amplitude at disk edge (parallel disks) - G elastic modulus - m torque (parallel disks) - M normalized torque (parallel disks) = 2m/R ³₀ - N ratio of viscous to elastic stresses (parallel plates) =µ A/ ₀ ratio of viscous to elastic stresses (parallel disks) =µ A _R/₀ - r normalized radial position (parallel disks) =r/R - r radial position (parallel disks) - R disk radius (parallel disks) - t normalized time = t — /2 - t time - _E elastic strain - _P plate strain (displacement of top plate or disk divided by distance between plates or disks) - _PR plate strain at disk edge (parallel disks) - ₀ yield strain - _E normalized elastic strain = _E/₀ - _P normalized plate strain = _P/₀ - _PR normalized plate strain at disk edge (parallel disks) = _PR/₀ - ₀ normalized plate strain amplitude (parallel plates) =A/ ₀ — normalized plate strain amplitude at disk edge (parallel disks) =A _R/₀ - phase shift between _P andT (parallel plates) — phase shift between _PR andM (parallel disks) - µ Bingham viscosity - stress - ₀ yield stress - T normalized stress =/ ₀ - frequency     

Schwere symmetrische Kreisel kleiner Nutationen

Zusammenfassung Zur Integration der Eulerschen Bewegungsgleichungen schwerer symmetrischer Kreisel werden der Winkel (t) (Abb. 1) durch (t)=₀+(t) ersetzt und in sämtlichen Reihenentwicklungen von abhängiger Funktionen die Potenzen höheren als zweiten Grades vernachlässigt. Dadurch ist es möglich, die Eulerschen Winkel (t), (t) und (t) durch elementare Formeln zu beschreiben und somit sind die wesentlichsten Erscheinungen im Bewegungsablauf der schweren symmetrischen Kreisel einfach zu übersehen.     

Transcendentally small transversality in the rapidly forced pendulum

The rapidly forced pendulum equation with forcing sin((t/), where =<₀^p,p = 5, for ₀, sufficiently small, is considered. We prove that stable and unstable manifolds split and that the splitting distanced(t) in the ( ,t) plane satisfiesd(t) = sin(t/) sech(/2) +O( ₀ exp(–/2)) (2.3a) and the angle of transversal intersection,, in thet = 0 section satisfies 2 tan/2 = 2S _s = (/2) sech(/2) +O(( ₀ /) exp(–/2)) (2.3b) It follows that the Melnikov term correctly predicts the exponentially small splitting and angle of transversality. Our method improves a previous result of Holmes, Marsden, and Scheuerle. Our proof is elementary and self-contained, includes a stable manifold theorem, and emphasizes the phase space geometry.