Rheology of molten polymers

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

---

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

     

Evolution of the Flow Field around a Circular Cylinder and a Sphere upon Instantaneous Start with a Supersonic Velocity

Evolution of the flow field around a circular cylinder and a sphere instantaneously starting with a constant supersonic velocity (M = 5 and Re = 10⁵) from the state at rest is studied by means of numerical integration of unsteady twodimensional Navier–Stokes equations.     

Unsteady forced convection heat/mass transfer from a flat plate

Numerical methods are used to investigate the transient, forced convection heat/mass transfer from a finite flat plate to a steady stream of viscous, incompressible fluid. The temperature/concentration inside the plate is considered uniform. The heat/mass balance equations were solved in elliptic cylindrical coordinates by a finite difference implicit ADI method. These solutions span the parameter ranges 10 Re 400 and 0.1 Pr 10. The computations were focused on the influence of the product (aspect ratio) × (volume heat capacity ratio/Henry number) on the heat/mass transfer rate. The occurrence on the plates surface of heat/mass wake phenomena was also studied.     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

Electrode layers in flows of an inviscid plasma with good conductivity

The structure of the electromagnetic electrode layers that are produced in flows across a magnetic field by a completely ionized and inviscid plasma with good conductivity and a high magnetic Reynolds number is examined in a linear approximation. Flow past a corrugated wall and flow in a plane channel of slowly varying cross section with segmented electrodes are taken as specific examples. The possibility is demonstrated of the formation of nondissipative electrode layers with thicknesses on the order of the Debye distance or electron Larmor radius and of dissipative layers with thicknesses on the order of the skin thickness, as calculated from the diffusion rate in a magnetic field [2].In plasma flow in a transverse magnetic field, near the walls, along with the gasdynamie boundary layers, which owe their formation to viscosity, thermal conductivity, etc. (because of the presence of electromagnetic fields, their structures may vary considerably from that of ordinary gasdynamic layers), proper electromagnetic boundary layers may also be produced. An example of such layers is the Debye layer in which the quasi-neutrality of the plasma is upset. No less important, in a number of cases, is the quasi-neutral electromagnetic boundary layer, in which there is an abrupt change in the frozen-in parameter k=B/p (B is the magnetic field and p is the density of the medium). This layer plays a special role when we must explicitly allow for the Hall effect and the related formation of a longitudinal electric field (in the direction of the veloeiryv of the medium). We will call this the magnetic layer. The magnetic boundary layer can be dissipative as well as noudissipative (see below). The dissipative magnetic layer has been examined in a number of papers: for an incompressible medium with a given motion law in [1], for a compressible medium with good conductivity in [2], and with poor conductivity in [3]. In the present paper, particular attention will be devoted to nondissipative magnetic boundary layers.     

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

Dielectric properties of the composite system polyurethane — sodium chloride

Summary The dielectric properties of the composite system polyurethane-sodium chloride have been measured at frequencies between 10^–4 Hz and 3 · 10⁵ Hz in the temperature range from –150 °C up to +90 dgC. Three dielectric loss mechanisms have been found; they are indicated by ₁, ₂ and. The activation energy of the ₁-transition is 35 kcal/mole, that of the-transition 6.7 kcal/mole. The ₂-loss peak was only observed at frequencies of 10³ Hz and higher, forming one broad peak with the ₁-loss peak at lower frequencies. In the composite materials, the- and ₂-loss peaks measured at fixed frequencies were found at the same temperature. The ₂-loss peak, however, was shifted to a lower temperature, due to the sodium chloride filler. Comparison of experimental data of and tan with theoretical predictions concerning the dielectric properties of composite systems showed only partial agreement. The difference mainly consisted in. the temperature shift in the tan-peak of the ₁-transition.

---

Zusammenfassung Die dielektrischen Eigenschaften des Verbundssystems Kochsalz-Polyurethankautschuk wurden im Frequenzgebiet zwischen 10^–4 Hz und 3.10⁵ Hz und im Temperaturbereich von –150 °C bis +90 °C gemessen. Es wurden drei dielektrische Verlustmaxima gefunden, die mit ₁, ₂ und angedeutet werden. Die Aktivierungsenergie des ₁-Überganges beträgt 35 kcal/Mol, die des-Überganges 6.7 kcal/Mol. Das ₂-Maximum konnte nur bei Frequenzen höher als 10³Hz vom ₁-Maximum gesondert erfaßt werden. Die Lage der ₂- und-Maxima war vom Füllgrad unabhängig. Das ₁-Maximum verschiebt sich mit steigendem Füllgrad zu niedrigeren Temperaturen. Die gemessenen Werte von und tan stimmen nur teilweise mit den Aussagen einer Theorie der dielektrischen Eigenschaften von Mischkörpern überein.

     

Approximate relationships for the determination of the pressure on the surface of a sphere or cylinder with an arbitrary mach number in the free stream

Numerical calculations have been made [1–4] of the pressure distribution over the surface of a sphere or cylinder during transverse flow in the range 0 /2, where is the angle reckoned from the stagnation point along the meridional plane, and on the basis of these results simple analytical equations have been proposed in order to determine the pressure for arbitrary Mach numbers M in the free stream. The gas is assumed to be ideal and perfect.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 185–188, March–April, 1985.     

Calculation of heat transfer accompanying hypersonic three-dimensional flow of air over slender blunt-nosed cones

The effective length method [1, 2] has been used to make systematic calculations of the heat transfer for laminar and turbulent boundary layers on slender blunt-nosed cones at small angles of attack ( + 5° in a separationless hypersonic air stream dissociating in equilibrium (half-angles of the cones 0 20°, angles of attack 0 15°, Mach numbers 5 M 25). The parameters of the gas at the outer edge of the boundary layer were taken equal to the inviscid parameters on the surface of the cones. Analysis of the results leads to simple approximate dependences for the heat transfer coefficients.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1981.     

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

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

Flow in porous media III: Deformable media

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

Base-pressure pulsations behind a cylinder and disk in a subsonic stream flow

Base-pressure fluctuations behind a long cylinder (l/d 5–10) and the disk (l/ d 0.0) is investigated experimentally in this paper. The spectral and correlation characteristics of the base-pressure fluctuations behind axisymmetric bodies at a Mach number M 1.0 are generalized on the basis of the data obtained and the results of other authors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 181–183, January–February, 1977.     

Stability and Transition on a Swept Cylinder in a Supersonic Flow

Results of experimental investigations of the evolution of natural disturbances and laminar–turbulent transition in a supersonic boundary layer on the attachment line of a circular cylinder with a sweep angle of 68° and a freestream Mach number M = 2 are presented. The experimental studies are supplemented by calculations of the mean flow and stability characteristics. Flow regimes in the boundary layer on the attachment line are determined by a hotwire technique as functions of the Reynolds number and height of twodimensional roughness elements. The results are compared with NASA (Ames) experiments.     

Asymptotic Nusselt and entropy generation numbers for non-Newtonian fluid flow

The first and the second laws of thermodynamics are applied to obtain general expressions for the asymptotic Nusselt (Nu) and entropy generation (Ns) numbers in ducts of circular cross-section or a channel made of two long parallel plates. The governing momentum and energy equations are simplified with reasonable assumptions and solved to obtain closed form of analytical solutions. Results of the asymptotic Nusselt and Bejan numbers are presented graphically as a function of fluid index (n), Peclet number (Pe), and group parameter (Br/).     

Application of the Chapman-Enskog method to the case of a binary two-temperature gas mixture

An interesting property of the flows of a binary mixture of neutral gases for which the molecular mass ratio =m/M1 is that within the limits of the applicability of continuum mechanics the components of the mixture may have different temperatures. The process of establishing the Maxwellian equilibrium state in such a mixture divides into several stages, which are characterized by relaxation times _i which differ in order of magnitude. First the state of the light component reaches equilibrium, then the heavy component, after which equilibrium between the components is established [1]. In the simplest case the relaxation times differ from one another by a factor of ^*.Here the mixture component temperature difference relaxation time _T /, where is the relaxation time for the light component. If 1, 1, so that _T ~1, then for the characteristic hydrodynamic time scale t~1 the relative temperature difference will be of order unity. In the absence of strong external force fields the component velocity difference is negligibly small, since its relaxation time _vt1.In the case of a fully ionized plasma the Chapman-Enskog method is quite easily extended to the case of the two-temperature mixture [3], since the Landau collision integral is used, which decomposes directly with respect to . In the Boltzmann cross collision integral, the quantity appears in the formulas relating the velocities before and after collision, which hinders the decomposition of this integral with respect to , which is necessary for calculating the relaxation terms in the equations for temperatures differing from zero in the Euler approximation [4] (the transport coefficients are calculated considerably more simply, since for their determination it is sufficient to account for only the first (Lorentzian [5]) terms of the decomposition of the cross collision integrals with respect to ). This led to the use in [4] for obtaining the equations of the considered continuum mixture of a specially constructed model kinetic equation (of the Bhatnagar-Krook type) which has an undetermined degree of accuracy.In the following we use the Boltzmann equations to obtain the equations of motion of a two-temperature binary gas mixture in an approximation analogous to that of Navier-Stokes (for convenience we shall term this approximation the Navier-Stokes approximation) to determine the transport coefficients and the relaxation terms of the equations for the temperatures. The equations in the Burnett approximation, and so on, may be obtained similarly, although this derivation is not useful in practice.     

Deformation to breaking of deep water gravity waves

The results of laboratory observations of the deformation of deep water gravity waves leading to wave breaking are reported. The specially developed visualization technique which was used is described. A preliminary analysis of the results has led to similar conclusions than recently developed theories. As a main fact, the observed wave breaking appears as the result of, first, a modulational instability which causes the local wave steepness to approach a maximum and, second, a rapidly growing instability leading directly to the breaking.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness =H/L - crest steepness =/L - trough steepness =/L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness =/F ₁ - crest rear steepness =/F ₂ - vertical asymmetry factor=F ₂/F ₁ (describing the wave asymmetry with respect to a vertical axis through the wave crest) - µ horizontal asymmetry factor=/H (describing the wave asymmetry with respect to a horizontal axis: SWL) - T ₀ wavemaker period - L ₀ theoretical wave length of a small amplitude sinusoïdal wave generated at T _inf0 ^sup–1 frequency - ₀ average wave height     

Hartmann Effect. Region of Existence and Oscillation Frequencies

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

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

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

Group Classification of Equations of the Form y″ = f(x,y)

The problem of classification of ordinary differential equations of the form y = f(x,y) by admissible local Lie groups of transformations is solved. Standard equations are listed on the basis of the equivalence concept. The classes of equations admitting a oneparameter group and obtained from the standard equations by invariant extension are described.     

Transverse Decay of Solutions in an Elastic Cylinder

In a semiinfinite cylinder composed of anisotropic linearised elastic material, loaded on the base and clamped along the lateral surface, it is known that the solution as measured, for example, by the strainenergy flux through a plane crossection decays longitudinally at most exponentially with respect to the axial distance from the base. There is, however, also a transverse radial decay of the solution, again measured for example by the strainenergy, occurring from the region close to the cylinder's axis to the region near the lateral surface, where the energy vanishes.This problem is considered in the present paper which discusses a circular semiinfinite cylinder and derives an estimate for the strainenergy contained in a cylindrical annulus at a given distance from the base and of variable height, and whose outer surface coincides with the lateral surface of the cylinder. It is shown that the strainenergy decays at most algebraically to zero as the inner radius of the annulus increases to that of the cylinder.Sommario. E'noto che in un cilindro semiinfinito composto da materiale elastico lineare anisotropo, caricato sulla base ed incastrato lungo la superficie laterale, la soluzione elastica, misurata, per esempio, dal flusso di energia di deformazione attraverso una sezione trasversale piana, decade con legge al più esponenziale con la distanza dalla base. C'è tuttavia, anche un decadimento radiale della soluzione, misurato, per esempio, dall'energia di deformazione che passa dalla regione vicina all'asse del cilindro a quella vicin alla superficie laterale dove l'energia si annullaQuesto problema è qui studiato. Si discute in particolare un cilindro circolare semiinfinito e si deduce una stima per l'energia di deformazione contenuta in un anello cilindrico ad una distanza assegnata dalla base e di altezza variabile, e la cui superficie esterna coincide con la superficie laterale del cilindro. Si dimostra che l'energia di deformazione decade al più con legge algebrica a zero quando il raggio interno del cilindro si avvicina a quello esterno.