Combined effect of external turbulence and strong acceleration on boundary layer transition in a nozzle

The effect of external turbulence on the boundary layer flow in a convergent-divergent nozzle with a high expansion ratio has been studied numerically. The external turbulence was simulated by the turbulent viscosity _e, for which we used the partial differential equation that serves to close the system of boundary layer equations [1–4]. It was found that there exists a critical value _cr such that for all _e< _cr the flow regime in the nozzle remains perfectly laminar, whereas for _ecr a laminar-turbulent transition takes place and the boundary layer in the supersonic part of the nozzle becomes turbulent. For postcritical values of _e the heat fluxes and friction losses are approximately an order greater than for precritical values. With increase in the Reynolds number, determined from the parameters in the nozzle throat, the value of _cr decreases; as the coordinate of the onset of boundary layer formation is displaced in the direction of flow the value of _cr increases.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 34–37, January–February, 1906.The authors are grateful to L. V. Gogish for participating in the discussion of the results.     

Effect of free-stream turbulence on surface friction in a turbulent boundary layer

The combined effect of the turbulence intensity , the turbulence scaleL, and the Reynolds number Re^** on the surface friction coefficientc _f in a turbulent boundary layer is studied. The dependence of the relative friction increment on the equivalent turbulence level _cq, which takes into account the simultaneous variation in ,L and Re^**, is determined. The threshold value _cq ^* below which the value ofc _f does not depend on _cq is found.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 65–75, March–April, 1995.     

Positively invariant regions for a problem in phase transitions

Positively invariant regions for the system v _t + p(W) _x = V _xx , W _t –V _x = W _xx are constructed where p < 0, w < , w > , p(w) = 0, w , > 0. Such a choice of p is motivated by the Maxwell construction for a van der Waals fluid. The method of an analysis is a modification of earlier ideas of Chueh, Conley, & Smoller [1]. The results given here provide independent L bounds on the solution (w, v).Dedicated to Professor James Serrin on the occasion of his sixtieth birthday     

Generation of vortices in a partially filled,rapidly rotating cylinder

We describe a system in which vortices are shed from a cylindrical free surface approximately centered in a rotating flow. Shedding is controlled by the parameter =2 g/ ² d, where g, , d denote gravity, rotation rate and the diameter of the free surface. We find vortex shedding for >0.162 and no vortex shedding for < 0.0847. The range depends on the aspect ratio L/d, where L is the column length, in a nonmonotonic fashion. These results are independent of viscosity and surface tension for small values of these parameters.Now at Martin Marietta, Orlando Aerospace, PO Box 5837, Mail Point 150, Orlando, FL 32855, USA     

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

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

Thin-Walled Beams: the Case of the Rectangular Cross-Section

In this paper we present an asymptotic analysis of the three-dimensional problem for a thin linearly elastic cantilever =×(0,l) with rectangular cross-section of sides and ², as goes to zero. Under suitable assumptions on the given loads, we show that the three-dimensional problem converges in a variational sense to the classical one-dimensional model for extension, flexure and torsion of thin-walled beams. Mathematics Subject Classifications (2000) 474K20, 74B10, 49J45.     

An absence of a certain class of periodic solutions in the Navier-Stokes equations

In the case of the 3D Navier-Stokes equations, it is proved that there exists a constant>0 with the following property: Every time-periodic solution with a period smaller than is necessarily a stationary solution. An explicit formula for is also provided.     

Lower semicontinuity of the attractor for a singularly perturbed hyperbolic equation

For a smooth, bounded domain R, n 3, and a real, positive parameter, we consider the hyperbolic equationu _tt +u _t –u=–f(u) –g in with Dirichlet boundary conditions. Under certain conditions onf, this equation has a global attractorA inH ₀ ¹ () ×L ²(). For=0, the parabolic equation also has a global attractor which can be naturally embedded into a compact setA ₀ inH ₀ ¹ () ×L ²(). If all of the equilibrium points of the parabolic equation are hyperbolic, it is shown that the setsA are lower semicontinuous at=0. Moreover, we give an estimate of the symmetric distance betweenA ₀ andA .     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

Application of dispersion theory to time domain reflectometry in soils

With time domain reflectometry (TDR) two dispersive parameters, the dielectric constant, , and the electrical conductivity, can be measured. Both parameters are nonlinear functions of the volume fractions in soil. Because the volume function of water ( _w) can change widely in the same soil, empirical equations have been derived to describe these relations. In this paper, a theoretical model is proposed based upon the theory of dispersive behaviour. This is compared with the empirical equations. The agreement between the empirical and theoretical aproaches was highly significant: the ( _w) relation of Topp et al. had a coefficient of determination r ² = 0.996 and the (_u) relation of Smith and Tice, for the unfrozen water content, _u, at temperatures below 0°C, had an r ² = 0.997. To obtain ( _w) relations, calibration measurements were performed on two soils: Caledon sand and Guelph silt loam. For both soils, an r ² = 0.983 was obtained between the theoretical model and the measured values. The correct relations are especially important at low water contents, where the interaction between water molecules and soil particles is strong.     

Discrete time semigroup transformations with random perturbations

Let (X, ) and (Y,C) be two measurable spaces withX being a linear space. A system is determined by two functionsf(X): X X and:X×YX, a (small) positive parameter and a homogeneous Markov chain {y _n } in (Y,C) which describes random perturbations. States of the system, say {x _n X, n=0, 1,}, are determined by the iteration relations:x _n+1 =f(x _n )+(x _n ,Y_n+1) forn0, wherex ₀ =x ₀ is given. Here we study the asymptotic behavior of the solutionx _n as 0 andn under various assumptions on the data. General results are applied to some problems in epidemics, genetics and demographics.Supported in part by NSF Grant DMS92-06677.Supported in part by NSF Grant DMS93-12255.     

Über den Einfluß von Querkrümmung und Dichte bei inhomogenen turbulenten Freistrahlen

Zusammenfassung Zur Berechnung turbulenter Strömungen wird das k--Modell im Ansatz für die turbulente Scheinzähigkeit erweitert, so daß es den Querkrümmungs- und Dichteeinfluß auf den turbulenten Transportaustausch erfaßt. Die dabei zu bestimmenden Konstanten werden derart festgelegt, daß die bestmögliche Übereinstimmung zwischen Berechnung und Messung erzielt wird. Die numerische Integration der Grenzschichtgleichungen erfolgt unter Verwendung einer Transformation mit dem Differenzenverfahren vom Hermiteschen Typ. Das erweiterte Modell wird auf rotationssymmetrische Freistrahlen veränderlicher Dichte angewendet und zeigt Übereinstimmung zwischen Rechnung und Experiment.

---

On the influence of transvers-curvature and density in inhomogeneous turbulent free jets

The prediction of turbulent flows based on the k- model is extended to include the influence of transverse-curvature and density on the turbulent transport mechanisms. The empirical constants involved are adjusted such that the best agreement between predictions and experimental results is obtained. Using a transformation the boundary layer equations are solved numerically by means of a finite difference method of Hermitian type. The extended model is applied to predict the axisymmetric jet with variable density. The results of the calculations are in agreement with measurements.

---

Bezeichnungen Wirbelabsorptionskoeffizient - c_i Massenkonzentration der Komponente i - c_D, c_L, c, c₁, c₂ Konstanten des Turbulenzmodells - d Düsendurchmesser - E bezogene Dissipationsrate - f bezogene Stromfunktion - f Korrekturfunktion für die turbulente Scheinzähigkeit - j turbulenter Diffusionsstrom - k Turbulenzenergie - k_i Schrittweite in -Richtung - K dimensionslose Turbulenzenergie - L turbulentes Längenmaß - M_i Molmasse der Komponente i - p Druck - allgemeine Gaskonstante - r Querkoordinate - r_0,5 Halbwertsbreite der Geschwindigkeit - r_0,5c Halbwertsbreite der Konzentration - T Temperatur - u Geschwindigkeitskomponente in x-Richtung - v Geschwindigkeitskomponente in r-Richtung - x Längskoordinate - y allgemeine Funktion - Y_i diskreter Wert der Funktion y - Relaxationsfaktor für Iteration - turbulente Dissipationsrate - transformierte r-Koordinate - kinematische Zähigkeit - Exponent - transformierte x-Koordinate - Dichte - _k, Konstanten des Turbulenzmodells - Schubspannung - allgemeine Variable - Stromfunktion - Turbulente Transportgröße Indizes 0 Strahlanfang - m auf der Achse - r mit Berücksichtigung der Krümmung - t turbulent - mit Berücksichtigung der Dichte - im Unendlichen - Schwankungswert oder Ableitung einer Funktion - – Mittelwert Herrn Professor Dr.-Ing. R. Günther zum 70. Geburtstag gewidmet     

Direct measurement of the density field using high speed differential interferometry

The optical method presented here is based on Wollaston prism differential interferometry using a white light source. It yields high-speed, instantaneous interferograms which directly display the isochoric lines of an unsteady, two-dimensional flow. The setup used is equivalent to conventional reference beam interferometric systems, but remains differential since one of the beams — the reference one — lies inside the undisturbed upstream flow. The method was used to obtain an instantaneous, high speed display of the unsteady flow around a cylinder spanning the test section.List of symbols d Cylinder diameter - t Time - t Interval of time between two consecutive interferograms - M Mach number - Density - P Pressure - P Pressure fluctuation - Angle of the prisms making up the birefringent system - Birefringent angle = () - gl Wavelength - n ₀ Ordinary index of the crystal - n _e Extraordinary index of the crystal - D Diameter of the spherical mirror - R Radius of the curvature of the spherical mirror - i Fringe width - 2g prism thickness - X Distance between the two beams interfering at the level of the spherical mirror The quantities relating to upstream conditions are indexed      

On a generalized nonlinearK-ɛ model for turbulence that models relaxation effects

In this paper, based on a similarity that exists between the constitutive relations for turbulent mean flow of a Newtonian fluid and that for the laminar flow of a non-Newtonian fluid, and making use of extended thermodynamics, we develop a generalized nonlinearK- model, whose approximate form includes the standardK- model and the nonlinearK- model of Speziale (1987) as special cases. Our nonlinearK- model, which is frame indifferent, can predict relaxation of the Reynolds stress, unlike most standardK- models. Also, our model is in keeping with that of Yakhotet al. (1992). Most interestingly, the linearized form of our model bears a striking resemblance to the model due to Yoshizawa and Nisizima (1993); however, it has been obtained from a totally different perspective.     

Extensional viscosity from entrance pressure drop measurements

Extensional rheological properties are important in characterization and processing of polymeric liquids. The use of entrance pressure drop to obtain extensional viscosity is particularly attractive because it can be applied to both low and high viscosity liquids using the Bagley correction obtained from a conventional capillary rheometer.Low density polyethylene of three different melt index values, including IUPAC-X (a different batch of IUPAC-A), and a high density polyethylene were tested using a commercial capillary rheometer. The entrance pressure drop (P _en ) was obtained with a zero-length orifice die with an abrupt contraction. The contraction ratio was 12:1. Predictions from several approximate analyses to calculate the uniaxial extensional viscosity _u (using an axisymmetric contraction) from P _en were compared. These comparisons are summarized in the appendices.Due to the transient nature of contraction flows, _u is also a function of the strain (). This was examined by comparing _u from P _en (Cogswell's analysis was chosen for convenience) with transient extensional viscosity ( _u ⁺) at different magnitudes of from fiber-windup technique (Padmanabhan et al., 1996). _u ⁺ at 3 was found to be close to _u from P _en (using Cogswell's analysis) for two LDPE samples that had fiber-windup data available. The magnitude of the strain in the contraction did not vary with strain rate.Dedicated to the memory of Tasos Papanastasiou     

The problems of nonlinear bending for orthotropic rectangular plate with four clamped edges

（黄家寅）（秦圣立）ＴＨＥＰＲＯＢＬＥＭＳＯＦＮＯＮＬＩＮＥＡＲＢＥＮＤＩＮＧＦＯＲＯＲＴＨＯＴＲＯＰＩＣＲＥＣＴＡＮＧＵＬＡＲＰＬＡＴＥＷＩＴＨＦＯＵＲＣＬＡＭＰＥＤＥＤＧＥＳ￥ＨｕａｎｇＪｉａｙｉｎ;ＱｉｎＳｈｅｎｇｌｉ（ＱｕｆｕＮｏｒｍａｌＵｎ...     

Theory of a laval nozzle for a two-phase mixture containing particles of small lag

A two-velocity and two-temperature model is considered for a continuous medium in relation to the flow of a mixture of gas and particles in the subsonic, transsonic, and supersonic parts of a Laval nozzle. It is assumed that the particles are small, and hence that the coefficients f and ^q, which define the interaction with the gas, are large (these coefficients are inversely proportional to the square of the particle radius for a Stokes mode of flow). This means that the velocity or thermal lag of the particles relative to the gas is small. The solution is sought as expansions with respect to the small parameters ₁=1/f and ₂=1/^q.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 89–100, July–August, 1973.     

Dirichlet problem of higher order elliptic equation with perturbation both in differential operator and in boundary

This paper is the continuation of article [7]. It gives further results about the asymptotic expression for the solution of higher order elliptic equation in the case of boundary perturbation combined with operator perturbation. When unperturbed problemA ₀ is not on the spectrum, the asymptotic expression for the solution of perturbation problemA may be expanded with respect to the small parameter . WhileA ₀ is on the spectrum, the asymptotic expression of the solution contains negative powers of the small parameter . The approximation of arbitrary order to the solution is considered and the recursive formula for the general term and the estimation of remainder term are given.     

Limit cycle analysis of a class of strongly nonlinear oscillation equations

The limit cycle of a class of strongly nonlinear oscillation equations of the form % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiqadwhagaWaaiabgUcaRmXvP5wqonvsaeHbbjxAHXgiofMCY92D% aGqbciab-DgaNjab-HcaOiaadwhacqWFPaqkcqWF9aqpcqaH1oqzca% WGMbGaaiikaiaadwhacaGGSaGabmyDayaacaGaaiykaaaa!50B8!\[\ddot u + g(u) = \varepsilon f(u,\dot u)\] is investigated by means of a modified version of the KBM method, where is a positive small parameter. The advantage of our method is its straightforwardness and effectiveness, which is suitable for the above equation, where g(u) need not be restricted to an odd function of u, provided that the reduced equation, corresponding to =0, has a periodic solution. A specific example is presented to demonstrate the validity and accuracy of our 09 method by comparing our results with numerical ones, which are in good agreement with each other even for relatively large .     

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.