A Strongly Degenerate Quasilinear Equation: the Parabolic Case

We prove the existence and uniqueness of entropy solutions of the Neumann problem for the quasilinear parabolic equation u_t=÷ a(u, Du), where a(z,)=f(z,), and f is a convex function of with linear growth as ||||, satisfying other additional assumptions. In particular, this class includes the case where f(z,)=(z)(), >0, and is a convex function with linear growth as ||||.     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

Filament stretching equation in Lagrangian coordinates

The general integral of the very simple equation ^21/n/() was found to describe the cross sectional area of filaments of isothermal power law fluids while in transient stretching where is time and is the initial location of fluid molecules at time = 0 given as the distance from a reference point fixed in space. Any such stretching transient given as a solution of the above equation is physically realizable subject to the restrictions > 0 and/ < 0.     

Analysis of slip and temperature jump coefficients in a binary gas mixture

The possibility of simplifying the formulas obtained by the Maxwell-Loyalka method for the velocity _u, temperature _T and diffusion _d slip coefficients and the temperature jump coefficient in a binary gas mixture with frozen internal degrees of freedom of the molecules is considered. Special attention is paid to gases not having sharply different physicochemical properties. The formulas are written in a form convenient for use without linearization in the thermal diffusion coefficient. They are systematically analyzed for mixtures of inert gases, N₂, O₂, CO₂, and H₂ at temperatures extending from room temperature to 2500°K. It is shown that for the molecular weight ratios m^* = m₂/m₁ considered the expressions for _u and can be radically simplified. With an error acceptable for practical purposes (up to 10%) it is possible to employ expressions of the same structural form as for a single-component gas: for _u if 1 m^* 6, and for if 1 m^* 3. When 1 m^* 2 the expression for _T can be simplified with a maximum error of 5%. Within the limits of accuracy of the method the expression for _t can be linearized in the thermal diffusion coefficient. An approximate expression convenient for practical calculations is proposed for _d Finally, the , _u, and _T for a single-component polyatomic gas with easy excitation of the internal degrees of freedom of the molecules are similarly analyzed; it is shown that these expressions can be considerably simplified.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 152–159, November–December, 1990.     

Mass transfer with chemical reaction in a laminar falling film

Zusammenfassung Es wird eine analytische Lösung für die Absorption in einem laminaren Rieselfilm mit homogener und heterogener chemischer Reaktion 1. Ordnung vorgestellt, wobei der Stofftransportwiderstand auf der Gasseite liegt. Die Lösung ist eine Funktion von drei dimensionslosen ParameternBi, und, welche die BiotZahl und einen homogenen bzw. heterogenen Reaktionsparameter darstellen. Es wird gezeigt, daß für feste Werte vonBi und die Absorptionsrate (bezogen auf die Breite 1 des Rieselfilms) über eine gewisse Länge (dimensionslos) des Rieselfilms unabhängig von ist, wenn, < 0,6 ist. Die laufende Länge wird von der Stelle aus gemessen, an der die Absorption beginnt. Für b 0,6 nimmt der FlußQ mit zu, erreicht aber einen Sättigungswert bei=10, wonachQ nurmehr sehr langsam anwächst. Jedoch für ein gegebenes und ohne Übergangswiderstand im Film (Bi ) nimmtQ mit für alle 0 zu.

---

Mass transfer with chemical reaction in a laminar falling film

An analytical solution is presented for gas absorption in a laminar falling film with first-order homogeneous and heterogeneous chemical reaction and external gas-phase mass transfer resistance. The solution depends on three dimensionless parametersBi, and, wich represent the Biot number, homogeneous and heterogeneous reaction parameters, respectively. It is shown that for fixed values ofBi and, the rate of gas absorption (per unit breadth) over a certain length; (dimensionless) along the falling film measured from the point where surface absorption begins is independent of if < 0.6. For 0.6, this fluxQ increases with but reaches a saturation value at=10 beyond whichQ increases very slowly. But for given and zero gas film resistance (Bi ),Q increases with for all 0.

     

Flow of viscoelastic fluid between confocal elliptic cylinders

Summary Using elliptic coordinates, the flow pattern of a fluid of grade four between two elliptic tubes is determined. A comparison between the position of the maximum of the axial velocity in the present case and in the case of two concentric circular tubes shows a basic difference. In the elliptic case the maximum is shifted towards the external wall, while in the case of concentric circular tubes the shift is in the direction of the internal wall. The secondary flow shows dissymmetry with reference to the intermediate line , which itself lies nearer to the external wall.

---

Zusammenfassung Unter Benutzung elliptischer Koordinaten wird die Strömung zwischen zwei elliptischen Rohren bestimmt. Ein Vergleich zwischen der Lage des axialen Geschwindigkeitsmaximums im vorliegenden Fall und im Fall zweier konzentrischer Kreisrohre ergibt einen grundsätzlichen Unterschied: Das Maximum ist im elliptischen Fall zur äußeren Wand hin verschoben, während die Verschiebung im Fall der konzentrischen Kreisrohre zur inneren Wand hin erfolgt. Die Sekundärströmung ist unsymmetrisch relativ zur mittleren Stromlinie , die selbst näher zur äußeren Wand liegt.

---

A planar domain representing the annular region - vector inx ₁,x ₂-plane - x _i rectangular coordinates - rectangular unit vectors - , elliptic coordinates - ₁, ₂ ellipses representing respectively the internal and external tubes - = ₂ – ₁ annular widthy = ( – ₁)/ - µ 1^st grade material constant - _i 2^nd grade material constants - _i 3^rd grade material constants - _i 4^th grade material constants - I unit tensor - T _E extra stress (T + pI) - V potential of body forces - material density = (p/) + V = –ax ₃ + () - a specific driving force - arbitrary scalar function - A _k Rivlin-Eriksen tensors - S stress scalar defined onA - t stress vector defined onA - P stress tensor defined onA - v axial velocity - v _i i ^th term in the approximation ofv - u velocity vector perpendicular to the axis 4( ₃ + ₄ + ₅ + 1/2₆) –2/µ(2 ₁ + ₂)( ₂ + ₃) - T stress tensor - p arbitrary hydrostatic pressure - u _i i ^th term in the approximation ofu - stream function definingu - _i i ^th term in the approximation of With 8 figures and 1 table     

Elastohydrodynamic lubrication in a plane slider bearing

Summary The present work deals with the case of a two-dimensional slider bearing with a rigid pad and an elastic bearing. Fluid viscosity is assumed to be only a pressure function. We determined the bearing deformation, the pressure distribution and the load capacity at different values of the inclination angle of the slider, with a numerical integration of the system consisting of the elasticity and Reynolds equations. The results show that, with an iso-viscous fluid, bearing elasticity causes a load capacity decrease. Instead bearing elasticity together with the variation of fluid viscosity due to pressure causes a load capacity greater than that of the iso-viscous case (=0).

---

Sommario Il presente lavoro studia il problema della coppia prismatica lubrificata con pattino rigido di allungamento infinito e cuscinetto deformabile; si suppone che la viscosità del fluido sia funzione della sola pressione. Il sistema di equazioni, costituito dall'equazione di Reynolds e dall'equazione dell'elasticità, è stato risolto numericamente, determinando la deformazione del cuscinetto, andamento della pressione e la capacità di carico per diversi valori dell'inclinazione del pattino. I risultati dimostrano che, con fluido isoviscoso, la deformabilità del cuscinetto determina una riduzione della capacità di carico. Se si considera, invece, effetto combinato dell'elasticità del cuscinetto e della variazione della viscosità del fluido, la capacità di carico risulta maggiore di quella che si ottiene con fluido isoviscoso (=0).

---

Nomenclature /L - /L - x/L - x/L - - ¯C CZ/h ₁ - E elasticity modulus - h film thickness - H elastic deformation of the bearing - h ₁ minimum film thickness - h ₂ inlet thickness - inclination of the pad - h Z/h ₁ - HZ/h ₁ - L pad length - viscosity - ₀ viscosity with no over-pressure - p over pressure - p P _ec-P _rc where:ec=elastic caserc=rigid case - P h ₁ ² /6₀VL - h ₂/h ₁=1+L/h ₁ - FV bearing velocity - W load capacity per unit width - Wh ₂ ¹ /6₀ VL ² - Z E h ₃ ¹ /12 ₀ VL ² A first version of this paper was presented at the 7th National AIMETA congress, held at Trieste, October 2–5, 1984. This work was supported by C.N.R.     

Abkühlungscharakteristik einer gut wärmeleitenden Platte,die bei der Bewegung durch ein Fluid gestreckt wird

Zusammenfassung Das Abkühlungsverhalten einer gut wärmeleitenden Platte, die durch ein Fluid bewegt und dabei nach einem beliebigen Potenzgesetz gestreckt wird, wurde untersucht. Der Impulstransport wird dabei durch eine Gleichung beschrieben, die der Falkner-Skan-Gleichung mit einem Parameter für die Dehnung ähnelt. Für das thermische Problem wurden bezüglich des Ortes ähnliche Lösungen gefunden, die sowohl für kleine als auch große Werte des Fließabstandes gültig sind. Die Konvergenz der Reihen für kleine-Werte konnte durch die Euler-Transformation verbessert werden, so daß die Ergebnisse für alle-Werte Gültigkeit haben.

---

The cooling of a low heat resistance stretching sheet moving through a fluid

The cooling characteristics of a low-heat-resistance sheet moving through a fluid have been studied for arbitrary power law stretching. The momentum transfer is governed by an equation that resembles the Falkner-Skan equation with stretching as a parameter. For the thermal problem, locally similar solutions have been obtained for small and large values of the streamwise distance. The convergence of low series has been improved by Euler transformation to yield results valid for all values of.

     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

Motion of large gas bubbles ascending in a liquid

An equation is derived for the ascent velocity of large gas bubbles in a liquid. This velocity is assumed to be governed by the propagation of a wavelike perturbation caused by the bubble in the liquid.Notation w bubble (or drop) velocity - specific gravity - dynamic viscosity - kinematic viscosity - r bubble (or drop) radius - surface tension - coefficient of friction - g gravitational acceleration - D bubble (or drop) diameter - p pressure - c propagation velocity of the wavelike perturbation - wavelength     

Conditions of simulation of stagnation point heat transfer from a high-enthalpy flow

The problem of local simulation of stagnation point heat transfer to a blunt body is solved within the framework of boundary layer theory on the assumption that the simulation subsonic high-enthalpy flow is in equilibrium outside the boundary layer on the model, while the parameters of the natural flow are in equilibrium at the outer edge of the boundary layer on the body. The parameters of the simulating subsonic flow are expressed in terms of the total enthalpyH ₀, the stagnation point pressurep _w and the velocityV ₁ for the natural free-stream flow in the form of universal functions of the dimensionless modeling coefficients=R _m ^* /R _b ^* ( .<1),=V ₁/2H ₀ ( .<1) whereR _m ^* and R _b ^* are the effective radii of the model and the body at their stagnation points. Approximate conditions for modeling the heat transfer from a high-enthalpy (including hypersonic) flow to the stagnation point on a blunt body by means of hyposonic (M1) flows, corresponding to the case ²1, are obtained. The possibilities of complete local simulation of hypersonic nonequilibrium heat transfer to the stagnation point on a blunt body in the hyposonic dissociated air jets of a VGU-2 100-kilowatt induction plasma generator [4, 5] are analyzed.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.1, pp. 172–180, January–February, 1993.     

Discrete time semigroup transformations with random perturbations

Let (X, ) and (Y,C) be two measurable spaces withX being a linear space. A system is determined by two functionsf(X): X X and:X×YX, a (small) positive parameter and a homogeneous Markov chain {y _n } in (Y,C) which describes random perturbations. States of the system, say {x _n X, n=0, 1,}, are determined by the iteration relations:x _n+1 =f(x _n )+(x _n ,Y_n+1) forn0, wherex ₀ =x ₀ is given. Here we study the asymptotic behavior of the solutionx _n as 0 andn under various assumptions on the data. General results are applied to some problems in epidemics, genetics and demographics.Supported in part by NSF Grant DMS92-06677.Supported in part by NSF Grant DMS93-12255.     

Heat transport in a parabolic tube of circular cross-section

The local temperature has been determined for a viscous liquid flowing through a paraboloidal tube. Wall temperature and inlet temperature have been considered constant. The liquid flow was considered as creeping flow and its velocity distribution was determined by solving the biharmonic differential equation of the stream function. The local temperature was evaluated numerically from the analytical results.

---

Wärmetransport im Paraboloidrohr

Zusammenfassung Es wird die lokale Temperatur in einer viskosen Strömung durch ein Paraboloidrohr bestimmt. Dabei wird konstante Wand- und Einlauftemperatur angenommen. Die kriechende Strömungsgeschwindigkeit wurde aus der Lösung der biharmonischen Differentialgleichung der Stromfunktion bestimmt. Die lokale Temperatur wurde aus den analytischen Ergebnissen für einige Paraboloidrohre numerisch bestimmt.

---

Nomenclature ₁ F ₁ confluent hypergeometric function - diffusivity - T(, , ) temperature - T _w temperature at the paraboloidal wall - T _i temperature at the inlet - u(, ) flow velocity of viscous liquid in -direction - volumetric flow - eigenvalues of confluent hypergeometric function - streamfunction - _o wall of paraboloidal tube - _o inlet of paraboloidal tube - , , paraboloidal coordinates     

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(, )     

Analysis of suddenly started laminar flow in the entrance region of a circular tube

Suddenly started laminar flow in the entrance region of a circular tube, with constant inlet velocity, is investigated analytically by using integral momentum approach. A closed form solution to the integral momentum equation is obtained by the method of characteristics to determine boundary layer thickness, entrance length, velocity profile, and pressure gradient.Nomenclature M(, , ) a function - N(, , ) a function - p pressure - p* p/1/2U ², dimensionless pressure - Q(, , ) a function - R radius of the tube - r radial distance - Re 2RU/, Reynolds number - t time - U inlet velocity, constant for all time, uniform over the cross section - u velocity in the boundary layer - u* u/U, dimensionless velocity - u ₁ velocity in the inviscid core - x axial distance - y distance perpendicular to the axis of the tube - y* y/R, dimensionless distance perpendicular to the axis - boundary layer thickness - * displacement thickness - /R, dimensionless boundary layer thickness - momentum thickness - absolute viscosity of the fluid - /, kinematic viscosity of the fluid - x/(R Re), dimensionless axial distance - density of the fluid - tU/(R Re), dimensionless time - _w wall shear stress     

Uniaxial extension of an isothermal fluid cylinder in the presence of inertia,surface tension and gravity

Effects of inertia, surface tension and gravity in the constant force stretching of isothermal cylindrical filaments of Newtonian, power-law and Maxwell-type fluids were analysed in Lagrangian coordinates. Solution for the purely gravitational extension of Newtonian fluid cylinder was found to be as simple as = 1 – C ₃ (1 – ) where designates the cross sectional area, the Lagrangian distance and the time. Analytical solutions were also available for the case of inertialess Newtonian and power-law fluids.A first-order backward differencing scheme and minimal computer time were sufficient to numerically analyse the constant force extension of Maxwell-type fluids in the presence of inertia, gravity and surface tension. Effects of inertia, surface tension and gravity on the severity of neck down occurring at either end of the filament are summarized in diagrams. The present approach is valid on any other constitutive model as far as there is a numerical scheme to analyse thehomogeneous extension of a cylinder of that particular fluid.     

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.     

Analysis of melt spinning transients in Lagrangian coordinates

Transients in melt spinning of isothermal power law and Newtonian fluids were found to be governed by an extremely simple partial differential equation ² ( ^1/n )/() = 0 in Lagrangian coordinates where is the cross-sectional area,n the power law exponent, the time and the the time at which a fluid molecule constituting the spinline left the spinneret. The general integral ^1/n =f() +g () of the above governing equation containing two arbitrary functions represents physically attainable spinline transients. Hitherto unknown analytical transient solutions of the above governing equation were obtained for the response of isothermal constant tension spinlines to a stepwise change in tension, spinneret hole area, extrusion speed or extrusion viscosity and for the starting transient in gravitational spinning. Linearized perturbation solutions and the stability limit of the spinline derived from the above new found nonlinear solutions were in agreement with previous findings and the above nonlinear response of the spinline to a step increase in the spinneret hole area was found to be equivalent to Orowan's tandem cylinder model of dent growth in filament stretching.     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

Mass transport in a converging-diverging nozzle

The local concentration has been determined for a viscous liquid flowing through a converging-diverging tube with constant wall- and initial concentration. The liquid flow was considered as creeping flow and its velocity distribution was determined by solving the biharmonic differential equation of the stream function. The mass transport was presented in form of an infinite series of Legendre functions, which rendered with the Galerkin condition a determinant of finite order for the determination of the eigenvalues. The local concentration was evaluated numerically for differently diverging tubes, of which the case of mass transport for viscous (linearized) flow through an orifice of circular cross-section presented a special case.

---

Stofftransport in einem konvergierenden-divergierenden Rohr

Znsammenfassung Es wird die lokale Konzentration eines in einer viskosen Strömung durch ein Venturirohr befindlichen Stoffes bei konstanter Wand-und Einlaßkonzentration bestimmt. Die viskose (kriechende) Strömungsgeschwindigkeit wurde aus der Lösung der biharmonischen Differentialgleichung der Stromfunktion bestimmt. Die Bestimmung des Massentransportes im Rohr wird mit einer Lösung in Form einer unendlichen Reihe von Legendrefunktionen bewältigt, wobei die Eigenwertgleichung eine Determinante unendlicher Ordnung ergibt. Die lokale Konzentration wurde für verschieden divergierende Rohre numerisch ausgewertet, wobei der Massentransport bei Strömung durch eine kreisförmige Öffnung als Sonderfall auftrat.

---

Nomenclature a radius of throat area - c (, ) concentration - c _w wall concentration - c _i initial concentration - D diffusion coefficient - P _n ⁰ Legendre function - u (, ) flow velocity of liquid (in-direction) - flow volume per time unit - _n roots ofP ⁰ ( ₀)=0 - stream function - , , oblate spheroidal coordinates - ₀ wall of the nozzle