Similarity Solutions for Thawing Processes with a Heat Flux Condition at the Fixed Boundary

Similarity solutions for a mathematical model for thawing in a saturated semi-infinite porous medium is considered when change of phase induces a density jump and a heat flux condition of the type is imposed on the fixed face x=0. Different cases depending on physical parameters are analysed and the explicit solution is obtained if and only if an inequality for the thermal coefficient q ₀ is verified. An improvement for the existence of a similarity solution for the same free boundary problem with a constant temperature on the fixed face x=0 is also obtained. Sommario. Vengono considerate soluzioni di similarità per un modello matematico di disgelo di un mezzo poroso saturo semi-infinito allorquando il cambiamento di fase induce un salto di densità ed una condizione di flusso di calore del tipo viene imposta sulla faccia fissa x=0. Si analizzano differenti casi dipendenti da parametri fisici e la soluzione esplicita viene ottenuta se e solo se risulta verificata una diseguaglianzo per il coefficiente termico q ₀. Si ottiene altresi un miglioramento della condizione di esistenza di una soluzione di similarità per lo stesso problema al contorno libero con temperatura costante sulla faccia fissa x=0.     

On the extended self similarity and the form function

In this work some investigations on the properties of the so calledform function which characterizes the scaling behavior of the small scales fluctuations in a turbulent flow are presented. The present analysis is based on previous experimental measurements in homogeneous and non-homogeneous grid-generated turbulence at low Re. The universality properties of the form function are investigated in the frame of the Extended Self Similarity (ESS) form of scaling.

---

Sommario Nel presente lavoro viene presentato uno studio riguardante le proprietà della cosiddetta funzione di forma che caratterizza le leggi di scala delle fluttuazioni di velocità in un flusso turbolento. L'analisi é basata su precedenti misure sperimentali effettuate in flussi turbolenti generati da griglie a bassi Re ed in condizioni omogenee e non omogenee. Le proprietá di universalitá della funzione di forma sono studiate nell'ambito della cosiddetta Extended Self-Similarity.

     

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.     

Stochastic response of structures with small geometric imperfections

Summary A probabilistic model of the geometric imperfections of a real structure is proposed, in order to provide a general theory of the stochastic response of structures in presence of small random deviations from the perfect scheme. The main statistical measures of the stochastic response are derived and an application to the study of a particular conservative elastic system is developed.

---

Sommario Si propone una teoria generale della risposta probabilistica di strutture, in presenza di piccole deviazioni aleatorie dei dati iniziali rispetto allo schema geometrico perfetto. Si deducono le principali proprietà statistiche della risposta della struttura a sollecitazioni esterne deterministiche, e si sviluppa una applicazione riguardante il comportamento aleatorio di un particolare sistema elastico conservativo.

---

List of symbols element of the sample space of events - _kn random variables modelling the structural imperfections - P(o) probability density of random variables - random imperfection of the unloaded structure - u additional displacement of the loaded structure - u_o deterministic fundamental solution for the perfect structure - difference between the additional displacement of the loaded structure and the deterministic fundamental solution for the perfect structure - V₁=u₁ buckling mode of the perfect structure - _i intrinsic coordinates of the structure - suitable measure of the magnitude of the random imperfections - scalar geometric variable representing the internal product - random imperfection divided by - single scalar variable denoting the magnitude of the prescribed loads - potential energy of the structure - potential energy of the perfect structure - difference between and - _c lowest critical load - _s real local maximum for the magnitude of the prescribed loads - _c divided by _S - E{} expected value of a random variable - ² variance of a random variable - , random variables defined by Eq. (21)     

A fourfold generalization of Peaucellier's inversion cell

New inversors are proposed which are generalizations of the well-known inversor of Peaucellier. It appears that a kite in the Peaucellier cell is replaceable by an arbitrary 4-bar linkage (abhk) whereas the direction and length of the straight line, produced by the inversor, can be manipulated through the particular choice of the relative polar coordinates of a vertex A of the triangular input link. Formulas are derived for practical inversors with a revolving input link. The ones selected are basically governed by the choice of two transmission angles, ₁ and ₃, by the length L of the acquired line, as well as by its direction represented by the angle /2- comprised between the line L and the frame.

---

Sommario Vengono proposti nuovi inversori quali generalizzazioni del ben noto inversore di Peaucellier. Si mostra che un quadrilatero isoscele nella cella di Peaucellier è sostituibile da un'arbitraria connessione a quattro barre, mentre la direzione e la lunghezza della retta, prodotta dall'inversore, possono essere manipolate con la particolare scelta delle relative coordinate polari di un vertice A del collegamento triangolare di input. Vengono derivate formule per inversori funzionali con un collegamento di input rotante. Quelli solezionati sono principalmente governati dalla scelta di due angoli di trasmissione, ₁ e ₃, dalla lunghezza L della linea ottenuta, così come dalla sua direzione rappresentata dall' angolo /2- compreso tra la linea L ed il riferimento.

     

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     

Second order gravitational field theories with hilbert gauge

Summary We consider, in the field-theoretical approach, a class of gravitational theories deducible by a variational principle in the unrenormalized pseudo-Euclidean space-time. At first order in the coupling constant f we require the theories to coincide with the Einstein one. Moreover we assume the Hilbert gauge which assure the exclusion of the vector component of the gravitational potential . To get the higher order consistency we substitute the most general energy-momentum tensor for the particle tensorT ^(p) in the field equations. Requiring the latter to be deducible by a variational principle varying the potentials , we get a Lagrangian which, varying the particle coordinates, gives the equations of motion. So we get a class of theories depending on 5 arbitrary parameters. To have observable quantities we have to renormalize. So we realize that, to satisfy the equivalence principle, we have to put one of the arbitrary parameters equal to zero. With this choice the class of theories coincides at second order with general relativity.

---

Sommario Si vuole ottenere una classe di teorie gravitazionali deducibili da un principio variazionale, nell'ambito della teoria dei campi e nello spazio-tempo pseudoeuclideo non-rinormalizzato. Si richiede che tali teorie coincidano, al primo ordine nella costante di accoppiamento f, con la teoria di Einstein. Si assume inoltre la gauge di Hilbert al fine di escludere la presenza della componente vettoriale del potenziale . Per ottenere la consistenza al secondo ordine delle equazioni di campo, si sostituisce, in queste ultime, al tensore della particellaT ^(p) il più generale tensore energia-quantità-di-moto . Imponendo alle equazioni di campo di essere deducibili mediante un principio variazionale ove si varino i potenziali , si ottiene una lagrangiana che, ove si varino le coordinate della particella di prova, dà le equazioni di moto. In tal modo si ottiene una classe di teorie dipendenti da 5 parametri arbitrari. Per un confronto con i dati sperimentali è necessario rinormalizzare, onde esprimere quantità osservabili. Si dimostra così che per soddisfare il principio di equivalenza al secondo ordine è necessario porre uno dei 5 parametri uguale a zero e che, con tale scelta, l'intera classe di teorie coincide, al secondo ordine, con la relatività generale.

---

---

Research sponsored by the CNR, Gruppi di ricerca Matematica     

Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

Instability of Buckley-Leverett flow in a heterogeneous medium

We study the simultaneous one-dimensional flow of water and oil in a heterogeneous medium modelled by the Buckley-Leverett equation. It is shown both by analytical solutions and by numerical experiments that this hyperbolic model is unstable in the following sense: Perturbations in physical parameters in a tiny region of the reservoir may lead to a totally different picture of the flow. This means that simulation results obtained by solving the hyperbolic Buckley-Leverett equation may be unreliable.Symbols and Notation f fractional flow function varying withs andx - value off outsideI - value off insideI - local approximation off around¯x - f ^–,f ⁺ values of - f _j ⁿ value off atS _j ⁿ andx _j - g acceleration due to gravity [ms^–2] - I interval containing a low permeable rock - k dimensionless absolute permeability - k ^* absolute permeability [m²] - k _c ^* characteristic absolute permeability [m²] - k _ro relative oil permeability - k _rw relative water permeability - L ^* characteristic length [m] - L ₁ the space of absolutely integrable functions - L the space of bounded functions - P _c dimensionless capillary pressure function - P _c ^* capillary pressure function [Pa] - P _c ^* characteristic pressure [Pa] - S similarity solution - S _j ⁿ numerical approximation tos(x_j, t_n) - S ₁, S₂,S ₃ constant values ofs - s water saturation - value ofs at - s ^L left state ofs (wrt. ) - s ^R right state ofs (wrt. ) - s s for a fixed value of in Section 3 - T value oft - t dimensionless time coordinate - t ^* time coordinate [s] - t _c ^* characteristic time [s] - t _n temporal grid point,t _n=n t - v ^* total filtration (Darcy) velocity [ms^–1] - W, , v dimensionless numbers defined by Equations (4), (5) and (6) - x dimensionless spatial coordinate [m] - x ^* spatial coordinate [m] - x _j spatial grid piont,x _j=j x - discontinuity curve in (x, t) space - right limiting value of¯x - left limiting value of¯x - angle between flow direction and horizontal direction - t temporal grid spacing - x spatial grid spacing - length ofI - parameter measuring the capillary effects - argument ofS - _o dimensionless dynamic oil viscosity - _w dimensionless dynamic water viscosity - _c ^* characteristic viscosity [kg m^–1s^–1] - _o ^* dynamic oil viscosity [kg m^–1s^–1] - _w ^* dynamic water viscosity [k gm^–1s^–1] - _o dimensionless density of oil - _w dimensionless density of water - _c ^* characteristic density [kgm^–3] - _o ^* density of oil [kgm^–3] - _w ^* density of water [kgm^–3] - porosity - dimensionless diffusion function varying withs andx - ^* dimensionless function varying with s andx ^* [kg^–1m³s] - _j ⁿ value of atS _j ⁿ andx _j This research has been supported by VISTA, a research cooperation between the Norwegian Academy of Science and Letters and Den norske stats oljeselskap a.s. (Statoil).     

Nonlinear Taylor stability of viscoelastic fluids

Using Stuart's energy method, the torque on the inner cylinder, for a second order fluid, in the supercritical regime is calculated. It is found that when the second normal stress difference is negative, the flow is more stable than for a Newtonian fluid and the torque is reduced. If the second normal stress difference is positive, then the flow is more stable and there is no torque reduction. Experimental data related to the present work are discussed.Nomenclature a amplitude of the fundamentals - A _ij ⁽¹⁾ , A _ij ⁽²⁾ first and second Rivlin-Ericksen tensors - d r ₂–r ₁ - D d/dx - E - F - g _ij metric tensor - G torque on the inner cylinder in the supercritical regime - h height of the cylinders - k ₀ /d ² - k ₁ /d ² - I ₁ - I ₂ - I ₃ - I ₄ - r ₁, r ₂ radii of inner and outer cylinders respectively - r ₀ 1/2(r ₁+r ₂) - R Reynolds number ₁ r ₁ d/ ₀ - R _c critical Reynolds number - T Taylor number r _{1 1} ² d ^{3 2}/ ₀ ² *) - T _c critical Taylor number - u ₁, v ₁, w ₁ Fundamentals of the disturbance - u _i , v _i , w _i , (i>1) harmonics - mean velocity (not laminar velocity) - u –u ₁/ar _{1 1} - v v ₁/Rar _{1 1} - x (r–r ₀)/d - , material constants - 0 viscosity - wave number d - density - ₁ angular velocity of inner cylinder - tilde denotes complex conjugate     

Darcy-Forchheimer natural,forced and mixed convection heat transfer in non-Newtonian power-law fluid-saturated porous media

The governing equation for Darcy-Forchheimer flow of non-Newtonian inelastic power-law fluid through porous media has been derived from first principles. Using this equation, the problem of Darcy-Forchheimer natural, forced, and mixed convection within the porous media saturated with a power-law fluid has been solved using the approximate integral method. It is observed that a similarity solution exists specifically for only the case of an isothermal vertical flat plate embedded in the porous media. The results based on the approximate method, when compared with existing exact solutions show an agreement of within a maximum error bound of 2.5%.Nomenclature A cross-sectional area - b _i coefficient in the chosen temperature profile - B ₁ coefficient in the profile for the dimensionless boundary layer thickness - C coefficient in the modified Forchheimer term for power-law fluids - C ₁ coefficient in the Oseen approximation which depends essentially on pore geometry - C _i coefficient depending essentially on pore geometry - C _D drag coefficient - C _t coefficient in the expression forK ^* - d particle diameter (for irregular shaped particles, it is characteristic length for average-size particle) - f _p resistance or drag on a single particle - F _R total resistance to flow offered byN particles in the porous media - g acceleration due to gravity - g _x component of the acceleration due to gravity in thex-direction - Grashof number based on permeability for power-law fluids - K intrinsic permeability of the porous media - K ^* modified permeability of the porous media for flow of power-law fluids - l _c characteristic length - m exponent in the gravity field - n power-law index of the inelastic non-Newtonian fluid - N total number of particles - Nux,_D,F local Nusselt number for Darcy forced convection flow - Nux,_D-F,F local Nusselt number for Darcy-Forchheimer forced convection flow - Nux,_D,M local Nusselt number for Darcy mixed convection flow - Nux,_D-F,M local Nusselt number for Darcy-Forchheimer mixed convection flow - Nux,_D,N local Nusselt number for Darcy natural convection flow - Nux,_D-F,N local Nusselt number for Darcy-Forchheimer natural convection flow - pressure - p exponent in the wall temperature variation - _Pe c characteristic Péclet number - _Pe x local Péclet number for forced convection flow - _Pe x modified local Péclet number for mixed convection flow - _Ra c characteristic Rayleigh number - _Ra x local Rayleigh number for Darcy natural convection flow - _Ra x local Rayleigh number for Darcy-Forchheimer natural convection flow - Re convectional Reynolds number for power-law fluids - Reynolds number based on permeability for power-law fluids - T temperature - T _e ambient constant temperature - T w,_ref constant reference wall surface temperature - T _w(X) variable wall surface temperature - T _w temperature difference equal toT w,_ref–T _e - T ₁ term in the Darcy-Forchheimer natural convection regime for Newtonian fluids - T ₂ term in the Darcy-Forchheimer natural convection regime for non-Newtonian fluids (first approximation) - T _N term in the Darcy/Forchheimer natural convection regime for non-Newtonian fluids (second approximation) - u Darcian or superficial velocity - u ₁ dimensionless velocity profile - u _e external forced convection flow velocity - u _s seepage velocity (local average velocity of flow around the particle) - u _w wall slip velocity - U c _M characteristic velocity for mixed convection - U c _N characteristic velocity for natural convection - x, y boundary-layer coordinates - x ₁,y ₁ dimensionless boundary layer coordinates - X coefficient which is a function of flow behaviour indexn for power-law fluids - effective thermal diffusivity of the porous medium - shape factor which takes a value of/4 for spheres - shape factor which takes a value of/6 for spheres - ₀ expansion coefficient of the fluid - _T boundary-layer thickness - T ₁ dimensionless boundary layer thickness - porosity of the medium - similarity variable - dimensionless temperature difference - coefficient which is a function of the geometry of the porous media (it takes a value of 3 for a single sphere in an infinite fluid) - ₀ viscosity of Newtonian fluid - ^* fluid consistency of the inelastic non-Newtonian power-law fluid - constant equal toX(2 ^2–n )/ - density of the fluid - dimensionless wall temperature difference     

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     

Numerical modelling of unbounded domains in coupled fleld problems

Summary Two- and three-field problems are often defined in domains which may be assumed as unbounded. The traditional approach for their numerical simulation, within the framework of the finite element method, is by simple truncation of the mesh at a finite boundary. This fact both results in a large number of degrees of freedom and causes often errors in the analysis, due to the difficulty of setting correct conditions at the finite boundary.This paper shows the possible errors of the ensuing numerical solution and points out the usefulness of the infinite elements to simulate the far field response. Three examples from the field of isothermal and nonisothermal consolidation are presented where the improvements in the numerical simulation obtained by the use of infinite elements are evidenced. These examples may be considered as representative for a series of other coupled problems involving partial differential equations with first order time derivatives.

---

Sommario Problemi di interazione fra due e tre campi sono spesso deflniti in domini che possono essere assunti come illimitati. Il modo tradizionale per la modellazione numerica di questi casi, nell'ambito del metodo degli elementi finiti, è quello di assumere frontiere fittizie in corrispondenza alle quali si devono imporre condizioni al contorno spesso di non facile valutazione. Questo modo di procedere comporta un numero di gradi di libertà assai elevato e può essere fonte di errori derivanti dall'imposizione di condizioni al contorno non corrette in corrispondenza della frontiera fittizia.Nel presente lavoro vengono evidenziati possibili errori delle conseguenti soluzioni numeriche e viene rimarcata l'utililità di elementi infiniti nella trattazione di questi problemi. In tre esempi di consolidazione isoterma e non, vengono messi in luce i miglioramenti della soluzione dovuti all'uso di elementi infiniti. Questi esempi sono rappresentativi di una più vasta classe di problemi accoppiati govemati da equazioni differenziali alle derivate parziali con derivate del primo ordine rispetto al tempo.

---

Notation b body force vector - c coefficient of consolidation - c strain independent vector defining the creep strain rate - C _s specific heat of the solid phase - C _w specific heat of the fluid phase - D _T tangential stiffness matrix - g gravity acceleration - k absolute permeability matrix - k coefficient of thermal diffusivity - k _s bulk modulus of the solid phase - k _w bulk modulus of the fluid - L differential operator which relates displacements to strains - m (1 1 1 0 0 0) ^T - p pore pressure - Q _e volumetric outflow of the fluid per unit volume of the solid - Q _h outflow of heat per unit volume of solid - t time variable - dimensionless time parameter - T temperature increase over an equilibrium state - boundary traction vector - u displacement vector - V ^a apparent velocity of the fluid - z elevation above some datum - _s thermal expansion coefficient of the solid phase - _w thermal expansion coefficient of the fluid - total strain vector of the soil skeleton - ₀ represents all other strains not directly associated with stress changes - thermal conductivity matrix of the soil - dynamic viscosity - _s density of the solid - _w density of the fluid - effective stress in the soil skeleton - porosity Paper presented at the First Italian Meeting on Computational Mechanics (Milan, June 1986).     

A technique for evaluation of three-dimensional behavior in turbulent boundary layers using computer augmented hydrogen bubble-wire flow visualization

A system is described which allows the recreation of the three-dimensional motion and deformation of a single hydrogen bubble time-line in time and space. By digitally interfacing dualview video sequences of a bubble time-line with a computer-aided display system, the Lagrangian motion of the bubble-line can be displayed in any viewing perspective desired. The u and v velocity history of the bubble-line can be rapidly established and displayed for any spanwise location on the recreated pattern. The application of the system to the study of turbulent boundary layer structure in the near-wall region is demonstrated.List of Symbols Reynolds number based on momentum thickness u /v - t⁺ nondimensional time - u shear velocity - u local streamwise velocity, x-direction - u ⁺ nondimensional streamwise velocity - v local normal velocity, -direction - x ⁺ nondimensional coordinate in streamwise direction - ⁺ nondimensional coordinate normal to wall - ₊ ^wire wire nondimensional location of hydrogen bubble-wire normal to wall - z ⁺ nondimensional spanwise coordinate - momentum thickness - v kinematic viscosity - _W wall shear stress     

Zur analytischen Beschreibung des Füllstoffeinflusses auf das Fließverhalten von Kunststoff-Schmelzen

From literature some representative equations have been compiled describing the influence of filler on the viscosity of polymer melts. By application of these on the experimental results obtained from GF-SAN it was found that the relative viscosity _R , i.e. the ratio of the viscosities of the filled and unfilled melt, shows a pronounced dependence on the shear rate but not on the shear stress. Defining _R with constant and not with constant (as it is usually done), an analytical approach is possible independent of Further the influence of pressure, temperature and filler content on the zero-shear viscosity of filled polymer melts may be expressed by a modified Arrhenius equation.

     

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     

Zum Mischen rheologisch inhomogener Stoffsysteme auf Einschneckenmaschinen

Zusammenfassung Das Mischen von Stoffen mit unterschiedlichen rheologischen Eigenschaften in Schneckenmaschinen ist in der Kunststoffauf- und -verarbeitung eine Standardaufgabe. Trotzdem gibt es hierfür kein zufriedenstellendes mathematisch-physikalisches Modell. Daher werden zunächst einfache Mischmodelle diskutiert. Auf der Basis dieser Modelle wird dann unter Berücksichtigung der Besonderheiten des Plastifizierextruderprozesses eine Mischgütebeziehung mathematisch formuliert. Die experimentelle Überprüfung erfolgt mit Hilfe der Grauwertanalyse extrudierter Zweistoffsysteme, bei denen ein Stoff mit Ruß eingefärbt war. Da der Mischprozeß hochgradig stochastisch ist, streuen die Meßergebnisse. Unter Berücksichtigung dieses Tatbestandes ist der theoretische Ansatz zufriedenstellend.

---

Mixing of polymer resins with different rheological properties is a usual demand in plastics processing using screw extruders. A mathematical model describing this processing problem sufficiently is not known, however. Therefore, simple mixing models will be discussed. Based on these, a concept for the calculation of mixing homogeneity will be presented, including the particular requirement of the plasticating screw process. An experimental investigation utilizes the grey-value analysis of extruded two-component materials, which in one phase is carbon-black filled. Considering the fact that the mixing process is highly random, the theoretical model leads to a good level of aggreement with the scattering measurement data.

---

b Schneckenkanalbreite - B Bandbreite der Grauwerte - c Konstante - mittlere Konzentration, bezogen auf die Grauwertbandbreite - h Höhe, Gangtiefe, Schneckenkanalhöhe - h ₀ Gangtiefe der Einzugszone - h ₁ Gangtiefe der Ausstoßzone - L Länge - gemittelte Schmelzebettlänge - n Exponent des Potenzfließgesetzes - s Standardabweichung der Grauwerte bezogen auf die Grauwertbandbreite - S Standardabweichung der Grauwerte - t Verweilzeit - t ₁ kürzeste Verweilzeit - mittlere Verweilzeit - ₀ Umfangsgeschwindigkeit - mittlere Geschwindigkeit - V Volumenstrom - w Dicke eines Kontrollelements - w Ausstreichdicke eines Kontrollelements - x Koordinate - Mittelwert der Grauwerte - y Koordinate - Scherdeformationswinkel - Scherdeformation - mittlere Scherdeformation - Schergeschwindigkeit - Viskosität - ₁ dimensionslose kürzeste Verweilzeit - dimensionsloser Volumenstrom - _LSM laminarer Schermischgrad - _{LSM, the} theoretischer laminarer Schermischgrad - _{LSM, exp} experimenteller laminarer Schermischgrad - ² Varianz der Verweilzeit im Schmelzebett - Schubspannung - Gangsteigungswinkel der Schnecke - ø Volumenanteil - dimensionslose Kennzahl     

The mathematically linear analysis of strain in continua

Summary One-one correspondence is postulated for two coordinate continua. One continuum is regarded as the initially undeformed state of a currently deformed continuum. The two continua are orthogonal trasformations each to the other. The square root of the quadratic metric, when the appropriate self-conjugate stretch dyadic is expressed in its principal form, gives the mathematically linear form to the analysis. The self-conjugate dyadic is expressed as =grad grad P in terms of a scalar potential function P. The physical and mathematical continuity of the strain dyadic is ensured by curl =0.A 4-vector analysis is evolved from a 4-vector quinternion analysis. The 4-vector analysis is of the same form as the usual 3-vector analysis that evolved from Hamilton's quaternion analysis.

---

Sommario Si postula una corrispondenza biunivoca per due continui coordinati. Un continuo viene considerato come lo stato inizialmente non deformato di un continuo attualmente deformato.I due continui sono trasformazioni ortogonali uno dell'altro. La radice quadrata della metrica quadratica, quando la relativa diade di dilatazione autoconiugata venga espressa nella sua forma principale dà la forma matematicamente lineare alla analisi.La diade coniugata viene espresse come =grad grad P in termini di una funzione potenziale scalare P. La continuità fisica e matematica della diade di deformazione è assicurata da rot =0.Una analisi tetravettoriale viene sviluppata da analisi di quinternioni tetravettoriali. L'analisi tetravettoriale è della stessa forma dell'analisi trivettoriale che viene vsiluppata dalla analisi quaternionica trivettoriale di Hamilton.

---

A generalisation that, when particularised, applies to the finite or infinitesimal straining of elastic bodies, incremental straining of plastic or fluidic bodies in 3-dimensional continua and to the space-time continuum as a 4-dimensional continuum.     

Evaluation of the performance characteristics of a thermal transient anemometer

The Goertler instability of a hypersonic boundary layer and its influence on the wall heat transfer are experimentally analyzed. Measurements, made in a wind tunnel by means of a computerized infrared (IR) imaging system, refer to the flow over two-dimensional concave walls. Wall temperature maps (that are interpreted as surface flow visualizations) and spanwise heat transfer fluctuations are presented. Measured vortices wavelengths are correlated to non-dimensional parameters and compared with numerical predictions from the literature.List of symbols c _p Specific heat coefficient at constant pressure of the free stream - F Input (true) image - F ₀ Fourier number - Restored image - G Recorded (degraded) image - G Goertler number based on the boundary layer thickness, as defined by Eq. (3) - H System transfer function - M Mach number - Pr Prandtl number - p ₀ Stagnation pressure - Exchanged net heat flux - Convective heat flux - Radiative heat flux - r Recovery factor - Re _m Unit Reynolds number - Re _x Local Reynolds number based on the distance from the leading edge - Re Local Reynolds number based on the boundary layer thickness - Curvature radius - St Stanton number, as defined by Eq. (7) - T _aw Adiabatic wall temperature - T _w Wall temperature - T ₀ Stagnation temperature - t Time - V Free stream velocity - x Streamwise spatial coordinate - y Normal-to-wall spatial coordinate - z Spanwise spatial coordinate - Thermal diffusivity coefficient - Disturbance wavenumber - Non dimensional wavenumber - Boundary layer thickness - Goertler number based on the vortices wavelength - Vortices wavelength - Free stream density - Disturbance total amplification, as defined by Eq. (3) - Disturbance (spatial) growth rate - Non-dimensional growth rate - Perturbation amplitude of a generic quantity - Perturbation amount     

Effect of nonlinear constitutive equation on vibrations of beams

Summary This paper is concerned with the analytical investigation of static and dynamic nonlinear behaviors of beams with different boundary conditions. While geometric type of nonlinearities on beams have been investigated extensively, material type nonlinearities have received very little attention. Therefore, material nonlinearities of the Ramberg-Osgood type are considered in this analysis. The use of Self-Generating functions for nonlinear beam problems is demonstrated for this type of nonlinearity. Transverse shear and rotatory inertia effects have been included in the formulation to study moderately thick beams. For all the cases investigated here nonlinear frequency ratios are calculated at various amplitudes of vibration and geometric parameters of beams. Numerical results indicate that the Ramberg-Osgood type nonlinearity produces softening-type responses. The study is limited to materials which are nonlinearly elastic and the effect of geometric nonlinearity is not considered in this paper.

---

Sommario Questo lavoro riguarda lo studio dei comportamenti non lineari statici e dinamici di travi con diverse condizioni a contorno. Mentre le non linearità di tipo geometrico sono state studiate estesamente, quelle di tipo non lineare hanno ricevuto un'attenzione molto ridotta. Perciò in questa ultima analisi si considerano non linearità materiali del tipo di Ramberg-Osgood. Si dimostra uso delle funzioni autogeneratrici nei problemi non lineari per le travi con questo tipo di non linearità. Nella formulazione dello studio di travi moderatamente spesse si sono inclusi effetti di taglio e inerzia rotatoria. Per tutti i casi qui studiati si calcolano i rapporti di frequenza non lineare per varie ampiezze di vibrazione e parametri geometrici delle travi. I risultati numerici indicano che la non linearità del tipo di Ramberg-Osgood produce risposte del tipo ammorbidimento. Lo studio si limita a materiali con nonlinearità elastica e non si considera nel lavoro l'effetto della non linearitá geometrica.

---

List of symbols T _s Transverse shear - k ₁ Shear connection factor - A, B, m Material constants - I Moment of inertia - G Shear modulus - Mass per unit length - w Lateral displacement - h Beam thickness - b Bredth of beam - t Time - R _i Rotatory inertia - a Area of cross-section of beam - q(x) Lateral load on beam - _x Stress - _x Strain - Length of beam - x Beam coordinate - - r Radius of gyration - w Nondimensional maximum deflectionw _max/r - q ₀ ^* Nondimensional load, (q ₀ ³/Al - Thickness parameter,h/ - (T _sK₁/Ga) - ₀ Linear frequency - Nonlinear frequency - w _max w measured at the point of maximum deflection.