An Existence Theorem for Generalized von Kármán Equations

Using techniques from formal asymptotic analysis, the first two authors have recently identified generalized von Kármán equations, which constitute a two-dimensional model for a nonlinearly elastic plate where only a portion of the lateral face is subjected to boundary conditions of von Kármán's type, the remaining portion being free. In this paper, we establish an existence theorem for these equations. To this end, we first reduce them to a single equation, which generalizes a cubic operator equation introduced by M.S. Berger and P. Fife. We then directly solve this equation, notably by adapting a crucial compactness method due to J.-L. Lions. Résumé. En utilisant les techniques de l'analyse asymptotique formelle, les deux premiers auteurs ont récemment identifié des équations de von Kármán généralisées, qui constituent un modèle bi-dimensionnel de plaque non linéairement élastique dont une partie seulement de la face latérale est soumise à des conditions aux limites de von Kármán, la partie restante étant libre. Dans cet article, on établit un théorème d'existence pour ces équations. À cette fin, elles sont d'abord réduites à une seule équation, qui généralise une équation faisant intervenir un opérateur cubique, introduite par M.S. Berger et P. Fife. On résout ensuite directement cette équation, en adaptant notamment une méthode cruciale de compacité due à J.-L. Lions.     

The chaotic response of a fluttering panel: The influence of maneuvering

The influence of maneuvering on the chaotic response of a fluttering buckled plate on an aircraft has been studied. The governing equations, derived using Lagrangian mechanics, include geometric non-linearities associated with the occurrence of tensile stresses, as well as coupling between the angular velocity of the maneuver and the elastic degrees of freedom. Numerical simulation for periodic and chaotic responses are conducted in order to analyze the influence of the pull-up maneuver on the dynamic behavior of the panel. Long-time histories phase-plane plots, and power spectra of the responses are presented. As the maneuver (load factor) increases, the system exhibits complicated dynamic behavior including a direct and inverse cascade of subharmonic bifurcations, intermittency, and chaos. Beside these classical routes of transition from a periodic state to chaos, our calculations suggest amplitude modulation as a possible new mode of transition to chaos. Consequently this research contributes to the understanding of the mechanisms through which the transition between periodic and strange attractors occurs in, dissipative mechanical systems. In the case of a prescribed time dependent maneuver, a remarkable transition between the different types of limit cycles is presented.Nomenclature a plate length - a _r u _r /h - D plate bending stiffness - E modulus of elasticity - g acceleration due to gravity - h plate thickness - j₁,j₂,j₃ base vectors of the body frame of reference - K spring constant - M Mach number - n 1 + ₀/g - N ₁ applied in-plane force - p–p aerodynamic pressure - P pa ⁴/Dh - q ₀/2 - Q _r generalized Lagrangian forces - R rotation matrix - R ₄ N, a ²/D - t time - kinetic energy - u plate deflection - u displacement of the structure - u _r modal amplitude - v₀ velocity - x coordinates in the inertial frame of reference - z coordinates in the body frame of reference - Ka/(Ka+Eh) - - elastic energy - 2qa ³/D - a/_mh - Poisson's ratio - material coordinates - air density - _m plate density - - _r prescribed functions - _r sin(r z/a) - angular velocity - a/v₀ - skew-symmetric matrix form of the angular velocity     

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.     

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     

Elastico-viscous boundary-layer flow on the surface of a sphere

Summary The Prandtl boundary-layer theory is extended for an idealized elastico-viscous liquid. The boundary layer equations are solved approximately by Kármán-Pohlhausen technique for the case of a sphere. It is shown that the increase in the elasticity of the liquid causes a shift in the point of separation towards the forward stagnation point.

---

Zusammenfassung Die Prandtlsche Grenzschicht-Theorie wird für eine idealisierte viskoelastische Flüssigkeit erweitert. Die Grenzschichtgleichungen werden für den Fall einer angeströmten Kugel näherungsweise mit Hilfe der Kármán-Pohlhausen-Methode gelöst. Es wird gezeigt, daß das Anwachsen der Flüssigkeitselastizität eine Verschiebung des Ablösepunktes auf den vorderen Staupunkt hin zur Folge hat.

---

Nomenclature b ^ik arbitrary contravariant tensor - D non-dimensional boundary layer thickness - g _ik metric tensor of a fixed coordinate system - K curvature at any point on the generating curve - K ₀ elastico-viscous parameter - p arbitrary hydrostatic pressure - p _ik stress tensor - p ^ik part of stress tensor associated with the change of shape of material - R radius of the sphere - r radius of any transverse cross-section of the sphere - t time - U potential velocity around the body - U stream-velocity at a large distance from the body - u, w velocity components along (x, z) directions respectively - x distance measured along a generating line from the forward stagnation point - z distance measured along a normal to the surface - non-dimensional elastico-viscous parameter - density of the liquid - boundary layer thickness - convected time derivative - ₀ limiting viscosity for very small changes in deformation velocity - angle measured along the transverse direction - x/R - v kinematic coefficient of viscosity - T _s shearing stress on the surface of the sphere With 2 figures and 1 table     

The derivation of dual extremum and complementary stationary principles in geometrical non-linear shell theory

Summary In non-linear elasticity dual extremum principles can be formulated for some class of elastic deformations, for which uniqueness of the solution is assured. These results are used in the present paper to derive extremal variational principles for geometrical non-linear shells with moderate rotations. Furthermore two complementary variational principles are considered, which are stationary principles without any extremum property. The proposed theorems are valid also for the special cases of linear plates and shells, for the non-linear von Kármán plate theory and for non-linear Donnell-Marguerre type shells.

---

Übersicht In der nichtlinearen Elastizitätstheorie lassen sich duale Extremalprinzipe für solche elastische Verformungen herleiten, für die Eindeutigkeit der Lösung gewährleistet ist. Diese Resultate werden in der vorliegenden Arbeit benutzt, um für geometrisch-nichtlineare Schalen mit moderaten Rotationen Extremalprinzipe zu erhalten. Darüber hinaus werden zwei komplementäre Variationstheoreme angegeben, die Stationaritätsprinzipe ohne Extremaleigenschaft sind. Die vorgeschlagenen Verfahren gelten auch für die Sonderfälle der linearen Platten- und Schalentheorie, für die nichtlineare von Kármánsche Plattentheorie sowie für die nichtlineare Donnell-Marguerresche Schalentheorie.

     

Transformation and some solutions of the equations of motion of an ideal gas

The equations of one-dimensional and plane steady adiabatic motion of an ideal gas are transformed to a new form in which the role of the independent variables are played by the stream function and the function introduced by Martin [1, 2], It is shown that the function retains a constant value on a strong shock wave (and on a strong shock for plane flows). For one-dimensional isentropic motions the resulting transformation permits new exact solutions to be obtained from the exact solutions of the equations of motion. It is shown also that the one-dimensional motions of an ideal gas with the equation of state p=f(t) and the one-dimensional adiabatic motions of a gas for which p=f() are equivalent (t is time, is the stream function). It is shown that if k=s=–1, m and n are arbitrary (m+n0) and =1, the general solution of the system of equations which is fundamental in the theory of one-dimensional adiabatic self-similar motions [3] is found in parametric form with the aid of quadratures. Plane adiabatic motions of an ideal gas having the property that the pressure depends only on a single geometric coordinate are studied.     

New trends in the analysis of masonry structures

The modern theory of masonry structures has been set up on the hypothesis of no-tension behaviour, with the aim of offering a reference model, independent of materials and building techniques employed. This hypothesis gives rise to inequalities which have to be satisfied by the stress tensor components and, as a dual aspect, to the kinematic behaviour characteristics of media which can be classified as lying between solids and fluids: the structure of the masonry material consists of particles reacting elastically only when in contact. An examination of the plane-stress problem leads us to define, within the prescribed domain under admissible loads, three different subdomains with null, regular, or non-regular principal stress tensors, respectively. As the boundaries of such subdomains are not known a priori, the problem can be classified as a free boundary value problem. The analysis concerns mainly the subdomains where the stress tensor is non-regular; and a non-regularity condition det =0 is added to the equilibrium equations. This condition makes the stress problem isostatic and leads to a violation of Saint-Venant's compliance conditions on strains. Hence there is a need to introduce a strain tensor, not related to the stress tensor, which can be decomposed into an extensional component and a shearing component; we prove that such strains, of the class ^c, are similar to those of the theory of plastic flow. From the point of view of computational analysis the anelastic strains are considered as given distortions; they are computed by means of the Haar-Kármán principle, modified for computational purposes by an idea of Prager and Hodge.

---

Sommario La moderna teoria delle strutture murarie, fondata sulla rigorosa non reagenza a trazione del materiale, ha lo scopo di fornire un modello di riferimento indipendente sia dalle caratteristiche del materiale sia dalle techniche costruttive impiegate. L'ipotesi di non reagenza a trazione si traduce in disuguaglianze che le componenti del tensore di stress devono verificare; dualmente il comportamento caratteristico cinematico può esser classificato di confine, come del resto la stessa statica, tra solidi e fluidi: la struttura ipotizzata del materiale muratura consiste di particelle che reagiscono solo se sono in contatto. L'esame del problema piano porta a definire all'interno del dominio di definizione tre differenti tipi di sub-regioni in cui lo stress è nullo, canonico, o singolare. Poiché le frontiere di queste sub-regioni non sono note a priori il problema può anche essere classificato di frontiera libera. L'analisi concerne fondamentalmente la sub-regione in cui il tensore è non regolare, perché deve verificare anche la condizione det =0. Ciò rende isostatico il problema e conduce anche alla violazione della condizione di integrabilità delle deformazioni. Questo passaggio può essere superato introducendo un tensore di deformazioni a tensioni nulle che si può decomporre in una componente estensionale ed in una componente di scorrimento; si dimostra che queste deformazioni sono equivalenti a quelle che intervengono nella Teoria del flusso plastico. Dal punto di vista computazionale le deformazioni anelastiche sono considerate come distorsioni impresse determinate attraverso il principio di Haar-Kármán modificato, per le techniche computazionali, su idee di Prager e Hodge.

     

Asymptotic analysis of the nonlinear equations for an infinite,rubber-like slab under an equilibrated vertical line load

An infinite slab of incompressible Rivlin-Saunders material of constant thickness 2H is subject to an equilibrated, radially varying, vertical body force, comprising a concentrated, downward line load and a smooth, upward, exponentially distributed load with a characteristic decay length R. The deformation is axisymmetric and described by three stretches and a shear strain (or, equivalently, four strains) and a rotation which satisfy three relatively simple compatibility conditions. Force equilibrium is satisfied identically by the introduction of three stress functions. The incompressibility constraint is used to eliminate the normal stretch. With the introduction of stress-strain relations, the field equations are reduced to a set of seven, first-order, quasilinear partial differential equations. The loads, the radial distance, and the unknowns are scaled by the small parameter =H/R. As 0, 11 possible sets of field equations are found, including linear plate theory, von Kármán plate theory, Föppl membrane theory, large-strain membrane theory, and Wu's large-stretch (asymptotic) membrane theory. Notably absent as limiting cases are thick plate theories.This work was supported by the National Science Foundation under grant MSM-8618657-02.     

Kinds of local self-similarity of the velocity field of prewall turbulent flows

An analysis of dimensionalities and an approach used by Millikan [1] in analysis of mean motion are applied to investigation of the pulsational motion of three types of prewall flows of an incompressible liquid, i.e., in a boundary layer with longitudinal flow around a plate, in a round tube, and in a flat channel. It is shown that with sufficiently large Reynolds numbers there exists an interval of distances from the wall x₂, within which the integral one-point correlations and the narrow-band one-point correlations _jk do not depend on x₂. In frequency space, there exists a hyperbolic interval in which _jk=A_jku ²f^-1. Here A_jk=const; u is the dynamic velocity; and f is the frequency. It is also shown that, from the point of view of the mean motion, a distinction must be made between Kármán turbulent flow with rather large Reynolds numbers and non-Kármán flow with small, but turbulent Reynolds numbers. In the latter case, the coefficients in the logarithmic profiles of the velocity and in the law of the resistance depend on the Reynolds number. The article gives an evaluation of the Reynolds number, which can be assumed to be rather large.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 35–42, November–December, 1976.The author considers it his pleasant duty to express his indebtedness to M. A. Kashina for furnishing experimental data.     

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.

     

On the minimum of extended dissipation in viscous nematodynamics

This paper develops the variational principle of minimum extended dissipation for slow (low Reynolds number) flows of nematic liquids as described by the five parametric Leslie–Ericksen–Parodi (LEP) constitutive equations. It is shown that the Eulers equations for minimizer of the extended dissipative functional are the Stokes equations for the LEP fluid. When the molecular (including magnetic) field is absent, the extended dissipative functional coincides with the true dissipative functional, whose Euler equations are the Stokes equations for the Ericksen fluid.     

An experimental investigation of separating flow on a convex surface

The mean and turbulent characteristics of an incompressible turbulent boundary layer developing on a convex surface under the influence of an adverse pressure gradient are presented in this paper.The turbulence quantities measured include all the components of Reynolds stresses, auto-correlation functions and power spectra of the three components of turbulence. The results indicate the comparative influence of the convex curvature and adverse pressure gradient which are simultaneously acting on the flow. The investigation provides extensive experimental information which is much needed for a better understanding of turbulent shear flows.Nomenclature a, b constants in equation for velocity defect profile (Fig. 6) - c _f skin-friction coefficient (= _w/F _1/2 U ₁ ² ) - E(k ₁) one-dimensional wave number spectra - f frequency in Hz - G Clauser's equilibrium parameter = (H–1)/H(c _f /2) - H shape parameter (= ₁/ ₂) - k ₁ wave number (=2f/U) - L _u, L _v, L _w length scales of u, v and w fluctuations - p _s static pressure on the measurement surface - p _w reference tunnel wall static pressure - q ² total turbulent kinetic energy - R radius of curvature of the convex surface - R() auto-correlation function - T _u, T _v, T _w time scales of u, v and w fluctuations - U local mean velocity - U ₁ local free stream velocity - U _* friction velocity - u, v, w velocity fluctuations in x, y and z directions respectively - X streamwise coordinate measured along the surface from A (Fig. 1b) - x streamwise coordinate measured along the surface reckoned from station 9 - y coordinate normal to the surface - z spanwise coordinate - ₁/ _w · dp/dx - - boundary layer thickness - ₁ displacement thickness - ₂ momentum thickness - ₃ energy thickness - kinematic viscosity - density - time delay - _w wall shear stress     

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.     

Boundary layer effect on thermal NO decomposition behind incident shock waves

The thermal decomposition of nitric oxide (diluted in Argon) has been measured behind incident shock waves by means of IR diode laser absorption spectroscopy. In two independent runs the diode laser was tuned to the=0 =1² _3/2 R(18.5)-rotational vibrational transition and the=1 =2² _3/2 R(20.5)-rotational vibrational transition of nitric oxide, respectively. These two transitions originating from the vibrational ground state (=0) and the first excited vibrational state (=1) were selected in order to probe the homogeneity along the absorption path. The measured NO decomposition could satisfactorily be described by a chemical reaction mechanism after taking into account boundary layer corrections according to the theory of Mirels. The study forms a further proof of Mirels' theory including his prediction of the laminar-turbulent transition. It also shows, that the inhomogeneities from the boundary layer do not affect the IR linear absorption markedly.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Nonlinear viscoelastic properties and change in entanglement structure of linear polymer

Simultaneous measurements of stress relaxation and differential dynamic modulus were made at 268 K over a time scale of 10 to 10⁴⁵ s for nearly monodisperse polybutadiene (M _w =2.2x10⁵, 1,2-structure 70%, M _e =3600) and also one having coarse cross-linking (M _c =29000). Static shear strain ranged from 0.1 to 2.0. In a long-time region (t> _k ), the relaxation modulus G (; t) could be expressed by the product G (0; t) h (y). The observed h() agreed well with the Doi-Edwards theory without use of IA approximation. Both the cured and uncured samples showed initial drop of the differential storage modulus G (), ; t) followed by gradual recovery, but did not attain the value before shearing G (, ; t) for the uncured sample showed smaller values than that for the cured one in the whole measured time scale at the higher strain, confirming the two origins of nonlinear viscoelasticity of well entangled polymer; induced chain anisotropy and induced decrement in entanglement density. G (, ; t) curves for the cured sample agreed well with the BKZ predictions. But the curves for the uncured sample agreed well with the BKZ prediction only at the time scale of t< _k . BKZ prediction showed significant upward deviations at t> _k . Such the differences are discussed in terms of the two origins.Dedicated to Prof. John D. Ferry on the occasion of his 85th birthday.     

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     

On the prediction of riblet performance with engineering turbulence models

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

Pulsatile flow in a tube with wall injection and suction

This paper studies similarity solutions for pulsatile flow in a tube with wall injection and suction. The Navier-Stokes equations are reduced to a system of three ordinary differential equations. Two of the equations represent the effects of suction and injection on the steady flow while the third represents the effects of suction and injection on pulsatile flow. Since the equations for steady flow have been studied previously, the analysis centers on the third equation. This equation is solved numerically and by the method of matched asymptotic expansions. The exact numerical solutions compare well with the asymptotic solutions.The effects of suction and injection on pulsatile flow are the following: a) Small values of suction can cause a resonance-like effect for low frequency pulsatile flow. b) The annular effect still occurs but for large injection or suction the frequency at which this effect becomes dominant depends on the cross-flow Reynolds number. c) The maximum shear stress at the wall is decreased by injection, but may be increased or decreased by suction.Nomenclature a radius of the tube - a ₀ ² i ² - A₀, B₀, C₀, D₀, E₀ constant coefficients appearing in the expression for pressure - b a non-dimensionalized length - b ₀ ² i ^{2 2} - b _k complex coefficients of a power series - B - C ₁, C ₂, D complex constants - d - D _1,2 - f() F(a ^1/2)/aV - f ₀,f ₁,... functions of order one used in asymptotic expansions of f() - F(r) rv _r - g() - G(r) a steady component of velocity in axial direction - h() 4/C₀ a ² H(a ^1/2) - h ₀,h ₁,h ₂,...;l ₀,l ₁,l ₂,... functions of order one used in asymptotic expansions for h() in outer regions - H(r) complex valued function giving unsteady component of velocity - H ₀, H ₁, H ₂, ... K ₀, K ₁, K ₂, ...; L ₀, L ₁, L ₂, ... functions of order one used in asymptotic expansions for h() in inner regions - i - J ₀, J ₁, Y ₀, Y ₁ Bessel functions of first and second kind - k - K Rk/2b ² - O order symbol - p pressure - p ₁(z, t) arbitrary function related to pressure - r radial coordinate - r ₀ (1+16 ^{4 4})^1/4 - R Va/, the crossflow Reynolds number - t time - u() G(r)/V - v _r radial velocity - v _z axial velocity - V constant velocity at which fluid is injected or extracted - z axial coordinate - ² a ²/4 - 4.196 - small parameter; =–2/R (Sect. 4); =–R/2 (Sect. 5); =2/R(Sect. 6) - r ²/a ² - * 0.262 - Arctan (4 ^{2 2}) - , inner variables - kinematic viscosity - b - * zero of g() - density - (r, t) arbitrary function related to axial velocity - frequency     

Zur Berechnung des Wärmeübergangs bei einphasiger freier Konvektion in der Nähe des kritischen Zustandes

Zusammenfassung Experimentelle Ergebnisse zum überkritischen Wärmeübergang weisen für den Zustandsbereich nahe dem kritischen Punkt zum Teil große Abweichungen von der für unterkritische Fluide bekannten Abhängigkeit des Wärmeübergangskoeffizienten von der Wärmestromdichte auf. Am Beispiel des Kältemittels RC318 (C₄F₈) wird gezeigt, daß auch diese Ergebnisse mit den bekannten Beziehungen zwischen der Nußelt-, der Grashof- und der Prandtl-Zahl beschrieben werden können, wenn der thermische Ausdehnungskoeffizient und die spezifische Wärme in Gr bzw. Pr durch Differenzenquotienten ersetzt und zwei zusätzliche Parameter zur Beschreibung der Dichteänderung innerhalb der beheizten Fluidzone eingeführt werden. Da ein Teil der in den Kennzahlen benötigten Stoffwerte von RC318 im interessierenden Zustandsbereich nicht bekannt ist, werden die fehlenden Stoffwerte mit Hilfe des allgemeinen Korrespondenzprinzips berechnet.

---

Calculation of free convective heat transfer near the critical state

For certain conditions free convective heat transfer from horizontal tubes to fluids near the critical state differs widely from the well-known dependency of heat transfer coefficient from heat flux. It is shown that experiments with refrigerant RC318 (C₄F₈) even for these conditions can be described by one of the often applied relationships between Nusselt and Rayleigh numbers, if the special form of density variation within the heated region of the fluid is taken into account. Most of the thermophysical properties of RC318 being unknown near the critical state, thermodynamic similarity considerations are used to calculate these data.

---

Formelzeichen F Korrekturfaktor - R individuelle Gaskonstante - T Temperatur - Z Realfaktor - Gr, Nu, Pr Kennzahlen - a Temperaturleitzahl - c_v, C_p spezifische Wärme - d Rohrdurchmesser - g Erdbeschleunigung - h spezifische Enthalpie - m Molekülmasse - p Druck - q Wärmestromdichte - u innere Energie - v spezifisches Volumen - Wärmeübergangskoeffizient - _k Riedel-Parameter - thermischer Ausdehnungskoeffizient - Realanteil - Differenz zwischen einer Zustandsgröße des Fluids an der Heizwand und außerhalb der beheizten Zone - Asymmetrieparameter - Viskosität, dynamische - Wärmeleitzahl - Viskosität, kinematische - Dichte Indizes-hochgestellt normierte Größe - * auf den Wert am kritischen Punkt normierte Größe - 0 im Zustand des idealen Gases Indizes-tiefgestellt B, + Bezugswert - f Fluid außerhalb der beheizten Zone - k am kritischen Punkt - W an der Wand Herrn Professor Dr.-Ing. H. Glaser, Stuttgart, zum 70. Geburtstag gewidmet.Die Autoren danken Herrn Prof. Dr. K. Bier für die unterstützung der Arbeit und für wertvolle Diskussionsbeiträge.