Stability characteristics of condensate films

The linear stability theory is used to study stability characteristics of laminar condensate film flow down an arbitrarily inclined wall. A critical Reynolds number exists above which disturbances will be amplified. The magnitude of the critical Reynolds number is in all practical situations so small that a laminar gravity-induced condensate film can be expected to be unstable. Several stabilizing effects are acting on the film flow; at an inclined wall these effects are due to surface tension, gravity and condensation mass transfer.

---

Zusammenfassung Mit Hilfe der linearen Stabilitätstheorie werden die Stabilitätseigenschaften laminarer Kondensatfilme an einer geneigten Wand untersucht. Es zeigt sich, daß Kondensatfilme in jedem praktischen Fall ein unstabiles Verhalten aufweisen. Der stabilisierende Einfluß von Oberflächenspannung, Schwerkraft und Stoffübertragung durch Kondensation bewkkt jedoch, daß Störungen in bestimmten Wellenlängenbereichen gedämpft werden.

---

Nomenclature c=c^*/u₀ complex wave velocity, celerity, dimensionless - c^*=c _r ^* + i c _i ^* complex wave velocity, celerity, dimensional - c_p specific heat at constant pressure - g gravitational acceleration - h_fg latent heat - k thermal conductivity of liquid - p^* pressure - p=p^*/u₀ ² dimensionless pressure - Pe=Pr Re^* Peclet number - Pr Prandtl number - Re^*=u₀ / Reynolds number (defined with surface velocity) - S temperature perturbation amplitude - t^* time - t=t^* u₀/ dimensionless time - T temperature - T_s saturation temperature - T_w wall temperature - T=T_s-T_w temperature drop across liquid film - u^*, v^* velocity components - u=u^*/u₀ dimensionless velocity components - v=v^*/u₀ dimensionless velocity components - u₀ surface velocity of undisturbed film flow - v _g ^* vapor velocity - x^*, y^* coordinates - x=x^*/ dimensionless coordinates - y=y^*/ dimensionless coordinates Greek Symbols =^* wave number, dimensionless - ^*=2 /^* wave number dimensional - ^* wave length, dimensional - =^*/ wave length, dimensionless - local thickness of undisturbed condensate film - kinematic viscosity, liquid - density, liquid - _g density vapor - surface tension - = (1 +) film thickness of disturbed film, Fig. 1 - stream function perturbation amplitude - angle of inclination Base flow quantities are denoted by^–, disturbance quantities are denoted by.     

The theory of a vibrating-rod densimeter

The paper presents a theory of a device for the accurate determination of the density of fluids over a wide range of thermodynamic states. The instrument is based upon the measurement of the characteristics of the resonance of a circular section tube, or rod, performing steady, transverse oscillations in the fluid. The theory developed accounts for the fluid motion external to the rod as well as the mechanical motion of the rod and is valid over a defined range of conditions. A complete set of working equations and corrections is obtained for the instrument which, together with the limits of the validity of the theory, prescribe the parameters of a practical design capable of high accuracy.Nomenclature A, B, C, D constants in equation (60) - A _j , B _j constants in equation (18) - a _j ⁺ , a _j ^– wavenumbers given by equation (19) - C _f drag coefficient defined in equation (64) - C _f ^/0 , C _f ^/1 components of C _f in series expansion in powers of - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - force per unit length - F _j ⁺ , F _j ^– constants in equation (24) - f, g functions of defined in equations (56) - G modulus of rigidity - I second moment of area - K constant in equation (90) - k, k constants defined in equations (9) - L half-length of oscillator - Ma Mach number - m _a mass per unit length of fluid a - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equation (17) - P power (energy per cycle) - P _a , P _b power in fluids a and b - p pressure - R radius of rod or outer radius of tube - R _c radius of container - R _i inner radius of tube - r radial coordinate - T tension - T _visc temperature rise due to heat generation by viscous dissipation - t time - v _r , v radial and angular velocity components - y lateral displacement - z axial coordinate - dimensionless tension - _a dimensionless mass of fluid a - _b dimensionless added mass of fluid b - _b dimensionless drag of fluid b - dimensionless parameter associated with - ₀ dimensionless coefficient of internal damping - dimensionless half-width of resonance curve - dimensionless frequency difference defined in equation (87) - spatial resolution of amplitude - R, , , _s , increments in R, , , _s , - dimensionless amplitude of oscillation - dimensionless axial coordinate - ratio of to - _a , _b ratios of to for fluids a and b - angular coordinate - parameter arising from distortion of initially plane cross-sections - _f thermal conductivity of fluid - dimensionless parameter associated with - viscosity of fluid - _a , _b viscosity of fluids a and b - dimensionless displacement - _j jth component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - density of fluid calculated on assumption that * - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - _{rr rr}, _r radial normal and shear stress components - spatial component of defined in equation (13) - _j jth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - phase angle - _r phase difference - _ra , _rb phase difference for fluids a and b - streamfunction - _j jth component defined in equation (22) - dimensionless frequency (based on ) - _a , _b dimensionless frequency in fluids a and b - _s dimensionless frequency (based on _s ) - angular frequency - ₀ resonant frequency in absence of fluid and internal damping - _r resonant frequency in absence of internal fluid - _ra , _rb resonant frequencies in fluids a and b - dimensionless frequency - dimensionless frequency when _a vanishes - dimensionless frequencies when _a vanishes in fluids a and b - dimensionless resonant frequency when _a , _b, _b and ₀ vanish - dimensionless resonant frequency when _a , _b and _b vanish - dimensionless resonant frequency when _b and _b vanish - dimensionless frequencies at which amplitude is half that at resonance     

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

Flow birefringence of polymer melts: Application to the investigation of time dependent rheological properties

Summary Transient stresses including normal stresses, which are developed in a polymer melt by a suddenly imposed constant rate of shear, are investigated by mechanical measurement and, indirectly, with the aid of the flow birefringence technique. For the latter purpose use is made of the so-called stress-optical law, which is carefully checked.It appears that the essentially linear model of the rubberlike liquid, as proposed byLodge, is capable of describing the behaviour of polymer melts rather well, if the applied total shear does not exceed unity. In order to describe also steady state values of the stresses successfully, one should extend measurements to extremely low shear rates.These statements are verified with the aid of a method which was originally designed bySchwarzl andStruik for the practical calculation of interrelations between linear viscoelastic functions. In the present paper dynamic shear moduli are used as reference functions.

---

Zusammenfassung Mit der Zeit anwachsende Spannungen, darunter auch Normalspannungen, wie sie sich nach dem plötzlichen Anlegen einer konstanten Schergeschwindigkeit in einer Polymerschmelze entwickeln, werden mit Hilfe mechanischer Messungen und indirekt mit Hilfe der Strömungsdoppelbrechung untersucht. Für den letzteren Zweck wird das sogenannte spannungsoptische Gesetz herangezogen, dessen Gültigkeit sorgfältig überprüft wird.Es ergibt sich, daß das im Wesen lineare Modell der gummiartigen Flüssigkeit, wie es vonLodge vorgeschlagen wurde, sich recht gut zur Beschreibung des Verhaltens von Polymerschmelzen eignet, solange der im ganzen angelegte Schub den Wert Eins nicht überschreitet. Um auch stationäre Werte der Spannungen in die Beschreibung erfolgreich einzubeziehen, sollte man die Messungen bis zu extrem niedrigen Schergeschwindigkeiten ausdehnen.Die gemachten Feststellungen werden mit Hilfe einer Methode verifiziert, die vonSchwarzl undStruik ursprünglich für die praktische Berechnung von Beziehungen zwischen Zustandsfunktionen entwickelt wurde, die dem linear viskoelastischen Verhalten entsprechen. In der vorliegenden Veröffentlichung dienen die dynamischen Schubmoduln als Bezugsfunktionen.

---

a _T shift factor - B _ij Finger deformation tensor - C stress-optical coefficient, (m²/N) - f (p _jl ) undetermined scalar function - G shear modulus, (N/m²) - G(t) time dependent shear modulus, (N/m²) - G() shear storage modulus, (N/m²) - G() shear loss modulus, (N/m²) - G _r reduced shear storage modulus, (N/m²) - G _r reduced shear loss modulus, (N/m²) - H() shear relaxation time spectrum, (N/m²) - k Boltzmann constant, (Nm/°K) - n _ik refractive index tensor - p undetermined hydrostatic pressure, (N/m²) - p _ij ,p _ik stress tensor, (N/m²) - p ₂₁ shear stress, (N/m²) - p ₁₁ –p ₂₂ first normal stress difference, (N/m²) - p ₂₂ –p ₃₃ second normal stress difference, (N/m²) - q shear rate, (s^–1) - t, t time, (s) - T absolute temperature, (°K) - T ₀ reference temperature, (°K) - x the ratiot/ - x position vector of a material point after deformation, (m) - x position vector of a material point before deformation, (m) - ₀, ₁ constants in eq. [37] - ₀, ₁ constants in eq. [37] - shear deformation - (t, t) time dependent shear deformation - _ij unity tensor - n flow birefringence in the 1–2 plane - (q) non-Newtonian shear viscosity, (N s/m²) - ^* () complex dynamic viscosity, (N s/m²) - | ^* ()| absolute value of complex dynamic viscosity, (N s/m²) - () real part of complex dynamic viscosity, (N s/m²) - () imaginary part of complex dynamic viscosity, (N s/m²) - (t — t) memory function, (N/m² · s) - v number of effective chains per unit of volume, (m^–3) - temperature dependent density, (kg/m³) - ₀ density at reference temperatureT ₀, (kg/m³) - relaxation time, (s) - integration variable, (s) - (x) approximate intensity function - ₁ (x) error function - extinction angle - _m orientation angle of the stress ellipsoid - circular frequency, (s^–1) - 1 direction of flow - 2 direction of the velocity gradient - 3 indifferent direction - t time dependence The present investigation has been carried out under the auspices of the Netherlands Organization for the Advancement of Pure Research (Z. W. O.).North Atlantic Treaty Organization Science Post Doctoral Fellow.Research Fellow, Delft University of Technology.With 11 figures and 2 tables     

On the cross-effects of coupled macroscopic transport equations in porous media

In this paper, the derivation of macroscopic transport equations for this cases of simultaneous heat and water, chemical and water or electrical and water fluxes in porous media is presented. Based on themicro-macro passage using the method of homogenization of periodic structures, it is shown that the resulting macroscopic equations reveal zero-valued cross-coupling effects for the case of heat and water transport as well as chemical and water transport. In the case of electrical and water transport, a nonsymmetrical coupling was found.Notations b mobility - c concentration of a chemical - D rate of deformation tensor - D molecular diffusion coefficient - D _ij ^eff macroscopic (or effective) diffusion tensor - electric field - E ₀ initial electric field - k _ij molecular tensor - j, j ^*, current densities - K _ij macroscopic permeability tensor - l characteristic length of the ERV or the periodic cell - L characteristic macroscopic length - L _ijkl coupled flows coefficients - n _i unit outward vector normal to - p pressure - q _t ,q _t ⁺ , heat fluxes - q _c ,q _c ⁺ , chemical fluxes - s specific entropy or the entropy density - S entropy per unit volume - t time variable - t _ij local tensor - T absolute temperature - v _i velocity - V ₀ initial electric potential - V electric potential - x macroscopic (or slow) space variable - y microscopic (or fast) space variable - _i local vectorial field - _i local vectorial field - electric charge density on the solid surface - , bulk and shear viscosities of the fluid - _ij local tensor - _ij local tensor - _i local vector - _ij molecular conductivity tensor - _ij ^eff effective conductivity tensor - homogenization parameter - fluid density - ₀ ion-conductivity of fluid - _ij dielectric tensor - _i ¹ , _i ² , _i ³ local vectors - ⁴ local scalar - _S solid volume in the periodic cell - _L volume of pores in the periodic cell - boundary between _S and _L - ^s rate of entropy production per unit volume - total volume of the periodic cell - _l volume of pores in the cell On leave from the Politechnika Gdanska; ul. Majakowskiego 11/12, 80-952, Gdask, Poland.     

Impacts of local dispersion and first-order decay on solute transport in randomly heterogeneous porous media

Stochastic subsurface transport theories either disregard local dispersion or take it to be constant. We offer an alternative Eulerian-Lagrangian formalism to account for both local dispersion and first-order mass removal (due to radioactive decay or biodegradation). It rests on a decomposition of the velocityv into a field-scale componentv , which is defined on the scale of measurement support, and a zero mean sub-field-scale componentv _s , which fluctuates randomly on scales smaller than. Without loss of generality, we work formally with unconditional statistics ofv _s and conditional statistics ofv . We then require that, within this (or other selected) working framework,v _s andv be mutually uncorrelated. This holds whenever the correlation scale ofv is large in comparison to that ofv _s . The formalism leads to an integro-differential equation for the conditional mean total concentration c which includes two dispersion terms, one field-scale and one sub-field-scale. It also leads to explicit expressions for conditional second moments of concentration cc. We solve the former, and evaluate the latter, for mildly fluctuatingv by means of an analytical-numerical method developed earlier by Zhang and Neuman. We present results in two-dimensional flow fields of unconditional (prior) mean uniformv . These show that the relative effect of local dispersion on first and second moments of concentration dies out locally as the corresponding dispersion tensor tends to zero. The effect also diminishes with time and source size. Our results thus do not support claims in the literature that local dispersion must always be accounted for, no matter how small it is. First-order decay reduces dispersion. This effect increases with time. However, these concentration moments c and cc of total concentrationc, which are associated with the scale below, cannot be used to estimate the field-scale concentrationc directly. To do so, a spatial average over the field measurement scale is needed. Nevertheless, our numerical results show that differences between the ensemble moments ofc and those ofc are negligible, especially for nonpoint sources, because the ensemble moments ofc are already smooth enough.     

On the prediction of riblet performance with engineering turbulence models

The paper reports the outcome of a numerical study of fully developed flow through a plane channel composed of ribleted surfaces adopting a two-equation turbulence model to describe turbulent mixing. Three families of riblets have been examined: idealized blade-type, V-groove and a novel U-form that, according to computations, achieves a superior performance to that of the commercial V-groove configuration. The maximum drag reduction attained for any particular geometry is broadly in accord with experiment though this optimum occurs for considerably larger riblet heights than measurements indicate. Further explorations bring out a substantial sensitivity in the level of drag reduction to the channel Reynolds number below values of 15 000 as well as to the thickness of the blade riblet. The latter is in accord with the trends of very recent, independent experimental studies.Possible shortcomings in the model of turbulence are discussed particularly with reference to the absence of any turbulence-driven secondary motions when an isotropic turbulent viscosity is adopted. For illustration, results are obtained for the case where a stress transport turbulence model is adopted above the riblet crests, an elaboration that leads to the formation of a plausible secondary motion sweeping high momentum fluid towards the wall close to the riblet and thereby raising momentum transport.Nomenclature c _f Skin friction coefficient - c _f Skin friction coefficient in smooth channel at the same Reynolds number - k Turbulent kinetic energy - K ⁺ k/ _w - h Riblet height - S Riblet width - H Half height of channel - Re Reynolds number = volume flow/unit width/ - Modified turbulent Reynolds number - R _t turbulent Reynolds numberk ²/ - P _k Shear production rate ofk, _t (U _i /x _j + U _j /x _i ) U _i /x _j - dP/dz Streamwise static pressure gradient - U _i Mean velocity vector (tensor notation) - U Friction velocity, _w/ where _w=–H dP/dz - W Mean velocity - W _b Bulk mean velocity through channel - y ⁺ yU /v. Unless otherwise stated, origin is at wall on trough plane of symmetry - Kinematic viscosity - _t Turbulent kinematic viscosity - Turbulence energy dissipation rate - Modified dissipation rate – 2(k ^1/2/x _j )² - Density - _k , Effective turbulent Prandtl numbers for diffusion ofk and      

Application of Rayleigh-Ritz method to forced convection turbulent heat transfer

An analytical model to predict heat transfer rates to an incompressible fluid in turbulent flow, with fully developed velocity profile, between a heated plate and a parallel, insulated plate is developed. The model employs van Driest's mixing length expression near the wall, a constant eddy diffusivitiy in the core and a constant turbulent Prandtl number. An approximate solution obtained by employing Rayleigh-Ritz method is shown to compare well with the exact solution obtained by numerical integration of the differential equations. The results are compared with the available experimental data and analytical solutions.

---

Anwendung der Rayleigh-Ritz-Methode auf die Wärmeübertragung bei erzwungener turbulenter Strömung

Zusammenfassung Es wird ein analytisches Modell zur Berechnung der Wärmeübertragung an ein inkompressibles Fluid in turbulenter Strömung mit voll ausgebildetem Geschwindigkeitsprofil zwischen einer beheizten Platte und einer dazu parallelen isolierten Platte angegeben. Das Modell verwendet van Driest's Ausdruck für die wandnahe Mischungslänge, eine konstante Wirbeldiffusivität im Kern und eine konstante turbulente PrandtlZahl. Eine Näherungslösung nach der Rayleigh-Ritz-Methode läßt sich gut mit der exakten Lösung vergleichen, die durch numerische Integration der Differentialgleichungen erhalten wurde. Die Ergebnisse werden mit verfügbaren Versuchswerten und analytischen Lösungen verglichen.

---

Nomenclature A⁺ dimensionless constant in van Driest formula - a⁺ dimensionless distance from the wall after which the eddy diffusivity of momentum is constant - b half-gap of passage - b⁺ dimensionless half-gap=bu^*/ - C_f skin friction coefficient - C_p constant pressure specific heat - d hydraulic mean diameter defined as 4xarea/perimeter=4b - h convective heat transfer coefficient - K⁺ dimensionless constant in van Driest formula - k fluid thermal conductivity - m mass flow rate of fluid - Nu Nusselt number hd/k - P pressure - Pr Prandtl number=/ - Pr_t turbulent Prandtl number=_m/ - q_w heat flux at wall - Re Reynolds number=v_md/ - T Temperature - u⁺ dimensionless velocity=V_x/u^* - u^* friction velocity= - V_x axial velocity - x axial distance from the entrance - x⁺ dimensionless distance=x/d - y distance from the heated wall - y⁺ dimensionless distance=yu^*/ Greek Symbols thermal molecular diffusivity - function equal to (_H+)/ - boundary layer thickness - _H eddy diffusivity of heat - _m eddy diffusivity of momentum - _m0 uniform eddy diffusivity of momentum in the core - dimensionless temperature - T-T_i/q_wd/k uniform heat flux - T-T_w/T_i-T_w uniform temperature - fluid kinematic viscosity - fluid density - fluid shearing stress - bulk mean temperature—fully developed region - fully developed transverse temperature profile Suffixes 1 fully developed - 2 in the entrance region - i at the inlet - m bulk mean value - w at the heated wall     

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     

Transient laminar free convection in horizontal cylinders

A numerical study of laminar natural convection inside uniformly heated, partially or fully filled horizontal cylinders is made. A coordinate transformation which simplifies the discretization of the equations of motion and energy is utilized. The resulting system of partial differential equations with their boundary conditions is solved using central differences for various Prandtl and Grashof numbers for two different grid sizes. The flow in completely filled cylinders for which experimental data are available is predicted. Close agreement between steady-state predictions and experiments is obtained for temperature and velocity profiles as well as for the streamline contours and isotherms. The technique is further demonstrated by solving the transient natural convection flow inside a partially filled horizontal cylinder with an adiabatic free surface and subjected to uniform wall heating.

---

Laminare freie Konvektion in horizontalen Zylindern

Zusammenfassung Es wurde eine numerische Berechnung der laminaren, freien Konvektion in gleichmäßig beheizten, teilweise oder ganz gefüllten, horizontalen Zylindern durchgeführt. Dabei wird eine Koordinatentransformation benützt, welche die Diskretisierung der Bewegungs- und der Energiegleichung vereinfacht. Das so resultierende System von partiellen Differentialgleichungen wird, zusammen mit seinen Randbedingungen, unter Verwendung einer Differenzenmethode für verschiedene Prandtl und Grashof-Zahlen sowie für zwei verschiedene Gittergrößen gelöst. Für den vollständig gefüllten Zylinder, für den experimentelle Daten verfügbar sind, wird die Strömung vorhergesagt. Dabei wird für stationäre Zustände gute Übereinstimmung zwischen Rechnung und Experiment erzielt. Dies gilt sowohl für den Verlauf der Stromlinien als auch für den der Isothermen. Das Verfahren wird weiterhin am Beispiel der Berechnung instationärer, freier Konvektion in einem partiell gefüllten, horizontalen Zylinder demonstriert, wobei eine adiabate, freie Oberfläche und gleichmäßige Beheizung der Wand angenommen sind.

---

Nomenclature g acceleration due to gravity, m/s² - Gr _R ^* modified Grashof number =gqR⁴/kv² - Gr _R Grashof number =gTR³/v² - H heat function vector, dimensionless - k thermal conductivity, W/mK - L(Y) cord length associated with coordinateY, dimensionless - Pr Prandtl number=v/ - q wall heat flux, W/m² - R radius, m - r(X, Y,Z) distance of a boundary point from the reference axis, dimensionless - S vector derived from the flow field solution, dimensionless - T temperature, K - T _w wall temperature, K - T reference temperature, K - t time, s - u, v velocity components inx, y directions, m/s - U, V dimensionless velocity components inX- and Y-direction normalized withU - U reference velocity=gqR²/k or gTR, m/s - V velocity vector, dimensionless - W vorticity vector, dimensionless - W vorticity, dimensionless - x, y, z cartesian coordinates, m - X, Y, Z cartesian coordinates normalized with a reference length, dimensionless Greek letters thermal diffusivity, m²/s - coefficient of thermal expansion, K^–1 - ,,, non-dimensional coordinates in the transformed domain - non-dimensional temperature =(T–T)_k/q^R or T–T/T_w–T - v kinematic viscosity, m²/s - non-dimensional time=v/R² Gr_Rt or v/R² G _R ^* t - angle measured from the bottom of the cylinder, rads - ^* angle measured from the axis on (– ) plane, rads - heat potential, dimensionless - angle of incidence of the heat flux vector, rads - non-dimensional stream function - vector potential, dimensionless - grid size, dimensionless - ² Laplacian operator - gradient vector     

Viscoelastic swirling flow with free surface in cylindrical chambers

The flow of a viscoelastic liquid driven by the steadily rotating bottom cover of a cylindrical cup is investigated. The flow field and the shape of the free surface are determined at the lowest significant orders of the regular domain perturbation in terms of the angular velocity of the bottom cap. The meridional field superposed on a primary azimuthal field shows a structure of multiple cells. The velocity field and the shape of the free surface are strongly effected by the cylinder aspect ratio and the elasticity of the liquid. The use of this flow configuration as a free surface rheometer to determine the first two Rivlin-Ericksen constants is shown to be promising.Nomenclature R, ,Z Coordinates in the physical domain D - , , Coordinates in the rest stateD ₀ - r, ,z Dimensionless coordinates in the rest stateD ₀ - Angular velocity - Zero shear viscosity - Surface tension coefficient - Density - Dimensionless surface tension parameter - ₁, ₂ The first two Rivlin-Ericksen constants - Stream function - Dimensionless second order meridional stream function - ^* Dimensionless second normal stress function - 2 Dimensionless sum of the first and second normal stress functions - N ₁,N ₂ The first and second normal stress functions - n Unit normal vector - D Stretching tensor - A _n nth order Rivlin-Ericksen tensor - S Extra-stress - u Velocity field - U Dimensionless second order meridional velocity field - V Dimensionless first order azimuthal velocity field - p Pressure - Modified pressure field - P Dimensionless second order pressure field - J Mean curvature - a Cylinder radius - d Liquid depth at rest - D Dimensionless liquid depth at rest - h Free surface height - H Dimensionless free surface height at the second order     

Flow in porous media II: The governing equations for immiscible,two-phase flow

The Stokes flow of two immiscible fluids through a rigid porous medium is analyzed using the method of volume averaging. The volume-averaged momentum equations, in terms of averaged quantities and spatial deviations, are identical in form to that obtained for single phase flow; however, the solution of the closure problem gives rise to additional terms not found in the traditional treatment of two-phase flow. Qualitative arguments suggest that the nontraditional terms may be important when / is of order one, and order of magnitude analysis indicates that they may be significant in terms of the motion of a fluid at very low volume fractions. The theory contains features that could give rise to hysteresis effects, but in the present form it is restricted to static contact line phenomena.Roman Letters (, = , , and ) A interfacial area of the- interface contained within the macroscopic system, 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 the averaging volume, 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² - g gravity vector, m²/s - H mean curvature of the- interface, m^–1 - H area average of the mean curvature, m^–1 - H – H , deviation of the mean curvature, m^–1 - I unit tensor - K Darcy's law permeability tensor, m² - K permeability tensor for the-phase, m² - K viscous drag tensor for the-phase equation of motion - K viscous drag tensor for the-phase equation of motion - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - n unit normal vector pointing from the-phase toward the-phase (n = –n ) - p _c p – P , capillary pressure, N/m² - p pressure in the-phase, N/m² - p intrinsic phase average pressure for the-phase, N/m² - p – p , spatial deviation of the pressure in the-phase, N/m² - r ₀ radius of the averaging volume, m - t time, s - v velocity vector for the-phase, m/s - v phase average velocity vector for the-phase, m/s - v intrinsic phase average velocity vector for the-phase, m/s - v – v , spatial deviation of the velocity vector for the-phase, m/s - V averaging volume, m³ - V volume of the-phase contained within the averaging volume, m³ Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m² - surface tension of the- interface, N/m - viscous stress tensor for the-phase, N/m² - / kinematic viscosity, m²/s     

Unsteady convective diffusion in a rotating parallel plate channel

The paper presents an exact analysis of the dispersion of a passive contaminant in a viscous fluid flowing in a parallel plate channel driven by a uniform pressure gradient. The channel rotates about an axis perpendicular to its walls with a uniform angular velocity resulting in a secondary flow. Using a generalized dispersion model which is valid for all time, we evaluate the longitudinal dispersion coefficientsK _i (i=1, 2, ...) as functions of time. It is shown thatK ₁=0 andK ₃,K ₄, ... decay rapidly in comparison withK ₂. ButK ₂ decreases with increasing (the dimensionless rotation parameter) for values of upto approximately =2.2. ThereafterK ₂ increases with further increase in and its value gets saturated for large values of (say, 500) and does not change any further with increase in . A physical explanation of this anomalous behaviour ofK ₂ is given.

---

Instationäre konvektive Diffusion in einem rotierenden Parallelplattenkanal

Zusammenfassung In dieser Untersuchung wird eine exakte Analyse der Ausbreitung eines passiven Kontaminierungsstoffes in einer zähen Flüssigkeit gegeben, die, befördert durch einen gleichförmigen Druckgradienten, in einem Parallelplattenkanal strömt. Der Kanal rotiert mit gleichförmiger Winkelgeschwindigkeit um eine zu seinen Wänden senkrechte Achse, wodurch sich eine Sekundärströmung ausbildet. Unter Verwendung eines generalisierten, für alle Zeiten gültigen Dispersionsmodells werden die longitudinalen DispersionskoeffizientenK _i (i=1, 2, ...) als Funktionen der Zeit ermittelt. Es wird gezeigt, daßK ₁=0 gilt und dieK ₃,K ₄, ... gegenüberK ₂ schnell abnehmen.K ₂ nimmt ab, wenn , der dimensionslose Rotationsparameter, bis etwa zum Wert 2,2 ansteigt. Danach wächstK ₂ mit bis auf einem Endwert an, der etwa ab =500 erreicht wird. Dieses anomale Verhalten vonK ₂ findet eine physikalische Erklärung.

---

List of symbols C solute concentration - D molecular diffusivity - K _i longitudinal dispersion coefficients - 2L depth of the channel - P ₀ dimensionless pressure gradient along main flow - Pe Péclet number - q velocity vector - Q _x,Q _y mass flux along the main flow and the secondary flow directions - dimensionless average velocity along the main flow direction - (x, y, z) Cartesian co-ordinates Greek symbols dimensionless rotation parameter - the inclination of side walls withx-axis - kinematic viscosity - fluid density - dimensionless time - angular velocity of the channel - dimensionless distance along the main flow direction - dimensionless distance along the vertical direction - dimensionless solute concentration - integral of the dispersion coefficientK ₂() over a time interval     

Injection moulding of thermoplastics: Freezing of variable-viscosity fluids.

The injection moulding of thermoplastics involves, during mould filling, flows of hot polymer melts into mould networks, the walls of which are so cold that frozen layers form on them. An analytical study of such flows is presented here for the case when the Graetz and Nahme numbers are large and the Pearson number is small. Thus the flows are developing and temperature differences due to heat generation by viscous dissipation are sufficiently large to cause significant variations in viscosity (but the difference between the entry temperature of the polymer to a specific part of the mould network and the melting temperature of the polymer is not). Br Brinkman number - Gz Graetz number - h half-height of channel or disc - h ^* half-height of polymer melt region in channel or disc - L length of channel or pipe - m viscosity shear-rate exponent - Na Nahme number - p pressure - P pressure drop - Pe Péclet number - Pn Pearson number - Q volumetric flowrate - r radial coordinate in pipe or disc - R radius of pipe - Re Reynolds number - R _i inner radius of disc - R _o outer radius of disc - R ^* radius of polymer melt region in pipe - T temperature - T _ad adiabatic temperature rise - T _e entry polymer melt temperature - T _m melting temperature of polymer - T _max maximum temperature - T ₀ reference temperature - T _w wall temperature - flow-average temperature rise - u _r radial velocity in pipe or disc - u _x axial velocity in channel - u _y transverse velocity in channel or disc - u _z axial velocity in pipe - w width of channel - x axial coordinate in channel or modified radial coordinate in disc - y transverse coordinate in channel or disc - z axial coordinate in pipe - thermal conductivity of molten polymer - thermal conductivity of frozen polymer - scaled dimensionless axial coordinate in channel or pipe or radial coordinate in disc - ₀ undetermined integration constant - heat capacity of molten polymer - viscosity temperature exponent - dimensionless transverse coordinate in channel or disc - ^* dimensionless half-height of polymer melt region in channel or disc - H ^* scaled dimensionless half-height of polymer melt region in channel or disc or radius of polymer melt region in pipe - dimensionless temperature - ^* dimensionless wall temperature - scaled dimensionless temperature - numerical constant - µ viscosity of molten polymer - µ ₀ consistency of molten polymer - dimensionless pressure gradient - scaled dimensionless pressure gradient - density of molten polymer - dimensionless radial coordinate in pipe or disc - _i dimensionless inner radius of disc - ^* dimensionless radius of polymer melt region in pipe - dimensionless streamfunction - scaled dimensionless streamfunction - dummy variable - streamfunction - similarity variable - similarity variable     

The Gibbs-Thompson relation within the gradient theory of phase transitions

This paper discusses the asymptotic behavior as 0+ of the chemical potentials associated with solutions of variational problems within the Van der Waals-Cahn-Hilliard theory of phase transitions in a fluid with free energy, per unit volume, given by ²¦¦²+ W(), where is the density. The main result is that is asymptotically equal to E/d+o(), with E the interfacial energy, per unit surface area, of the interface between phases, the (constant) sum of principal curvatures of the interface, and d the density jump across the interface. This result is in agreement with a formula conjectured by M. Gurtin and corresponds to the Gibbs-Thompson relation for surface tension, proved by G. Caginalp within the context of the phase field model of free boundaries arising from phase transitions.     

Vorticity characteristics of the turbulent intermediate wake

Measurements of the lateral components _j (j=2 and 3) of the vorticity fluctuation vector have been made, using a vorticity probe consisting of two X-wires, in the intermediate wake of a circular cylinder. The effect of the spatial resolution of the probe on the measurement of _j has been studied. As the spatial resolution impairs, the variance and flatness factor of _j decrease whereas the skewness of _j increases. Reasonably accurate values of _j ² can be obtained by applying spectral corrections for the spatial resolution effect.Near the beginning of the intermediate wake, the variance of ₂ is larger than that of ₃ due to the significant contribution from ribs which connect consecutive spanwise roll vortices. This difference decreases with downstream distance. Also, the presence of the rolls is reflected by a local extremum in the skewness of ₃ on each side of the wake centerline. The magnitude of the extremum decreases with downstream distance.The support of the Australian Research Council is gratefully acknowledged.     

Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

Magnetohydrodynamic natural convection heat transfer from radiate vertical porous surfaces

Magnetohydrodynamic natural convection heat transfer from radiate vertical surfaces with fluid suction or injection is considered. The nonsimilarity parameter is found to be the conductive fluid injection or suction along the streamwise coordinate = V{4x/² g(T _w – T )}^1/4. Three dimensionless parameters had been found to describe the problem: the magnetic influence number N = B ² _y /V ², the radiation-conduction parameter R _d = k _R /4aT ³ , and the Gebhart number Ge _x = gx/cp to represent the effect of the viscous dissipation. It is found that increasing the magnetic field strength causes the velocity and the heat transfer rates inside the boundary layer to decrease. Its apparent that increasing the radiation-conduction parameter decreases the velocity and enhances the heat transfer rates. The Gebhart number, i.e, the viscous dissipation had no effect on the present problem.Nomenclature a Stefan-Boltzmann constant - B _y Magnetic field flux density Wb/m² - Cf _x Local skin friction factor - c _p Specific heat capacity - f Dimensionless stream function - Ge _x Gebhart number, gx/cp - g Gravitational acceleration - k Thermal Conductivity - L Length of the plate - N Magnetic influence number, B ² _y /V ² - p Pressure - Pr Prandtl number - q _r Radiative heat flux - q _w (x) Local surface heat flux - Q _w (x) Dimensionless Local surface heat flux - R _d Planck number (Radiation-Conduction parameter), k _R /4aT ³ - T Temperature - T Free stream temperature - T _w Wall temperature - u, v Velocity components in x- and y-directions - V Porous wall suction or injection velocity - V _w Porous wall suction or injection velocity - x, y Axial and normal coordinates - Thermal diffusivity Greek symbols _R Roseland mean absorption coefficient, 4/3R _d - Coefficient of thermal expansion - Nonsimilarity parameter, V{4x/² g(T _w – T )}^1/4 - Peseudo-similarity variable - Dimensionless temperature - _w Ratio of surface temperature to the ambient temperature, T _w /T - Dynamice viscosity - Kinemtic viscosity - Fluid density - Electrical conductivity - _w Local wall shear stress - Dimensional stream function     

Nonstationary vibration of a rotating shaft with nonlinear spring characteristics during acceleration through a critical speed (A critical speed of a summed-and-differential harmonic oscillation)

Nonstationary vibration of a flexible rotating shaft with nonlinear spring characteristics during acceleration through a critical speed of a summed-and-differential harmonic oscillation was investigated. In numerical simulations, we investigated the influence of the angular acceleration , the initial angular position of the unbalance _n and the initial rotating speed on the maximum amplitude. We also performed experiments with various angular accelerations. The following results were obtained: (1) the maximum amplitude depends not only on but also on _n and : (2) when the initial angular position _n changes. the maximum amplitude varies between two values. The upper and lower bounds of the maximum amplitude do not change monotonously for the angular acceleration: (3) In order to always pass the critical speed with finite amplitude during acceleration. the value of must exceed a certain critical value.Nomenclature O-xyz rectangular coordinate system - , ₁, ₁ inclination angle of rotor and its projections to thexy- andyz-planes - I _r polar moment of inertia of rotor - I diametral moment of inertia of rotor - i _r ratio ofI _r toI - dynamic unbalance of rotor - directional angle of fromx-axis - c damping coefficient - spring constant of shaft - N _nt ,N _nt nonlinear terms in restoring forees in ₁ and ₁ directions - ₄ representative angle - a small quantity - V. V _u .V _N potential energy and its components corresponding to linear and nonlinear terms in the restoring forees - directional angle - _n coefficients of asymmetrical nonlinear terms - _n coefficients of symmetrical nonlinear terms - coefficients of asymmetrical nonlinear terms experessed in polar coordinates - coefficients of symmetrical nonlinear terms expressed in polar coordinates - rotating speed of shaft - t time - _n initial angular position of att=0 - p natural frequency - p ₁.p _t natural frequencies of forward and backward precessions - , ₁, ₁ total phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - , ₁, ₁ phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - P, R _t ,R _b amplitudes of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - difference between phases ( = _f – _u) - acceleration of rotor - initial rotating speed - t _t ,r _b amplitudes of nonstationary oscillation during acceleration - (r _t )_max, (r _b )_max maximum amplitudes of nonstationary oscillation during acceleration - (r ₁ ¹ )_max, (r _b ¹ )_max maximum value of angular acceleration of non-passable case - ₀ critical value over which the rotor can always pass the critical speed - p ₁,p ₂,p ₃,p ₄ natural frequencies of experimental apparatus     

Translatory accelerating motion of a circular disk in a viscous fluid

The exact solution of the equation of motion of a circular disk accelerated along its axis of symmetry due to an arbitrarily applied force in an otherwise still, incompressible, viscous fluid of infinite extent is obtained. The fluid resistance considered in this paper is the Stokes-flow drag which consists of the added mass effect, steady state drag, and the effect of the history of the motion. The solutions for the velocity and displacement of the circular disk are presented in explicit forms for the cases of constant and impulsive forcing functions. The importance of the effect of the history of the motion is discussed.Nomenclature a radius of the circular disk - b one half of the thickness of the circular disk - C dimensionless form of C ₁ - C ₁ magnitude of the constant force - D fluid drag force - f(t) externally applied force - F() dimensionaless form of applied force - F ₀ initial value of F - g gravitational acceleration - H() Heaviside step function - k magnitude of impulsive force - K dimensionless form of k - M a dimensionless parameter equals to (1+37#x03C0;_s/4_f) - S displacement of disk - t time - t ₁ time of application of impulsive force - u velocity of the disk - V dimensionless velocity - V ₀ initial velocity of V - V _t terminal velocity - parameter in (13) - parameter in (13) - (t) Dirac delta function - ratio of b/a - () function given in (5) - dynamical viscosity of the fluid - kinematic viscosity of the fluid - _f fluid density - _s mass density of the circular disk - dimensionless time - _i dimensionless form of t _i - dummy variable - dummy variable