On forced vibration of anisotropic cylinders

Summary The dynamic response of a circular cylinder with thick walls of transverse curvilinear isotropy subjected to a uniformly distributed pressure varying periodically with time is analyzed by means of the Laplace transformation, and the exact solution is obtained in closed form. The previously obtained solutions for forced vibrations with isotropy, and free vibrations with transverse curvilinear isotropy are included as special cases of the general results reported here.Nomenclature t time - r, , z cylindrical coordinates - _ii components of normal strain - _ii components of normal stress - u radial displacement - c _ij elastic constant - mass density - c ² c ₁₁/ - ² c ₂₂/c ₁₁ - a, b inner, outer radius of the cylinder - , A, B constants - forced angular frequency - function defined by (9) - p, real, complex variables - constant defined by (14) - real number - , Lamé elastic constants - J (x) Bessel function of first kind - Y (x) Bessel function of second kind - I (x) modified Bessel function of first kind - K (x) modified Bessel function of second kind     

Critical layer formation in a stratified medium with mean shear flows

The problem of the excitation of internal waves with a given wave number k and frequency in a stratified medium with shear flows is considered. The internal wave field of the form v(z)exp(–it+ikx) established as t in a medium without dissipation has a singular point at the level z=z₀ (critical level), at which the flow velocity U(z) coincides with the phase velocity /k. Dissipative effects (viscosity and heat conduction) smooth out this singularity. An exact solution of the model equation describing as t and zz₀ the field excited by oscillating sources activated at t=0 is constructed with allowance for dissipation. This makes it possible to describe the limiting steady-state field, determine the critical layer as the neighborhood of the critical level in which dissipation effects are important, and to estimate its width and the rate of convergence to the limiting steady-state regime. The asymptotic behavior of the fields is examined for Ri1, where Ri is the Richardson number. It is shown that when the well-known Miles stability condition Ri>1/4 is satisfied there are no natural oscillations with a critical level.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 82–93, May–June, 1990.     

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.

     

Effect of measuring volume length on two-component laser velocimeter measurements in a turbulent boundary layer

The effects of finite measuring volume length on laser velocimetry measurements of turbulent boundary layers were studied. Four different effective measuring volume lengths, ranging in spanwise extent from 7 to 44 viscous units, were used in a low Reynolds number (Re=1440) turbulent boundary layer with high data density. Reynolds shear stress profiles in the near-wall region show that u v strongly depends on the measuring volume length; at a given y-position, u v decreases with increasing measuring volume length. This dependence was attributed to simultaneous validations on the U and V channels of Doppler bursts coming from different particles within the measuring volume. Moments of the streamwise velocity showed a slight dependence on measuring volume length, indicating that spatial averaging effects well known for hot-films and hot-wires can occur in laser velocimetry measurements when the data density is high.List of symbols time-averaged quantity - u wall friction velocity, ( _w /)^1/2 - v kinematic viscosity - d _p pinhole diameter - l _eff spanwise extent of LDV measuring volume viewed by photomultiplier - l ⁺ non-dimensional length of measuring volume, l _eff u /v - y ⁺ non-dimensional coordinate in spanwise direction, y u /v - z ⁺ non-dimensional coordinate in spanwise direction, z u /v - U ⁺ non-dimensional mean velocity, /u - u instantaneous streamwise velocity fluctuation, U &#x2329;U - v instantaneous normal velocity fluctuation, V–V - u RMS streamwise velocity fluctuation, u ^21/2 - v RMS normal velocity fluctuation, v ^21/2 - Re Reynolds number based on momentum thickness, U 0/v - R _uv cross-correlation coefficient, u v/u v - R₁₂(0, 0, z) two point correlation between u and v with z-separation, <u(0, 0, 0) v (0, 0, z)>/<u(0, 0, 0) v (0, 0, 0)> - N rate at which bursts are validated by counter processor - T Taylor time microscale, u (dv/dt²)^–1/2     

Rheology of molten polymers

Summary TheCross equation describes the flow of pseudoplastic liquids in terms of an upper and a lower Newtonian viscosity corresponding to infinite and zero shear, and ₀, and of a third material constant related to the mechanism of rupture of linkages between particles in the intermediate, non-Newtonian flow regime, Calculation of of bulk polymers is important, since it cannot be determined experimentally. The equation was applied to the melt flow data of two low density polyethylenes at three temperatures.Using data in the non-Newtonian region covering 3 decades of shear rate to extrapolate to the zero-shear viscosity resulted in errors amounting to about onethird of the measured ₀ values. The extrapolated upper Newtonian viscosity was found to be independent of temperature within the precision of the data, indicating that it has a small activation energy.The ₀ values were from 100 to 1,400 times larger than the values at the corresponding temperatures.The values of were large compared to the values found for colloidal dispersions and polymer solutions, but decreased with increasing temperature. This shows that shear is the main factor in reducing chain entanglements, but that the contribution of Brownian motion becomes greater at higher temperatures.

---

Zusammenfassung Die Gleichung vonCross beschreibt das Fließverhalten von pseudoplastischen Flüssigkeiten durch drei Konstante: Die obereNewtonsche Viskosität (bei sehr hohen Schergeschwindigkeiten), die untereNewtonsche Viskosität ₀ (bei Scherspannung Null), und eine Materialkonstante, die vom Brechen der Bindungen zwischen Partikeln im nicht-Newtonschen Fließbereich abhängt. Die Berechnung von ist wichtig für unverdünnte Polymere, wo man sie nicht messen kann.Die Gleichung wurde auf das Fließverhalten der Schmelzen von zwei handelsüblichen Hochdruckpolyäthylenen bei drei Temperaturen angewandt. Die Werte von ₀, durch Extrapolation von gemessenen scheinbaren Viskositäten im Schergeschwindigkeitsbereich von 10 bis 4000 sec^–1 errechnet, wichen bis 30% von den gemessenen ₀-Werten ab. Die Aktivierungsenergie der war so klein, daß die-Werte bei den drei Temperaturen innerhalb der Genauigkeit der Extrapolation anscheinend gleich waren. Die ₀-Werte waren 100 bis 1400 mal größer als die-Werte.Im Verhältnis zu kolloidalen Dispersionen und verdünnten Polymerlösungen war das der Schmelzen groß, nahm aber mit steigender Temperatur ab. Deshalb wird die Verhakung der Molekülketten hauptsächlich durch Scherbeanspruchung vermindert, aber der Beitrag derBrownschen Bewegung nimmt mit steigender Temperatur zu.

     

A turbulence model for the heat transfer near stagnation point of a circular cylinder

A one-equation low-Reynolds number turbulence model has been applied successfully to the flow and heat transfer over a circular cylinder in turbulent cross flow. The turbulence length-scale was found to be equal 3.7y up to a distance 0.05 and then constant equal to 0.185 up to the edge of the boundary layer (wherey is the distance from the surface and is the boundary layer thickness).The model predictions for heat transfer coefficient, skin friction factor, velocity and kinetic energy profiles were in good agreement with the data. The model was applied for Re 250,000 and Tu0.07.Nomenclature µ,C _D Constants in the turbulence kinetic energy equation - C ₁,C ₂ Constants in the turbulence length-scale equation - Skin friction coefficient atx - D Cylinder diameter - F Dimensionless flow streamwise velocityu/u _e - k Turbulence kinetic energy =1/2 the sum of the squared three fluctuating velocities - K Dimensionless turbulence kinetic energyk/u _e ^/2 - I Dimensionless temperature (T–T _w )/(T –T _w ) - l Turbulence length-scale - l _e Turbulence length-scale at outer region - Nu _D Nusselt number - p Pressure - Pr Prandtl number - Pr _t Turbulent Prandtl number - Pr _k Constant in the turbulence kinetic energy equation - R Cylinder radius - Re _D Reynolds number u D/µ - Re _x Reynolds number u x/µ - R _K Reynolds number of turbulence - T Mean temperature - T Mean temperature at ambient - T _s Mean temperature at surface - Tu Cross flow turbulence intensity, - u Mean flow streamwise velocity - u Fluctuating streamwise velocity - u _e Mean flow velocity at far field distance - u Mean flow velocity at ambient - u* Friction velocity - v Mean velocity normal to surface - V Dimensionless mean velocity normal to surface - x,x ₁ Distance along the surface - y Distance normal to surface - Dimensionless pressure gradient parameter - Boundary layer thickness atu=0.9995u _e - Transformed coordinate iny direction - Fluid molecular viscosity - _t Turbulent viscosity - _eff + _t - µ Fluid molecular viscosity at ambient - Kinematic viscosity/ - Density - Density at ambient - _w Wall shear stress - _w,0 Wall shear stress at zero free stream turbulence     

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

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

Influence of Porous Media and Airflow Rate on the Fate of NAPLs Under Air Sparging

Nonequilibrium air–water mass transfer experiments using a laboratoryscale singleair channel setup were conducted to investigate the influence of porous media and air velocity on the fate of nonaqueous phase liquids (NAPLs) under air sparging conditions. Benzene was used as a NAPL while silica sand 30/50 (dp₅₀=0.305mm, uniformity coefficient, UC=1.41) and silica sand 70/100 (dp₅₀=0.168mm, UC=1.64) were used as porous media. Air velocities ranged from 0 to 1.4cm/s. Mass transfer coefficients for the dissolution of NAPLs were estimated by numerical methods using a twodimensional dissolution–diffusion–volatilization model. The study showed that the presence of advective airflow in air channels controlled the spreading of the dissolved phase but the overall removal efficiency was independent of airflow rate. Removal efficiencies and dissolution rates of the NAPL were found to be strongly affected by the mean particle size of the porous media during air sparging. More than 50% reduction in the removal rate of benzene was found when silica sand 70/100 was used instead of silica sand 30/50. Mass transfer coefficients for the dissolution of benzene NAPL were estimated to be 0.0041cm/min for silica sand 70/100 and 0.227cm/min for silica sand 30/50. Increasing the air velocity from 0.6 to 1.4cm/s for silica sand 30/50 did not result in a higher removal rate. Quantitative estimation of the dissolution rates of benzene NAPL indicated that the dissolution rates (between 0.227 and 0.265cm/min) were similar in magnitude for the same porous media but different air flow rates. Based on the visualization study, air sparging may be used to control the spreading of the dissolved phase even when the glob of NAPL is several centimeters away from the air–water interface of the air channels.     

Hartmann Effect. Region of Existence and Oscillation Frequencies

Data of experiments on determining the region of existence for autooscillations (Hartmann effect) with a frequency of approximately 450 Hz with variation in the distance from the nozzle to the resonator. Results of these experiments differ from the wellknown results of Hartmann and his followers obtained for shallow resonators. It is shown that the region of autooscillations exists for large distances between the nozzle and the resonator. The results obtained are explained using the modern knowledge of the gasdynamic structure of a supersonic underexpanded jet. It is shown that in determining the frequency of lowfrequency oscillations, it suffices to allow for the resonator length and its added mass.     

Indentation of the semi-infinite elastic solid by a concave rigid punch

A solution is obtained for the relationship between load, displacement and inner contact radius for an axisymmetric, spherically concave, rigid punch, indenting an elastic half-space. Analytic approximations are developed for the limiting cases in which the ratio of the inner and outer radii of the annular contact region is respectively small and close to unity. These approximations overlap well at intermediate values. The same method is applied to the conically concave punch and to a punch with a central hole. , , . , . . .     

A non-dualistic unified field theory of gravitation,electromagnetism and spin

The wisdom of classicalunified field theories in the conceptual framework of Weyl, Eddington, Einstein and Schrödinger has often been doubted and in particular there does not appear to be any empirical reason why the Einstein-Maxwell (E-M) theory needs to be geometrized. The crux of the matter is, however not whether the E-M theory is aesthetically satisfactory but whether it answers all the modern questions within the classical context. In particular, the E-M theory does not provide a classical platform from which the Dirac equation can be derived in the way Schrödinger's equation is derived from classical mechanics via the energy equation and the Correspondence Principle. The present paper presents a non-dualistic unified field theory (UFT) in the said conceptual framework as propounded by M. A. Tonnelat. By allowing the metric formds ²=g dx ^v x ^v and the non-degenerate two-formF=(1/2> _l) dx ^vdx ^vto enter symmetrically into the theory we obtain a UFT which contains Einstein's General Relativity and the Born-Infeld electrodynamics as special cases. Above all, it is shown that the Dirac equation describing the electron in an external gravito-electromagnetic field can be derived from the non-dualistic Einstein equation by a simple factorization if the Correspondence Principle is assumed.     

On the Cρ-Smoothness of an Invariant Torus of a Countable System of Difference Equations Defined on an m-Dimensional Torus

Using the method of truncation, we establish sufficient conditions of differentiability of order 2 with respect to the angular variable for an invariant torus of a countable system of difference equations with deviation of the discrete argument.     

A new interpretation of viscosity and yield stress in dense slurries: Coal and other irregular particles

Illinois coal was ground and wet-sieved to prepare three powder stocks whose particle-size distributions were characterized. Three suspending fluids were used (glycerin, bromonaphthalene, Aroclor), with viscosities _s that differed by a factor of 100 and with very different chemistries, but whose densities matched that of the coal. Suspensions were prepared under vacuum, with coal volume fractions that ranged up to 0.46. Viscosities were measured in a cone-and-plate over a shear rate range 10^–3–10² s^–1. Reduced viscosity _r = /_s is correlated in the high-shear limit ( ) with/ _M, where _M is the maximum packing fraction for the high-shear microstructure, to reveal the roles of size distribution and suspending fluid character. A new model that invokes the stress-dependence of _M is found to correlate _r well under non-Newtonian conditions with simultaneous prediction of yield stress at sufficiently high; a critical result is that stress and not governs the microstructure and rheology. Numerous experimental anomalies provide insight into suspension behavior.     

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     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

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.     

The Conjugacy of Stochastic and Random Differential Equations and the Existence of Global Attractors

We consider stochastic differential equations in d-dimensional Euclidean space driven by an m-dimensional Wiener process, determined by the drift vector field f₀ and the diffusion vector fields f₁,...,f_m, and investigate the existence of global random attractors for the associated flows . For this purpose is decomposed into a stationary diffeomorphism given by the stochastic differential equation on the space of smooth flows on R^d driven by m independent stationary Ornstein Uhlenbeck processes z¹,...,z^m and the vector fields f₁,...,f_m, and a flow generated by the nonautonomous ordinary differential equation given by the vector field (_t/x)^–1[f₀(_t)+ _i=1 ¹ f_i(_t)z _t ⁱ ]. In this setting, attractors of are canonically related with attractors of . For , the problem of existence of attractors is then considered as a perturbation problem. Conditions on the vector fields are derived under which a Lyapunov function for the deterministic differential equation determined by the vector field f₀ is still a Lyapunov function for , yielding an attractor this way. The criterion is finally tested in various prominent examples.     

The asymptotic theory of semilinear perturbed telegraph equation and its application

I.IntroductionConsiderthefollowingsemilinearperturbedtelegraphequationuII-u., P'u==sj(t,x,u,ul,u,,e)(-ooo)(l.l)u(o,x)=u,(x,e)(-ooo,u=u(t,x),fuoandulsatisfycertainsmoo…     

Mixed convection from a horizontal plate to fluids of any Prandtl number

Laminar mixed convection over a horizontal plate with uniform wall temperature or uniform wall heat flux is analyzed by introducing proper buoyancy parameters and transformation variables for fluids of any Prandtl number between 0.001 and 10,000. Both cases of buoyancy assisting and opposing flow conditions are investigated. For the buoyancy-assisting case, the obtained numerical results are very accurate over the entire range of mixed convection intensity from pure forced convection limit to pure free convection limit. For the buoyancy-opposing case, solutions are obtained from the forced convection limit to the point of breakdown.

---

Mischkonvektion an einer horizontalen Platte für Fluide mit beliebiger Prandtl-Zahl

Zusammenfassung Es wurde laminare Mischkonvektion an einer horizontalen Platte mit einheitlicher Wandtemperatur oder einheitlicher Wandwärmestromdichte bei Einführung zweckmäßiger Auftriebsparameter und Transformationsvariablen für Fluide mit beliebiger Prandtl-Zahl zwischen 0,001 und 10 000 untersucht. Es wurden die Fälle der Strömung entgegen und in Richtung der Auftriebskraft untersucht. Für den Fall der Strömung in Richtung der Auftriebskraft wurden sehr genaue numerische Ergebnisse für den gesamten Bereich der gemischten Konvektion von rein erzwungener Konvektion bis zu rein freier Konvektion erhalten. Für den Fall der Strömung entgegen der Auftriebsrichtung wurden Lösungen für erzwungene Konvektion bis zum Umkehrpunkt erhalten.

---

Nomenclature C _f local friction coefficient - f reduced stream function - g gravitational acceleration - Gr local Grashof number for UWT,g (T _w –T )x ³/ ² - Gr* local Grashof number for UHF,g q _w x ⁴/k ² - m =10 for UWT; and =6 for UHF - n =5 for UWT; and =3 for UHF - Nu local Nusselt number - p pressure - Pr Prandtl number,/ - q _w wall heat flux - Ra local Rayleigh number for UWT,Gr Pr - Ra* local Rayleigh number for UHF,Gr*Pr - Re local Reynolds number,u x/ - T fluid temperature - T _w wall temperature - T free-stream temperature - u velocity component inx-direction - u free-stream velocity - v velocity component iny-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - =0 for UWT; and =1 for UHF - buoyancy parameter, =( Ra)^1/5/( Re)^1/2 for UWT; and =( Ra*)^1/6/( Re)^1/2 for UHF - pseudo-similarity variable, (y/x) - dimensionless temperature, =(T–T )/(T _w –T ) for UWT; and =(T–T )/(q _w x/k) for UHF - =[( Re)^1/2+( Ra)^1/5] for UWT; and =[( Re)^1/2+( Ra*)^1/6] for UHF - dynamic viscosity - kinematic viscosity - /(1+) - dimensionless pressure - density - Pr/(1+Pr) - _w wall shear stress,(u/y) _y=0 - stream function - Pr/(1+Pr)^1/3     

Similar solutions and the asymptotics of filtration equations

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