On the effects of a periodic, spanwise variation of suction upon asymptotic suction flow

Summary The effects of superposing streamwise vorticity, periodic in the lateral direction, upon two-dimensional asymptotic suction flow are analyzed. Such vorticity, generated by prescribing a spanwise variation in the suction velocity, is known to play an important role in unstable and turbulent boundary layers. The flow induced by the variation has been obtained for a freestream velocity which (i) is steady, (ii) oscillates periodically in time, (iii) changes impulsively from rest. For the oscillatory case it is shown that a frequency can exist which maximizes the induced, unsteady wall shear stress for a given spanwise period. For steady flow the heat transfer to, or from a wall at constant temperature has also been computed.Nomenclature (x, y, z) spatial coordinates - (u, v, w) corresponding components of velocity - (, , ) corresponding components of vorticity - t time - stream function for v and w - v _w mean wall suction velocity - nondimensional amplitude of variation in wall suction velocity - characteristic wavenumber for variation in direction of z - T temperature - P pressure - density - coefficient of kinematic viscosity - coefficient of thermal diffusivity - (/v _w)² - frequency of oscillation of freestream velocity - nondimensional amplitude of freestream oscillation - /v _w ² - z z - y –v _w y/ - v _w ² t/4 - /v _w - U ₀ characteristic freestream velocity - u/U ₀ - coefficient of viscosity - _w wall shear stress - Prandtl number (/) - q heat transfer to wall - T _w wall temperature - T (T _w–T)/(T _w–)     

Nicht-isotherme Rieselfilmströmung einer hochviskosen Flüssigkeit mit stark temperaturabhÄngiger ViskositÄt

Zusammenfassung Es werden Geschwindigkeitsverteilungen und Filmdickenabnahmen von nichtisothermen NEWTONschen und nicht-NEWTONschen (Potenzansatz) Rieselfilmen mit temperaturanhÄngiger ViskositÄt berechnet, wobei die Temperaturverteilung im Film als linear vorausgesetzt wird. An dicken Rieselfilmen mit Re=10^–4... 10^–2 sind Geschwindigkeitsprofile, Filmdicken und OberflÄchentemperaturen gemessen und daraus die thermische EinlauflÄnge bestimmt worden. Ausgehend von der Penetrationstheorie für eine endlich dicke Platte kann man für diese EinlauflÄnge eine Approximationsformel erhalten, die für Strömungen mit Re < 1000 verwendet werden kann.

---

Non-isothermal filmflow of a highly viscous liquid, the viscosity strongly depending on temperature

Velocity distributions and film thicknesses of nonisothermal NEWTONIAN and non-NEWTONIAN (power-law) falling films are computed assuming that the temperature across the film varies linearly. Experimental studies on thick falling films of Re=10^–4...10^–2 had been carried out to measure velocities, film thickness and surface temperature and to calculate the thermal entrance length. One can get for this entrance length a approximation formula which is valid for flows with RePr <1000 by applying the results for the thermal penetration into a material finite plate.

---

Bezeichnungen B dimensionsloser Temperaturkoeffizient - ¯B [K] Temperaturkoeffizient (ln)/(1/T) - c_p [J/kgK] spezif. WÄrme bei konst. Druck - Fo FOURIER-Zahl - g [m/s²] Erdbeschleunigung - H dimensionslose Filmdicke - h [m] Filmdicke - m [Pas^2–n] ViskositÄtskoeffizient im Potenzansatz von OSTWALD-DE WAELE - Nu NUSSELT-Zahl - n Flüssigkeitsexponent im Potenzansatz von OSTWALD-DE WAELE - Pr PRANDTL-Zahl (Gl.3.5) - q [W/m²] WÄrmestromdichte - Re REYNOLDS-Zahl (Gl.3.4) - T [K] Temperatur - t [s] Zeit - U dimensionslose Geschwindigkeit (X-Komponente) - u [m/s] Geschwindigkeitskomponente in x-Richtung - X dimensionslose Koordinate (X=x/h₀) - x [m] LÄnge, Koordinate - Y dimensionslose Koordinate (Y=y/h₀) - y [m] Höhe, Koordinate - [W/m²K] WÄrmeübergangskoeffizient - Plattenneigungswinkel gegen Horizontale - [s^–1] Schergeschwindigkeit - dimensionslose Temperatur (Gl.3.3) - [m²/s] TemperaturleitfÄhigkeit (Gl.3.3) - [W/mK] WÄrmeleitfÄhigkeit - [Pas] ViskositÄt - [kg/m³] spezif. Dichte - [Pa] Schubspannung Indizes a scheinbar (apparent) - 0 bei x=0, auch: isotherm - P auf die Penetrationszeit bezogen - s an der OberflÄche - T bei linearer Temperaturdifferenz T - w an der Wand - 99 auf =0,99 bezogen - gemittelt, Mittelwert - thermisch ausgebildet, bei x - proportional - ¯t ungefÄhr - kleiner oder gleich ungefÄhr     

Analytical solution for transient laminar fully developed free convection in vertical channels

Using Green's function method, analytical solutions for transient fully developed natural convection in open-ended vertical circular and two-parallel-plate channels are presented. Different fundamental boundary conditions for these two configurations have been investigated and the corresponding fundamental solutions are obtained. These fundamental solutions may be used to obtain solutions satisfying more general thermal boundary conditions. In terms of the obtained unsteady temperature and velocity profiles, the transient volumetric flow rate, mixing cup emperature and local nusselt number are estimated.Zusammenfassung Für oben und unten offene vertikale Kanäle mit Kreisquerschnitt bzw. als Parallelplattenanordnung werden unter Verwendung der Methode der Greenschen Funktionen analytische Lösungen für die nichtstationäre, vollausgebildete, natürliche Konvektion gefunden und zwar unter Zugrundelegung verschiedener Fundamental-Randbedingungen bezüglich beider Konfigurationen. Die so ermittelten Fundamentallösungen können zur Gewinnung von Lösungen für allgemeine Randbedingungen dienen. Der zeitlich veränderliche Volumenstrom, die Mischtemperatur und die Nusselt-Zahl werden mit Bezug auf die erhaltenen nichtstationären Profile für Temperatur und Geschwindigkeit näher analysiert.

---

Analytische Lösung für die nichtstationäre vollausgebildete laminare freie Konvektion invertikalen Kanälen

Nomenclature a local heat transfer coefficient based on the area of the heat transfer surface,q/(T _w –T ₀)=±(T/y)_w/(T_w–T₀), minus and plus signs apply respectively for heating and cooling in case of parallel-plate channel and vice versa in case of a tube - average heat transfer coefficient over the channel - c _p specific heat of fluid at constant pressure - f volumetric flow rate, for circular channels and or two-parallel-plate channels - F dimensionless volumetric flow rate,f/(2lvGr^*) for circular channels forfw/(lvGr ^*) for two-parallel-plate channels - g gravitational body force per unit mass (acceleration) - G Green's function - Gr Grashof number,±g(T _w–T₀)w³/v² in case of an isothermal boundary of±gqw ⁴/2kv² in case of a uniform heat flux (UHF) on the heat transfer boundary, the plus and minus signs apply to upward (heating) and downward (cooling) flows, respectively. ThusGr is a positive number in both cases. - Gr ^* modified Grashof number,wGr/l - h heat gained or lost by fluid from the entrance up to a particular elevation in the channel, ₀ fc _p(T _m–T ₀) for all cases - J ₀ Bessel function of zero order - k thermal conductivity of fluid - l height of channel - L dimensionless height of channel,1/Gr ^* - Nu local Nusselt number,|a| w/k - average Nusselt number, - p pressure of fluid inside the channel at any cross-section - p pressure defect at any point,p–p _s - p ₀ pressure of fluid at the channel entrance - p _s hydrostatic pressure, ₀ gz where the minus and plus signs are for upward (heating) and downward (cooling) flows, respectively - p dimensionless pressure defect at any point(pw ⁴)/(₀ l ²² Gr ²) - Pr Prandtl number,c _p/k - q heat flux at the heat transfer surface,q=±k(T/y)_w where the minus and plus signs are, respectively, for cooling and heating in case of circular pipe and vice versa in case of a parallel-plate channel - Ra Rayleigh number,GrPr - Ra ^* modified Rayleigh number,Gr ^*Pr - t time - T fluid temperature at any point - T _m mixing-cup (mixed-mean) temperature over any cross section, for circular channels, and for two-parallelplate channels - T ₀ initial and channel-inlet fluid temperature - T _w temperature of the heat-transfer wall - u axial velocity component at any point - U dimensionless axial velocity,uw ²/(lvGr^*) - w radius of circular tube or width (between plates) of parallel-plate channel - y radial or transverse coordinate - y dimensionless radial or transverse coordinate,y/w - z axial coordinate - Z dimensional axial coordinate,z/(lGr ^*) Greek symbols constant appears in Eq. (8) - parameter appears in Eq. (9) which equals the integration of with respect to or volumetric coefficient of thermal expansion - _n eigenvalues - parameter appears in Eq. (7) - _n eigenvalues - parameter appears in Eq. (12) - _n eigenvalues - parameter appears in Eq. (9) - dimensionless temperature,(T–T ₀)/(T_w–T₀) in case of an isothermal heat transfer boundary and(T–T ₀)/(qw/2k) for UHF boundary - _m dimensionless mixing cup temperature,(T _m–T₀)/(T_w–T₀) in case of an isothermal heat transfer boundary and(T _m–T₀)/(qw/2k) for UHF boundary - _w dimensionless temperature of the heat-transfer wall, equals unity in case of an isothermal heat transfer boundary and(T _w–T₀)/(qw/2k) for a UHF boundary - _n eigenvalues - dynamic viscosity of fluid - kinematic viscosity of fluid, /₀ - fluid density at temperatureT,₀[1–(T–T ₀)] - ₀ fluid density atT ₀ - demensionless time,tk/(cw²)     

Diffusion of moisture in deformable porous media. Anisotropic fibrous materials

This paper presents a study on the deformation of anisotropic fibrous porous media subjected to moistening by water in the liquid phase. The deformation of the medium is studied by applying the concept of effective stress. Given the structure of the medium, the displacement of the solid matrix is not taken into account with respect to the displacement of the liquid phase. The transport equations are derived from the model proposed by Narasimhan. The transport coefficients and the relation between the variation in apparent density and effective stress are obtained by test measurements. A numerical model has been established and applied for studying drip moistening of mineral wool samples capable or incapable of deformation.Nomenclature D mass diffusion coefficient [L²t^–1] - e void fraction - g gravity acceleration [Lt^–2] - J mass transfer density [ML^–2t^–1] - K hydraulic conductivity [Lt^–1] - K _s hydraulic conductivity of the solid phase [Lt^–1] - K ^* hydraulic conductivity of the deformable porous medium [Lt^–1] - P pressure of moistening liquid [ML^–1 t^–2] - S degree of saturation - t time [t] - V speed [Lt^–1] - X horizontal coordinate [L] - Z vertical coordinate measured from the bottom of porous medium [L] - z z-coordinate [L] Greek Letters porosity - ₁ total hydric potential [L] - _g gas density [ML^–3] - ₁ liquid density [ML^–3] - ₀ apparent density [ML^–3] - _s density of the solid phase [ML^–3] - density of the moist porous medium [ML^–3] - external load [ML^–1t^–2] - effective stress [ML^–1t^–2] - bishop's parameter - matrix potential or capillary suction [L] Indices g gas - 1 moistening liquid - p direction perpendicular to fiber planes - s solid matrix - t direction parallel to fiber planes - v pore Exponent * movement of solid particles taken into account     

An analytical solution for the temperature profiles during injection molding,including dissipation effects

An analytical solution is obtained for the stationary temperature profile in a polymeric melt flowing into a cold cavity, which also takes into account viscous heating effects. The solution is valid for the injection stage of the molding process. Although the analytical solution is only possible after making several (at first sight) rather stringent assumptions, the calculated temperature field turns out to give a fair agreement with a numerical, more realistic approach. Approximate functions were derived for both the dissipation-independent and the dissipation-dependent parts which greatly facilitate the temperature calculations. In particular, a closed-form expression is derived for the position where the maximum temperature occurs and for the thickness of the solidified layer.The expression for the temperature field is a special case of the solution of the diffusion equation with variable coefficients and a source term.Nomenclature a thermal diffusivity [m²/s] - c specific heat [J/kg K] - D channel half-height [m] - L channel length [m] - m 1/ - P pressure [Pa] - T temperature [°C] - T _W wall temperature [°C] - T _i injection temperature [°C] - T _A Br independent part of T - T _B Br dependent part of T - T ^core asymptotic temperature - v _z() axial velocity [m/s] - W channel width [m] - x cross-channel direction [m] - z axial coordinate [m] - (x) gamma function - (a, x) incomplete gamma function - M(a, b, x) Kummer function - small parameter - () temperature function - thermal conductivity [W/mK] - viscosity [Pa · s] - ₀ consistency index - power-law exponent - density [kg/m] - similarity variable Dimensionless variables Br Brinkman number - Gz Graetz number -      

Maximum density effects on thermal instability of horizontal laminar boundary layers

Linear stability theory is used to investigate the onset of longitudinal vortices in laminar boundary layers along horizontal semi-infinite flat plates heated or cooled isothermally from below by considering the density inversion effect for water using a cubic temperature-density relationship. The analysis employs non-parallel flow model incorporating the variation of the basic flow and temperature fields with the streamwise coordinate as well as the transverse velocity component in the disturbance equations. Numerical results for the critical Grashof number Gr _L ^* =Gr _X ^* /Re _X< ^{Emphasis>/3/2} are presented for thermal conditions corresponding to –0.5 ₁–2.0 and –0.8 ₂1.2.Nomenclature a wavenumber, 2/ - D operator, d/d - F (f–f)/2 - f dimensionless stream function - g gravitational acceleration - G eigenvalue, Gr _L/Re_L - Gr _L Grashof number based on L - Gr _X Grashof number based on X - L characteristic length, (X/U)^1/2 - M number of divisions in y direction - P pressure - Pr Prandtl number, / - p dimensionless pressure, P/( ² /Re _L) - Re _L, Re_X Reynolds numbers, (U L/)=Re _X< ^1/2 and (U), respectively - T temperature - U, V, W velocity components in X, Y, Z directions - u, v, w dimensionless perturbation velocities, (U, V, W)/U - X, Y, Z rectangular coordinates - x, y, z dimensionless coordinates, (X, Y, Z)/L - thermal diffusivity - coefficient of thermal expansion - ₁, ₂ temperature coefficients for density-temperature relationship - similarity variable, Y/L=y - dimensionless temperature disturbance, /T - dimensionless wavelength of vortex rolls, 2/a - ₁, ₂ thermal parameters defined by equation (12) - kinematic viscosity - density - dimensionless basic temperature, (T _b T )/T - –1 - T temperature difference, (T _w–T ) - * critical value or dimensionless disturbance amplitude - prime, disturbance quantity or differentiation with respect to - b basic flow quantity - max value at a density maximum - w value at wall - free stream condition     

Thermal stability of composite superconducting tape under the effect of a two-dimensional dual-phase-lag heat conduction model

Thermal stability of composite superconducting tape subjected to a thermal disturbance is numerically investigated under the effect of a two-dimensional dual-phase-lag heat conduction model. It is found that the dual-phase-lag model predicts a wider stable region as compared to the predictions of the parabolic and the hyperbolic heat conduction models. The effects of different design, geometrical and operating conditions on superconducting tape thermal stability were also studied.a conductor width, (m) - A conductor cross sectional area of, (m²) - As conductor aspect ratio, (a/b) - b conductor thickness, (m) - Bi Biot number - B dimensionless disturbance Intensity - C heat capacity, (J m^–3 K^–1) - D disturbance energy density, (W m^–3) - f volume fraction of the stabilizer in the conductor - g(T) steady capacity of the Ohmic heat source, (W m^–3) - g_max Ohmic heat generation with the whole current in the stabilizer, (W m^–3) - G_max dimensionless maximum Joule heating - h convective heat transfer coefficient, (W m^–2 K^–1) - J current density, (A m^–2) - k thermal conductivity of conductor, (W m^–1 K^–1) - q conduction heat flux vector, (W m^–2) - Q dimensionless Joule heating - R relaxation times ratio (_T/2_q) - t rime, (s) - T temperature, (K) - T_c critical temperature, (K) - T_c1 current sharing temperature, (K) - T_i initial temperature, (K) - T_o ambient temperature, (K) - x, y co-ordinate defined in Fig. 1, (m) - thermal diffusivity (m² s^–1) - dimensionless time - _i dimensionless duration time - dimensionless y-variable - _o superconductor dimensionless thickness - dimensionless temperature - _c1 dimensionless current sharing temperature - 1 dimensionless maximum temperature - dimensionless disturbance energy - numerical tolerance - x width of conductor subjected to heat disturbances, (m) - y thickness of conductor subjected to heat disturbances, (m) - dimensionless x-variable - _o superconductor dimensionless width - stabilizer electrical resistivity, () - _q relaxation time of heat flux, (s) - _T relaxation time of temperature gradient, (s) - i initial - sc current sharing - max maximum - o ambient     

Linearized analysis of magnetohydrodynamic entrance flow in a vertical channel with combined thermal convection

A linearized analysis is presented for the magnetohydrodynamic entrance flow with combined forced and free convection in a vertical, constant wall temperature parallel-plate channel. Numerical results are obtained for slug velocity profile at the entrance and for various Hartmann and Grashof Numbers. The results agree well with the finite difference numerical solutions obtained elsewhere. They demonstrate that the velocity development and pressure gradient in the channel entrance region are greatly influenced by the Hartmann Number and the Grashof Number. Increasing Hartmann Number decreases velocity entrance length while increasing Grashof Number increases it. Thermal development is also found to be dependent on the above mentioned parameters, but to a relatively minor extent.Nomenclature A _m constant defined by equation (23) - B ₀ applied magnetic field - C _n constant defined by equation (13) - E ₀ constant electric field - e nondimensional electric field parameter, E ₀/U _mB₀ - Gr Grashof Number, gL ³(T _w–T ₀)/ ² - L half-width of the channel - M Hartmann Number, B ₀ L(/)^1/2 - Nu Nusselt Number, (/y) _y=1/( _w– _m) - P pressure - Pr Prandtl Number, / - p nondimensional pressure parameter, (P–P ₀+ ₀ gX)/P ₀ U _m ² - Re Reynolds Number, U _m L/ - T temperature - T ₀ inlet temperature - T _w wall temperature - U velocity, X direction - U _m average velocity, (1/L) ₀ ^L UdY - u nondimensional form of U, U/U _m - u ₀(y) nondimensional inlet velocity - V velocity, Y direction - v nondimensional form of V, VL/ - X coordinate, axial direction - x nondimensional form of X, vX/L ² U _m - Y coordinate perpendicular to the channel - y nondimensional form of Y, Y/L - thermal diffusivity - _m eigenvalue defined by equation (25) - thermal expansion coefficient - _m eigenvalue defined by equation (24) - stretching factor, weighting function - nondimensional form of T, (T–T ₀)/(T _w–T ₀) - _m mean nondimensional temperature, ₀ ¹ udy - kinematic viscosity - magnetic permeability - mass density - electrical conductivity     

Mixed convection on vertical plate for fluids of any prandtl number

A mixed convection parameter=(Ra) ^1/4/(Re)^1/2, with=Pr/(1+Pr) and=Pr/(1 +Pr)^1/2, is proposed to replace the conventional Richardson number, Gr/Re², for combined forced and free convection flow on an isothermal vertical plate. This parameter can readily be reduced to the controlling parameters for the relative importance of the forced and the free convection,Ra ^1/4/(Re ^1/2 Pr ^1/3) forPr 1, and (RaPr)^1/2/(RePr ^1/2 forPr 1. Furthermore, new coordinates and dependent variables are properly defined in terms of, so that the transformed nonsimilar boundary-layer equations give numerical solutions that are uniformly valid over the entire range of mixed convection intensity from forced convection limit to free convection limit for fluids of any Prandtl number from 0.001 to 10,000. The effects of mixed convection intensity and the Prandtl number on the velocity profiles, the temperature profiles, the wall friction, and the heat transfer rate are illustrated for both cases of buoyancy assisting and opposing flow conditions.

---

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

Zusammenfassung Für die kombinierte Zwangs- und freie Konvektion an einer isothermen senkrechten Platte wird ein Mischkonvektions-Parameter=( Ra) ^1/4 (Re)^1/2, mit=Pr/(1 +Pr) und=Pr/(1 +Pr)^1/2 vorgeschlagen, den die gebräuchliche Richardson-Zahl, Gr/Re², ersetzen soll. Dieser Parameter kann ohne weiteres auf die maßgebenden Kennzahlen für den relativen Einfluß der erzwungenen und der freien Konvektion reduziert werden,Ra ^1/4/(Re ^1/2 Pr ^1/3) fürPr 1 und (RaPr)^1/4/(RePr)^1/2 fürPr 1. Weiterhin werden neue Koordinaten und abhängige Variablen als Funktion von definiert, so daß für die transformierten Grenzschichtgleichungen numerische Lösungen erstellt werden können, die über den gesamten Bereich der Mischkonvektion, von der freien Konvektion bis zur Zwangskonvektion, für Fluide jeglicher Prandtl-Zahl von 0.001 bis 10.000 gleichmäßig gültig sind. Der Einfluß der Intensität der Mischkonvektion und der Prandtl-Zahl auf die Geschwindigkeitsprofile, die Temperaturprofile, die Wandreibung und den Wärmeübergangskoeffizienten werden für die beiden Fälle der Strömung in und entgegengesetzt zur Schwerkraftrichtung dargestellt.

---

Nomenclature C _f local friction coefficient - C _p specific heat capacity - f reduced stream function - g gravitational acceleration - Gr local Grashoff number,g T _w –T )x³/v² - Nu local Nusselt number - Pr Prandtl number,v/ - Ra local Rayleigh number,g T _w –T x ³/( v) - Re local Reynolds number,u x/v - Ri Richardson number,Gr/Re ² - T fluid temperature - T _w wall temperature - T free stream temperature - u velocity component in thex direction - u free stream velocity - v velocity component in they direction - x vertical coordinate measuring from the leading edge - y horizontal coordinate Greek symbols thermal diffusivity - thermal expansion coefficient - mixed convection parameter (Ra)1/4/Re)^1/2 - pseudo-similarity variable,(y/x) - ₀ conventional similarity variable,(y/x)Re ^1/2 - dimensionless temperature, (T–T T _W –T - unified mixed-flow parameter, [(Re) ^1/2 + (Ra)^1/4] - dynamic viscosity - kinematic viscosity - stretched streamwise coordinate or mixed convection parameter, [1 + (Re)^1/2/(Ra) ^1/4]^–1=/(1 +) - density - Pr/(1 + Pr) _w wall shear stress - stream function - Pr/(l+Pr)^1/3 This research was supported by a grand from the National Science Council of ROC     

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

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 boundary conditions at the macroscopic level

We study the problem of the boundary conditions specified at the boundary of a porous domain in order to solve the macroscopic transfer equations obtained by means of the volume-averaging method. The analysis is limited to the case of conductive transport but the method can be extended to other cases. A numerical study enables us to illustrate the theoretical results in the case of a model porous medium. Roman Letters _sf interfacial area of the s-f interface contained within the macroscopic system m² - A _sf interfacial area of the s-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - d _s , d _f microscopic characteristic length m - g vector that maps to _s, m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l REV characteristic length, m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for the f-phase at the f-s interface - n _e outwardly directed unit normal vector at the dividing surface - T ^* macroscopic temperature field obtained by solving the macroscopic equation (3), K - V averaging volume, m³ - V _s , V _f volume of the considered phase within the averaging volume, m³ - volume of the macroscopic system, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² Greek Letters _s , _f volume fraction - ratio of thermal conductivities - _s , _f thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding to T ^* , K - spatial deviation temperature K - error on the temperature due to the macroscopic boundary conditions, K - spatial average - ^s , ^f intrinsic phase average     

General viscosity-composition relationship for dispersions,solutions and binary liquid systems

An equation for the viscosity of a mixture of two imaginary Newtonian liquids is derived. In the derivation the mathematical assumption is used that the effective activation energy for viscous flow of a binary liquid mixture is a linear combination of the reciprocals of the activation energy of the components. It contains two dependent fitting constants and has the same structure as the Mooney equation for dispersions of spherical solid particles, the Huggins equation for polymer solutions and is identical to an equation by Hoffmann and Rother, when written in the variables that the last authors used.As a consequence it can be shown that the viscosity of binary liquid mixtures, liquid resion solutions, dispersions of solid spherical particles and polymer solutions can be described very well by one and the same equation, up to the highest concentrations.It has further been found that the viscosity of dispersions of non-spherical particles, solutions of solids in organic solvents and solutions of electrolytes and non-electrolytes in water can also be described by this formula. The equation permits the construction of a straight line on which all liquids can be plotted.An algebraic analysis of the equation shows that each series of viscosity composition data can be placed in one of three rheological groups independent of the type of fraction that is used to characterize the composition.Seventy-four binary systems, covering a wide range of liquids have been used to show the applicability of the developed equation.It has been found that in most cases the data are best described by splitting them into two regions, each with its own set of dependent constants. General symbol for the fraction or concentration of the component with the higher viscosity determining the composition of a binary mixture [—] - _v Volume fraction of the component with the higher viscosity [—] - _w Weight fraction of the component with the higher viscosity [—] - _mw Molecular weight fraction of the component with the higher viscosity [—] - c Concentration of the component with the higher viscosity [g/cm³] - E ₂,E ₁,E Activation energy for viscous flow referring to the component with the higher viscosity, the lower viscosity and the viscosity of the binary mixtures, respectively [J] - ₂, ₁, Experimental parameter (with the dimension of energy) referring to the component with the higher viscosity, the component with the lower viscosity and to the binary mixtures, respectively [J] - ₁, ₂ Viscosity of the component with the lower and the higher viscosity, respectively [Pa · s] - Viscosity of a binary mixture [Pa · s] - [] The usual intrinsic viscosity of the component with the highest viscosity [cm³/g] - _r / ₁ [—] - _{sp r} – 1 [—] - [] -intrinsic viscosity [—] - [] _v Volume intrinsic viscosity [—] - [] _w Weight intrinsic viscosity [—] - [] _c Concentration intrinsic viscosity, identical to [] [cm³/g] - T _e Temperature at which the two liquids have the same viscosity [K] - _e Viscosity at temperatureT _e [Pa · s] - P ₁,P ₂ Density of the component with the lower and the higher viscosity, respectively - R Gas constant [J · Mol^–1 · K^–1]     

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     

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.     

Free convection heat transfer from an isothermal plate with arbitrary inclination

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

---

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

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

---

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

Parameter determination for thermodynamical models of the pseudoelastic behaviour of shape memory alloy by resistivity measurements and infrared thermography

Two thermodynamical models of pseudoelastic behaviour of shape memory alloys have been formulated. The first corresponds to the ideal reversible case. The second takes into account the hysteresis loop characteristic of this shape memory alloys.Two totally independent techniques are used during a loading-unloading tensile test to determine the whole set of model parameters, namely resistivity and infrared thermography measurements. In the ideal case, there is no difficulty in identifying parameters.Infrared thermography measurements are well adapted for observing the phase transformation thermal effects.Notations 1 austenite 2 martensite - ₍₎ Macroscopic infinitesimal strain tensor of phase - ₍₂₎ ^f Traceless strain tensor associated with the formation of martensite phase - Macroscopic infiniesimal strain tensor - Macroscopic infinitesimal strain tensor deviator - _f Trace - Equivalent strain - ^pe Macroscopic pseudoelastic strain tensor - x Distortion due to parent (austenite =1)product (martensite =2) phase transformation (traceless symmetric second order tensor) - M Total mass of a system - M() Total mass of phase - V Total volume of a system - V() Total volume of phase - z=M(2)/M Weight fraction of martensite - 1-z=M(1)/M Weight fraction of austenite - u ₀ ^* () Specific internal energy of phase (=1,2) - s ₀ ^* () Specific internal entropy of phase - Specific configurational energy - Specific configurational entropy - ₀ ^f (T) Driving force for temperature-induced martensitic transformation at stress free state ( ₀ ^f T) = T ^*–Ts ^*) - Kirchhoff stress tensor - Kirchhoff stress tensor deviator - Equivalent stress - Cauchy stress tensor - Mass density - K Bulk moduli (K ₀=K) - L Elastic moduli tensor (order 4) - E Young modulus - Energetic shear (₀ = ) - Poisson coefficient - M _s ^o (M _F ^o ) Martensite start (finish) temperature at stress free state - A _s ^o (A _F ^o ) Austenite start (finish) temperature at stress free state - C _v Specific heat at constant volume - k Conductivity - Pseudoelastic strain obtained in tensile test after complete phase transformation (AM) (unidimensional test) - ₀ Thermal expansion tensor - r Resistivity - 1MPa 10⁶ N/m ² - ⁽⁾ Specific free energy of phase - _n Specific free energy at non equilibrium (R model) - _n ^eq Specific free energy at equilibrium (R model) - _n ^v Volumic part of ^eq - Specific free energy at non equilibrium (R _L model) - _conf Specific coherency energy (R _L model) - _c Specific free energy at constrained equilibria (R _L model) - _it (T) Coherency term (R _L model)     

Two-phase geothermal flows with conduction and the connection with Buckley-Leverett theory

Two-phase flows of boiling water and steam in geothermal reservoirs satisfy a pair of conservation equations for mass and energy which can be combined to yield a hyperbolic wave equation for liquid saturation changes. Recent work has established that in the absence of conduction, the geothermal saturation equation is, under certain conditions, asymptotically identical with the Buckley-Leverett equation of oil recovery theory. Here we summarise this work and show that it may be extended to include conduction. In addition we show that the geothermal saturation wave speed is under all conditions formally identical with the Buckley-Leverett wave speed when the latter is written as the saturation derivative of a volumetric flow.Roman Letters C(P, S,q) geothermal saturation wave speed [ms^–1] (14) - c _t (P, S) two-phase compressibility [Pa^–1] (10) - D(P, S) diffusivity [m s^–2] (8) - E(P, S) energy density accumulation [J m^–3] (3) - g gravitational acceleration (positive downwards) [ms^–2] - h _w (P),h _w (P) specific enthalpies [J kg^–1] - J _M (P, S,P) mass flow [kg m^–2 s^–1] (5) - J _E (P, S,P) energy flow [J m^–2s^–1] (5) - k absolute permeability (constant) [m²] - k _w (S),k _s (S) relative permeabilities of liquid and vapour phases - K formation thermal conductivity (constant) [Wm^–1 K^–1] - L lower sheetC<0 in flow plane - m, c gradient and intercept - M(P, S) mass density accumulation [kg m^–3] (3) - O flow plane origin - P(x,t) pressure (primary dependent variable) [Pa] - q volume flow [ms^–1] (6) - S(x, t) liquid saturation (primary dependent variable) - S ^*(x,t) normalised saturation (Appendix) - t time (primary independent variable) [s] - T temperature (degrees Kelvin) [K] - T _sat(P) saturation line temperature [K] - TdT _sat/dP saturation line temperature derivative [K Pa^–1] (4) - T _c ,T _D convective and diffusive time constants [s] - u _w (P),u _s (P),u _r (P) specific internal energies [J kg^–1] - U upper sheetC > 0 in flow plane - U(x,t) shock velocity [m s^–1] - x spatial position (primary independent variable) [m] - X representative length - x, y flow plane coordinates - z depth variable (+z vertically downwards) [m] Greek Letters _P , _S remainder terms [Pa s^–1], [s^–1] - double-valued saturation region in the flow plane - h =h _s –h _w latent heat [J kg^–1] - = _w – _s density difference [kg m^–3] - line envelope - =D _K /D ₀ diffusivity ratio - porosity (constant) - _w (P), _s (P), _t (P, S) dynamic viscosities [Pa s] - v _w (P),v _s (P) kinematic viscosities [m²s^–1] - v ₀ =kh/KT kinematic viscosity constant [m² s^–1] - ₀ =v ₀ dynamic viscosity constant [m² s^–1] - _w (P), _s (P) density [kg m^–3] Suffixes r rock matrix - s steam (vapour) - w water (liquid) - t total - av average - 0 without conduction - K with conduction     

Existence and Uniqueness of Time-Periodic Physically Reasonable Navier-Stokes Flow Past a Body

Let be a three-dimensional exterior domain of class C^2,, 0<<1. Assume that a Navier-Stokes liquid is moving in under the action of a body force F that is time-periodic of period T, and that the velocity of the liquid is zero at spatial infinity. In this paper we show that, if F satisfies suitable conditions, and its norm, in appropriate function spaces, is sufficiently small, there is at least one time-periodic strong solution. Furthermore, the velocity field v of such a solution decays to zero for large |x| as |x|^–1 and its spatial gradient decays as |x|^–2, both uniformly in time. In addition, the pressure p decays like |x|^–2 and its gradient like |x|^–3, for almost all t[0,T]. In the special case where F is time-independent, these solutions are also time-independent and coincide with Finns physically reasonable solutions [4]. Moreover, we show that our time-periodic solutions are unique in a very large class, namely, the class of time-periodic weak solutions satisfying the energy inequality and with corresponding pressure fields verifying mild summability conditions in ×[0,T].     

Effects of MHD free convection and mass transfer on the flow past an oscillating infinite coaxial vertical circular cylinder

The effects of MHD free convection and mass transfer are taken into account on the flow past oscillating infinite coaxial vertical circular cylinder. The analytical expressions for velocity, temperature and concentration of the fluid are obtained by using perturbation technique.

---

Einwirkungen von freier MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden unendlichen koaxialen vertikalen Zylinder

Zusammenfassung Die Einwirkungen der freien MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden, unendlichen, koaxialen, vertikalen Zylinder wurden untersucht. Die analytischen Ausdrücke der Geschwindigkeit, Temperatur und Fluidkonzentration sind durch die Perturbationstechnik erhalten worden.

---

Nomenclature C _p Specific heat at constant temperature - C the species concentration near the circular cylinder - C _w the species concentration of the circular cylinder - C the species concentration of the fluid at infinite - * dimensionless species concentration - D chemical molecular diffusivity - g acceleration due to gravity - Gr Grashof number - Gm modified Grashof number - K thermal conductivity - Pr Prandtl number - r _a ,r _b radius of inner and outer cylinder - a, b dimensionless inner and outer radius - r,r coordinate and dimensionless coordinate normal to the circular cylinder - Sc Schmidt number - t time - t dimensionless time - T temperature of the fluid near the circular cylinder - T _w temperature of the circular cylinder - T temperature of the fluid at infinite - u velocity of the fluid - u dimensionless velocity of the fluid - U ₀ reference velocity - z,z coordinate and dimensionless coordinate along the circular cylinder - coefficient of volume expansion - * coefficient of thermal expansion with concentration - dimensionless temperature - H ₀ magnetic field intensity - coefficient of viscosity - _e permeability (magnetic) - kinematic viscosity - electric conductivity - density - M Hartmann number - dimensionless skin-friction - frequency - dimensionless frequency