Asymptotic analysis of equilibrium states for rotating turbulent flows

The equilibrium states of homogeneous turbulence simultaneously subjected to a mean velocity gradient and a rotation are examined by using asymptotic analysis. The present work is concerned with the asymptotic behavior of quantities such as the turbulent kinetic energy and its dissipation rate associated with the fixed point (/kS)=0, whereS is the shear rate. The classical form of the model transport equation for (Hanjalic and Launder, 1972) is used. The present analysis shows that, asymptotically, the turbulent kinetic energy (a) undergoes a power-law decay with time for (P/)<1, (b) is independent of time for (P/)=1, (c) undergoes a power-law growth with time for 1<(P/)<(C ₂–1), and (d) is represented by an exponential law versus time for (P/)=(C ₂–1)/(C ₁–1) and (/kS)>0 whereP is the production rate. For the commonly used second-order models the equilibrium solutions forP/,II, andIII (whereII andIII are respectively the second and third invariants of the anisotropy tensor) depend on the rotation number when (P/kS)=(/kS)=0. The variation of (P/kS) andII versusR given by the second-order model of Yakhot and Orzag are compared with results of Rapid Distortion Theory corrected for decay (Townsend, 1970).     

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.     

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.     

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.     

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.     

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

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

Ü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     

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     

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.     

Zur Variationsformulierung nichtlinearer Randwertprobleme

Übersicht Es werden verschiedene Bedingungen aufgestellt, die es erlauben, die durch die beiden (Systeme von) nichtlinearen DifferentialgleichungenA (u, ) = q, B (u, ) = und Randbedingungen zusammen mit den nichtlinearen algebraischen Relationenq = C(u, ), = D(u, ) beschriebene Aufgabe durch äquivalente Variationsprobleme zu ersetzen. Dabei zeigt sich ein enger Zusammenhang mit den in der Festkörpermechanik wohlbekannten Prinzipien der virtuellen Verschiebungen und der virtuellen Kräfte. Die auf systematischem Weg konstruierten Variationsfunktionale enthalten viele in der Physik bekannte Funktionale als Sonderfälle, insbesondere jene, die in der Elastomechanik nach Green, Castigliano, Hellinger, Reißner, Hu und Washizu benannt werden.

---

Summary In this paper there are established various conditions which allow a variational formulation of the problem described by the two (systems of) nonlinear differential equationsA(u, ) = q, B(u, ) = and boundary conditions together with the nonlinear algebraic relationsq = C(u, ), = D(u, ). Besides a close relationship is revealed to the principles of virtual displacements and virtual forces which are wellknown in solid mechanics. The systematically constructed variational functional contain many functionals in physics as special cases, mainly those of Green, Castigliano, Hellinger, Reißner, Hu and Washizu in elastomechanics.

     

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.     

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     

Rheology of concentrated coagulating suspensions in non-aqueous media

The rheology of concentrated coagulating suspensions is analysed on the basis of the following model: (i) at low shear rates, the shear is not distributed homogeneously but limited to certain shear planes; (ii) the energy dissipation during steady flow is due primarily to the overcoming of viscous drag by the suspended particles during motion caused by encounters of particles in the shear planes. This model is called the giant floc model.With increasing shear rate the distance between successive shear planes diminishes, approaching the suspended particles' diameter at average shear stresses of 88–117 Pa in suspensions of 78 µm particles (glass ballotini coated by a hydrophobic layer) in glycerol — water mixtures, at solid volume fractions between 0.35 and 0.40. Smaller particles form a more persistent coagulation structure. The average force necessary to separate two touching 78 µm particles is too large to be accounted for by London-van der Waels forces; thus coagulation is attributed to bridging connections between polymer chains protruding from the hydrophobic coatings.The frictional ratio of the glass particles in these suspensions is of the order of 10. Coagulation leads to build-up of larger structural units at lower shear rates; on doubling the shear rate the average distance between the shear planes decreases by a factor of 0.81 to 0.88. A inter-shear plane distance - A Hamaker constant - b radius of primary particles - f frictional ratio - F _A attractive force between two particles - g acceleration due to gravity - H distance between the surfaces of two particles - K proportionality constant in power law - l fraction of distance by which a moving particle entrains its neighbours - l effective length of inner cylinder in the rheometer - M torque experienced by inner cylinder during measurements - n exponent in power law - n ₀ ,n ₁ ,n ₂ constants in extended power law - NC _hex number of contacts, per mm², between particles in adjacent layers with an average degree of occupation, assuming a hexagonal arrangement of the particles within the layers - NC _cub asNC _hex, but with a cubical arrangement - p () d increase of slippage probability when the shear stress increases from to + d - q average coordination number of a particle in a coagulate - R _i radius of inner cylinder of rheometer - R _u radius of outer cylinder of rheometer - t _i time during which particlei moves - t ₀ time during which a particle bordering a shear plane moves from its rectilinear course, on meeting another particle - u angle between the direction of motion, and the line connecting the centers of two successive particles bordering a shear plane - V _A attractive energy between two particles - x, y, z Cartesian coordinates:x — the direction of motion;y — the direction of the velocity gradient - y ₀ ,z ₀ y, z value of a particle meeting another particle, when both are far removed from each other - y ₀ spread iny ₀ values - —2/n - ₀ capture efficiency - shear rate - average shear rate calculated for a Newtonian liquid - _i distance by which particlei moves - ₀ distance by which a particle bordering a shear plane moves from its rectilinear course, when it encounters another particle - square root of area occupied by a particle bordering a shear plane, in this plane - _c energy dissipated during one encounter of two particles bordering a shear plane - _p energy dissipated by one particle - energy dissipated per unit of volume and time during steady flow - viscosity - _app calculated as if the liquid is Newtonian - ₀ viscosity of suspension medium - _PL lim - [] intrinsic viscosity - _diff - _{diff, rel diff}/ ₀ - standard deviation of distribution ofy ₀ values - shear stress - _n average shear stress at the highest values applied - mass average particle diameter - _n number average particle diameter - solid volume fraction - _eff effective solid volume fraction in Dougherty-Krieger relation - _max maximum solid volume fraction permitting flow - _i angular velocity of inner cylinder in rheometer during measurements     

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      

A class of nonlinear,completely integrable abstract wave equations

IfL is a positive self-adjoint operator on a Hubert spaceH, with compact inverse, the second-order evolution equation int,u+Lu+u _H ² u=0 has an infinite number of first integrals, pairwise in involution. It follows from this that no nontrivial solution tends weakly to 0 inH ast. Under an additional separation assumption on the eigenvalues ofL, all trajectories (u,u) are relatively compact inD(L ^1/2)×H. Finally, if all the eigenvalues are simple, the set of initial values of quasi-periodic solutions is dense in the ball B=(u ₀,u ₀ )D(L ^1/2)×H; L^1/2 u _{0 H} ² +u ² < for sufficiently small.     

Invariant Manifolds for Nonautonomous Systems with Application to One-Step Methods

In this paper we study the existence of invariant manifolds for a special type of nonautonomous systems which arise in the study of discretization methods. According to [10], a one-step scheme of step-size for an autonomous system can be interpreted as the -flow of a perturbed nonautonomous system. The perturbation is rapidly forced in the sense that it is periodic with respect to time with period . Assuming a saddle node for the autonomous system, we prove that these rapidly forced perturbations have center manifolds which exist in a uniform neighborhood and which converge to a center manifold of the autonomous system as tends to zero. Our results are applied to obtain a smooth continuation as well as estimates of the well known center manifolds for one-step schemes. They also form the basis for studying saddle-node homoclinic orbits under discretization.     

Numerical analysis of the transient conjugated heat transfer in a circular duct with a power-law fluid

We treat numerically in this paper, the transient analysis of a conjugated heat transfer process in the thermal entrance region of a circular tube with a fully developed laminar power-law fluid flow. We apply the quasi-steady approximation for the power-law fluid, identifying the suitable time scales of the process. Thus, the energy equation in the fluids is solved analytically using the well-known integral boundary layer technique. This solution is coupled to the transient energy equation for the solid where the transverse and longitudinal heat conduction effects are taken into account. The numerical results for the temporal evolution of the average temperature of the tube wall, _av, is plotted for different nondimensional parameters such as conduction parameter, , the aspect ratios of the tube, and ₀ and the index of power-law fluid, n.     

Interpolation formulas for maxwell viscosity of certain metals as a function of shear-strain intensity and temperature

The purpose of this study is the construction of interpolation formulas for the dependence of Maxwell viscosity, a quantity which is the reciprocal of shear-strain relaxation time , on shear-strain intensity and temperature for several metals: iron, aluminum, copper, and lead. This function was interpolated in various temperature and deformation velocity ranges in accordance with available experimental data for iron (0 10⁷ sec^–1, 200 ° T 1500 °); aluminum (0 10⁷ sec^–1, 300 ° T 900 °); copper (0 10⁵ sec^–1, 300 ° T 1300 °); lead (0 10⁶ sec^–1, 90 ° T 400 °); temperatures in °K.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 114–118, July–August, 1974.