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.     

Effect of Streamwise Streaky Structures on Turbulization of a Circular Jet

Experimental investigations of the influence of streamwise streaky structures on turbulization of a circular laminar jet are described. The qualitative characteristics of jet evolution are studied by smoke visualization of the flow pattern in the jet and by filming the transverse and longitudinal sections of the jet illuminated by the laser sheet with image stroboscopy. It is shown that the streaky structures can be generated directly at the nozzle exit, and their interaction with the Kelvin–Helmholtz ring vortices leads to emergence of azimuthal beams ( structures) by a mechanism similar to threedimensional distortion of the twodimensional Tollmien–Schlichting wave at the nonlinear stage of the classical transition in nearwall flows. The effect of the jetexhaustion velocity and acoustic action on jet turbulization is considered.     

LDA measurements in low Reynolds number turbulent boundary layer

LDA measurements of the mean velocity in a low Reynolds number turbulent boundary layer allow a direct estimate of the friction velocity U from the value of /y at the wall. The trend of the Reynolds number dependence of / is similar to the direct numerical simulations of Spalart (1988).     

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.     

Transport in ordered and disordered porous media III: Closure and comparison between theory and experiment

In this paper we examine the closure problem associated with the volume averaged form of the Stokes equations presented in Part II. For both ordered and disordered porous media, we make use of a spatially periodic model of a porous medium. Under these circumstances the closure problem, in terms of theclosure variables, is independent of the weighting functions used in the spatial smoothing process. Comparison between theory and experiment suggests that the geometrical characteristics of the unit cell dominate the calculated value of the Darcy's law permeability tensor, whereas the periodic conditions required for thelocal form of the closure problem play only a minor role.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface associated with the local closure problem, m² - A _p surface area of a particle, m² - b vector used to represent the pressure deviation, m^–1 - B ⁰ B+I, a second order tensor that maps v _m ontov - B second-order tensor used to represent the velocity deviation - d _p 6V _p/A_p, effective particle diameter, m - d a vector related to the pressure, m - D a second-order tensor related to the velocity, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor calculated on the basis of a spatially periodic model, m² - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p characteristic length for the volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - characteristic length (pore scale) for the-phase - _i i=1, 2, 3 lattice vectors, m - weighting function - m(-y) , convolution product weighting function - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - n unit normal vector pointing from the-phase toward the -phase - p pressure in the-phase, N/m² - p _m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p – p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function - r position vector, m - r position vector locating points in the-phase, m. - V averaging volume, m³ - B volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v _m superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - v traditional superficial volume averaged velocity, m/s - v – v _m , spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the -phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * , weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Theory of plane flows of a readily conducting plasma in a pipe

Simplified equations are obtained describing slowly changing plane flows of a readily conducting quasineutral inviscid plasma in a pipe. The practically interesting case of flow in a channel with solid metal ideally conducting walls (electrodes) is analyzed. When the gas pressure is large by comparison with the magnetic pressure ( 1), the field and current distribution is determined by gas dynamic factors, and the solid electrodes perturb the longitudinal electric field in a skin of the flow, symmetrically on the two sides of the flow, leading to attenuation of the longitudinal electric field near the input to the pipe; we also consider problems in the motion of the plasma under ideal and under poor conductivity. In the converse limiting case ( 1), it is shown that as the motion of the plasma in the pipe accelerates near the anode, there is observed an increase in the intensity of the electric field which is sharply inhomogeneous in the transverse direction. The possibility of the plasma breaking away from the anode (the limiting regime) is indicated, this being accompanied by a divergence between the electron velocity and the velocity of the ions. A criterion is obtained for the breakaway of the plasma, and its possible connection with the occurrence of pre-anode explosions is noted. It is shown that for 1, Joule losses are small by comparison with the power in the charge and the magnitude of the losses is independent of the conductivity of the plasma.Translated from Zhurnal Prikladnoi Mekhaniki i Teknicheskoi Fiziki, No. 4, pp. 9–19, July–August, 1970.     

Messung des ersten Normalspannungskoeffizienten von Kunststoffschmelzen im Bogenspalt

Zusammenfassung Bei einer stationären Schichtenströmung in einem Bogenspalt (azimutale Druckströmung im Ringspalt) bildet sich zwischen Innen- und Außenwand eine Druckdifferenz aus, deren Größe ein Maß für den 1. Normalspannungskoeffizienten der elastischen Flüssigkeit im Spalt ist. Die Strömung läßt sich zur Messung des 1. Normalspannungskoeffizienten verwenden. Der Schergeschwindigkeitsbereich der Messung liegt, wie bei der Kapillarrheometrie zur Bestimmung der Viskosität, zwischen 1 und 1000 s^–1. Die Auswertung der Messungen ist wegen des inhomogenen Scherfeldes relativ kompliziert. In der Arbeit wird ein besonders wirkungsvolles numerisches Auswerteverfahren hergeleitet und auf bestehende Messungen angewendet. Eine Besonderheit des Auswerteverfahrens ist die Freiheit der Wahl des Approximationsansatzes für die Viskositätskurve, während analytische Verfahren meist an einen bestimmten Ansatz gebunden sind. Außerdem braucht, im Gegensatz zu anderen derartigen Verfahren, die Position des schubspannungsfreien Stromfadensr ₀ nicht bestimmt zu werden.

---

Summary The stress in steady viscometric flow of molten polymers is determined by the viscosity and by the two normal stress coefficients ₁ and ₂. The paper describes a method of measuring ₁ by means of steady circumferential shear flow in an annulus. The cylinders are stationary and the fluid flows due to a circumferential pressure gradient. The radial normal stresses at the outer and at the inner wall are different from each other. The pressure-differencep is a measure for the 1. normal stress coefficient of the viscoelastic fluid. Due to the inhomogeneous shear field, the evaluation of ₁ fromp measurements is quite complicated. A powerful numerical method of evaluation has been developed and applied to existing data. The method is not restricted to a special empirical formula for the flow curve (as an analytical method would be) and does not require the knowledge of the positionr ₀ of the stress-free stream line.

---

a Pa s² Stoffparameter des Ansatzes des 1. Normalspannungskoeffizienten, s. Gl. [8] - AR _i — Koeffizient des Druckgefälles in-Richtung (Programm PFEIL) - AU _i — Koeffizient für Integration nach Simpson-Regel (Programm PFEIL) - b s² Stoffparameter des Ansatzes des 1. Normalspannungskoeffizienten - B _i — Koeffizient auf der rechten Seite des linearen Gleichungssystems (Programm PFEIL) - c — Exponent des Ansatzes des 1. Normalspannungskoeffizienten - CL _i CM _i CR _i — Koeffizienten der dimensionslosen Geschwindigkeit in dem linearen Gleichungssystem (Programm PFEIL) - F ₁,F ₂,F ₃ — Ableitungen der Summe der Fehlerquadrate nacha, b undc - G _k — Gewichtsfaktor - h m Spaltweite,r _a –r _i - H — dimensionslose Spaltweite, (r _a –r _i )/r _a - l m Länge des Bogenspaltes, 0,75(r _a +r _i ) - m — Exponent des Potenzansatzes, s. Gl. [13] - n — Dämpfungskonstante - N ₁ Pa 1. Normalspannungsdifferenz, _rr – - N ₂ Pa 2. Normalspannungsdifferenz - p Pa Druck - p Pa Druckgradient in-Richtung - P — dimensionsloser Druckgradient in-Richtung, s. Gl. [14] - p, p _k Pa Normalspannungsdifferenz zwischen Innen- und Außenwand im Bogenspalt, (– p + _rr ) _a – (–p + _rr ) _i - Q — Summe der Fehlerquadrate - r, R= r/r _a m, — Radiusvektor (Koordinate in Gradientenrichtung) - r ₀,R ₀=r ₀/r _a m, — Radius des neutralen Fadens - R — dimensionslose radiale Schrittweite - T, °C Temperatur bzw. Bezugstemperatur - v ms^–1 Geschwindigkeitskomponente in-Richtung - V ,V _,i — dimensionslose Geschwindigkeitskomponente in-Richtung - V _a ,V _k — dimensionslose Geschwindigkeit an der Außen- bzw. Innenwand - v _r ,v _z ms^–1 Geschwindigkeitskomponenten inr-undz-Richtung - ms ^–1 mittlere Geschwindigkeit in-Richtung - z m Koordinate in der indifferenten Richtung - K^–1 Temperaturkoeffizient der Viskosität - s^–1 Schergeschwindigkeit - s^–1 kritische Schergeschwindigkeit der Viskositätskurve, s. Gl. [13] - s^–1 Bezugsschergeschwindigkeit, - — dimensionslose Schergeschwindigkeit - — dimensionslose kritische Schergeschwindigkeit, - Pa s Viskosität - ₀ Pa s Nullviskosität - Pa s Bezugsviskosität, - — Radienverhältnis,r _i /r _a - ₁ Pa s ² 1. Normalspannungskoeffizient - Pa s² mittlerer 1. Normalspannungskoeffizient - ₂ Pa s² 2. Normalspannungskoeffizient - — Koordinate in Strömungsrichtung - Pa Spannung - a an der Außenwand - i, an der Innenwand - i laufender Index inr-Richtung - k Nummer des Meßpunktes - n Anzahl der Meßpunkte - n _i nord für Programm PFEIL - s _i süd für Programm PFEIL Mit 9 Abbildungen und 2 Tabellen     

RADAR ECHO DELAY IN THE(Ω,Aμν)-FIELD THEORY

The present work aims to consider the.fourth test of general relativity theory by Shapiro.using radar echo delay in Yu’s(Ω,A_μν)-field theory.     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters 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² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - 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 and radius of a capillary tube, m - 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² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

Waves of finite amplitude in a hot plasma

A study is made of a plane shock wave of arbitrary strength propagating in a hot rarefied plasma across the magnetic field. The question of the propagation of nonstationary waves of finite but small amplitude under these conditions is examined.Fairly detailed studies have been made of waves of finite amplitude in a cold rarefied plasma. The profile of such waves is formed as the result of nonlinear and dispersion effects, the dispersion effects being caused by electron inertia and plasma anisotropy. If the gas-kinetic pressure of the plasma is taken into account, then dispersion effects appear which are associated with the fact that the Larmor radius of the ions is finite. Stationary waves of small but finite amplitude propagating across the magnetic field in a hot plasma (when the gas-kinetic pressure p is comparable with the magnetic pressure H²/87) have been treated in [1, 2]. In [1] an isolated rarefaction wave was found in a hot plasma, instead of the compression wave characteristic of a cold plasma, and a qualitative picture of the shock wave structure was given. In [2] a study was made of a small-amplitude shock wave with the finite size of the ion Larmor radius taken into account. The present paper investigates the structure of shock waves of arbitrary strength which propagate across the magnetic field in a fairly hot rarefied plasma, and also examines nonstationary waves of finite but small amplitude excited in a plasma by a magnetic piston acting over a limited time interval.Notation p gas-kinetic pressure - H magnetic field - u, v macroscopic velocities along the x and y axes - density - m_e(m_i) mass of electron (ion) - plasma conductivity - _H ion-cyclotron frequency - V_A Alfvèn velocity - c velocity of light - adiabatic exponent - V specific volume - _0e(_0i) electron (ion) plasma frequency - S₀ velocity of sound. In conclusion the author thanks R. Z. Sagdeev and N. N. Yanenko for discussing the paper, and also R. N, Makarov for helping with the numerical computations.     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

Correlations between porous medium flow data and turbulent tube flow results for aqueous hydroxypropylguar solutions (Part 2)

Turbulent tube flow and the flow through a porous medium of aqueous hydroxypropylguar (HPG) solutions in concentrations from 100 wppm to 5000 wppm is investigated. Taking the rheological flow curves into account reveals that the effectiveness in turbulent tube flow and the efficiency for the flow through a porous medium both start at the same onset wall shear stress of 1.3 Pa. The similarity of the curves = ( _w ) and = ( _w ), respectively, leads to a simple linear relation / =k, where the constantk or proportionality depends uponc. This offers the possibility to deduce (for turbulent tube flow) from (for flow through a porous medium). In conjunction with rheological data, will reveal whether, and if yes to what extent, drag reduction will take place (even at high concentrations).The relation of our treatment to the model-based Deborah number concept is shown and a scale-up formula for the onset in turbulent tube flow is deduced as well.     

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.     

Der Spannungs- und Verschiebungszustand drehsymmetrischer Membrane mit beliebig gekrümmter Meridiankurve

Zusammenfassung Die Einführung von Zylinderkoordinaten (x, r, ) in die Gleichgewichtsbedingungen der Schnittkräfte bzw. in die Beziehungen zwischen Verzerrung und Verschiebungen am differentialen Schalenabschnitt ermöglicht die Berechnung des Spannungs- und Verschiebungszustandes von drehsymmetrischen Membranen mit beliebig gekrümmter Meridiankurve auf die Integration einer einfachen, linearen partiellen Differentialgleichung zweiter Ordnung für eine charakteristische FunktionF bzw. zurückzuführen. Eine geschlossene Lösung und damit eine Darstellung der Schnittkräfte und Verschiebungen durch explizite Formeln ist bei harmonischer Belastung cosn für zwei Funktionsgruppen=x ² und=x ^–3 möglich. Im Sonderfall der drehsymmetrischen und der antimetrischen Belastung mitn=0 undn=1 gelten die Gleichungen der Schnitt- und Verschiebungsgrößen für eine beliebige Meridianfunktion=(). Die Betrachtungen der Randbedingungen offener Schalen bei harmonischer Belastung geben über die infinitesimalen Deformationen einer drehsymmetrischen Membran mit überall negativer Krümmung Aufschluß.     

Exact Solutions of the Hydrodynamic Equations Derived from Partially Invariant Solutions

The paper proposes a heuristic approach to constructing exact solutions of the hydrodynamic equations based on the specificity of these equations. A number of systems of hydrodynamic equations possess the following structure: they contain a reduced system of n equations and an additional equation for an extra function w. In this case, the reduced system, in which w = 0, admits a Lie group G. Taking a certain partially invariant solution of the reduced system with respect to this group as a seed:rdquo; solution, we can find a solution of the entire system, in which the functional dependence of the invariant part of the seed solution on the invariants of the group G has the previous form. Implementation of the algorithm proposed is exemplified by constructing new exact solutions of the equations of rotationally symmetric motion of an ideal incompressible liquid and the equations of concentrational convection in a plane boundary layer and thermal convection in a rotating layer of a viscous liquid.     

Numerical computation of high Rayleigh number natural convection and prediction of hot radiator induced room air motion

The two-dimensional stationary turbulent buoyant flow and heat transfer in a cavity at high Rayleigh numbers was computed numerically. The k– turbulence model was used. The time-averaged equations for momentum, energy and continuity, which are coupled to the turbulence equations, were solved using a finite difference formulation. In order to validate the computer code, a comparison exercise was carried out. The test results are in good agreement with the internationally accepted benchmark solution. Grid-refinement shows the necessity of a very fine grid at high Rayleigh numbers with especially small grid-distances in the near-wall region. The computed boundary layer velocity profiles are in excellent agreement with available experimental data. The local heat transfer in the turbulent part of the boundary layers is predicted 20% too high. Computations were carried out for the natural convective flow in a room induced by a hot radiator and a cold window. Various radiator configurations and types of thermal boundary conditions were applied including thermal radiation interaction between surfaces.Nomenclature a thermal diffusivity (m²/s) - C constant in _t expression - D cavity dimensions (m) - g acceleration of gravity (m/s²) - G _k production/destruction of k by buoyancy (kg/ms³) - h enthalpy (J/kg) - IX index of grid point - k turbulent kinetic energy (m²/s²) - m dimensionless stratification parameter - Nu overall Nusselt number - Nu _y local Nusselt number - NX total number of grid points - p pressure (N/m²) - P _k production of k by shear stress (kg/ms³) - Q heat flux through wall (W/m) - Ra overall Rayleigh number - Ra _y local Rayleigh number - Re _t turbulent Reynolds number - S source term in -equation (kg/ms⁴) - S source term for - T _c, T _h temperatures of cold and hot walls (K) - T _s (y) stratification temperature on vertical mid-line (K) - T ₀ mean cavity temperature (K) - u, v horizontal and vertical velocity components (m/s) - u ₀ Brunt-Vaisälä velocity scale (m/s) - x, y horizontal and vertical coordinates (m) - non-linearity parameter for grid - coefficient of thermal expansion (l/K) - jet angle (°) - diffusivity for - S dissipation rate for turbulent kinetic energy (m²/s³) - variable to be solved - thermal conductivity (W/mK) - , _t kinematic and eddy viscosities (m²/s) - stream function (kg/ms) - density (kg/m³) - _k, , _t constants in k– model     

The mirror relations and nonlinear viscoelasticity of polymer melts

An attempt is made to incorporate into a quasilinear viscoelastic constitutive equation of the Boltzmann superposition type the two mirror relations of Gleissle, as well as his relation between the steady-state first normal-stress difference and the shear viscosity curve. It is shown that the three relations can hold separately within this constitutive model, but not simultaneously, because they require a different nonlinear strain measure, namelyS ₁₂ () = – a ( – 1) (a = 0 for 1,a = 1 for 1) for the mirroring of the viscosities,S ₁₂ () = – a (–k ²/) (a = 0 for k, a = 1 for k) for the mirroring of the first normal-stress coefficients, and for the third relation. Here denotes the shear strain and erf the error function. Experimental data on melts of a low-density polyethylene, a high-density polyethylene and a polypropylene show that the mirror relations are passable approximations, but that the third relation meets reality surprisingly close if the right value ofk is used.     

Dehnungsverhalten von Polyäthylen-Schmelzen

Zusammenfassung In einem Dehnungsrheometer werden Spannungs-Dehnungs-Diagramme von Polyäthylen-Schmelzen bei 150 °C und bei konstanter Dehnungsgeschwindigkeit gemessen ( zwischen 0,001 und 1 sec^–1). Weiterhin wird der reversible (elastische) Dehnungsanteil bestimmt. Messungen mit einem Dehnungstester für Kunststoff-Schmelzen ergänzen die Ausführungen.Die Ergebnisse zeigen deutlich, daß bei Dehnung mit zunehmender Verformungsgeschwindigkeit die Dehnungsviskosität nicht abnimmt, im Gegensatz zu dem bekannten strukturviskosen Verhalten bei Scherung.Bei Dehnungen bis zu=1 kann das Verhalten unabhängig von beschrieben werden, wenn als viskoelastische Materialfunktion die Dehnungs-Spannviskosität betrachtet wird. In diesem Bereich von gilt dabei die BeziehungT(t)=3 _s (t) mit _s (t) als zeitabhängige Scherviskosität im linear-viskoelastischen Bereich.Bei größeren Dehnungen und nicht zu kleinen Dehnungsgeschwindigkeiten zeigt verzweigtes Polyäthylen eine zusätzliche starke Spannungszunahme. In dem Bereich dieser zusätzlichen Verfestigung ist das Verhalten im wesentlichen eine Funktion der Dehnung und fast unabhängig von . Die zusätzliche Verfestigung scheint eine Folge der Verzweigungsstruktur des verzweigten Polyäthylens zu sein, da bei Linear-PE ein derartiger Verlauf des Spannungs-Dehnungs-Diagramms nicht beobachtet wird.Die Betrachtung des reversiblen Dehnungsanteils _R zeigt bei der ausführlich untersuchten Schmelze I (verzweigtes PE) drei verschiedene Bereiche: Unterhalb einer Grenzdehnungsgeschwindigkeit ist _R =0, unterhalb einer Versuchszeitt ^** ist _R =. Im dazwischenliegenden Bereich treten elastische und viskose Dehnungsanteile auf,= _R + _V , wobei für niedrige gilt, daß _R lg . Die Grenze wird der Frequenz der thermisch aktivierten Platzwechsel zugeordnet,t ^** erscheint als Zeit, innerhalb der die Verhakungen wie fixierte Vernetzungen wirken.In dem Anhang wird der Einfluß der Grenzflächenspannung zwischen PE-Schmelze und Silikonöl auf die Ergebnisse der Dehnungsversuche diskutiert.

---

Summary Stress-strain relations for different PE melts are measured at 150 °C in an extensional rheometer under the condition of a constant extensional strain rate ( between 0,001 and 1 sec^–1). Further, the recoverable (elastic) portion _R of the total strain ( in Hencky's measure) is determined and additional measurements with a tensile tester for polymer melts are described.The results show clearly that in extension there is no decrease of the tensile viscosity with increasing deformation rate, in contrast to the well-known pseudoplastic behaviour in shear. Within total strains<1 the tensile behaviour can be described independently from by means of a viscoelastic material function called stressing viscosity . In this range of the relation _T (t)=3 _s (t) holds, where _s (t) is the stressing viscosity in shear in the linear viscoelastic range. For larger tensile strains and not too small branched PB shows a remarkable increase in stress. This work-hardening behaviour is mainly a function of and almost independent from . This additional hardening seems to be due to the branches in branched PE, because linear PE does not show this effect.The discussion of the recoverable tensile strain _R gives three regions of tensile rate: Below a critical there is _R =0. At times shorter thant ^** the equation _R = is valid. Within these limits both elastic and viscous portions of the total strain= _R + _V exist. may correlate with the frequency of the thermally activated position changes of the molecular segments.t ^** is assumed to be the time for the entanglements to act as fixed cross-links.In the appendix the influence of the interface tension between PE melt and silicone oil on the results of the tensile experiments is discussed.

---

---

Vorgetragen auf der Deutschen Rheologen-Tagung, Berlin, 11.-13. Mai 1970.

---

An der Weiterentwicklung des Dehnungsrheometers, an der Durchführung und Auswertung der Messungen waren die HerrenB. Kienle, F. Landmesser, M, Reuther undF. Scherr beteiligt. Herr Dr.F.Ramsteiner und HerrH. Schroeck haben sich um die Herstellung der Stränge aus Linear-PE bemüht. Herr Dr.W. Ball besorgte die GPC-Messungen und Herr Dr.P. Simak die Ultrarot-Untersuchung. Den vorgenannten Herren sei für ihre Hilfe beim Zustandekommen dieser Arbeit gedankt. Herrn Dr.H. Baur danke ich für wertvolle Diskussionen.