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.     

The shear gage: For reliable shear modulus measurements of composite materials

A special strain gage called the shear gage was developed for composite materials testing with notched shear specimens. The shear-gage records the average shear strain across the entire test section between the notches of the losipescu and compact shear specimens rather than just sampling the shear strain over a small region in the center of the test section. Hence, the shear stress/strain response is obtained by dividing the average shear stress (load divided by the cross-sectional area between the notches) by the average shear strain. By placing gages on both faces of the specimen, accurate and repeatable shear-modulus measurements can be made without prior knowledge of the shear strain or stress distributions. This scheme essentially integrates the shear strain through the entire test section. Knowledge of other material properties is not required to accurately determine shear modulus values. The shear gage was tested on a variety of composite and isotropic materials resulting in more reliable shear modulus determination and less scatter than previously possible.     

Fractal porous media III: Transversal Stokes flow through random and Sierpinski carpets

The transversal Stokes flow of a Newtonian fluid through random and Sierpinski carpets is numerically calculated and the transversal permeability derived. In random carpets derived from site percolation, the average macroscopic permeability varies as (- _c)^3/2, close to the critical porosity _c. This exponent is found to be slightly different from the conductivity exponent. Results for Sierpinski carpets are presented up to the fourth generation. The Carman equation is not verified in these two model porous media.     

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.     

Singular Limit and Homogenization for Flame Propagation in Periodic Excitable Media

This paper is concerned with a class of singular equations modelling the combustion of premixed gases in periodic media. The model involves two parameters: the period of the medium |L| and a singular parameter related to the activation energy. The existence of pulsating travelling fronts for fixed and |L| was proved by Berestycki & Hamel in [BH]. In the present paper, we investigate the behaviour of such solutions when More precisely, we establish that pulsating travelling fronts behave like travelling waves, when the period |L| is small and . We also study the convergence of the solution, as goes to zero (and |L| is fixed), toward a solution of a free boundary problem.     

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).     

Numerical validation of the shear compression specimen. Part II: Dynamic large strain testing

Part I of this work addressed quasi-static loading of the shear compression specimen (SCS), which has been especially developed to investigate the shear dominant response of materials at various strain rates. The stress and strain states were characterized numerically. Approximations were presented to reduce the measured load,P, and displacement,d, into equivalent stress and strain . This paper addresses dynamic loading of the SCS. Several simulations were made for representative materials, whose stress-strain behavior is assumed to be rate-independent. The results show that stress wave loading induces strong oscillations in theP-d curve. However, the curve remains smooth in the gage section. The oscillations are about the quasistatic load values, so that with suitable filtering of the dynamicP-d curves, the quasi-static ones are readily recovered. Consequently, the approach that was developed for quasi-static loading of the SCS is now extended to dynamic loading situations. The average strain rate is rather constant and scales linearly with the prescribed velocity. As the plastic modulus becomes smaller, the strain rate reaches higher values. Friction at the end pieces of the specimen is also investigated, and shown to have a small overall influence on the determined mechanical characteristics. This paper thereby confirms the potential of the SCS for large strain testing of materials, using a unified approach, over a large range of strain rates in a seamless fashion.     

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.     

Ü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     

A New Shear-Compression-Specimen for Determining Quasistatic and Dynamic Polymer Properties

A new design of the shear compression specimen (SCS) for investigating the viscoelastic shear response of polymers is presented. The specimen consists of a polymer gage section with two metal ends that remain essentially rigid during deformation. Two closed-form analytic models are developed to predict the average stress and strain in the gage section from the deformation-load histories. This new SCS design and its analytic models are thoroughly evaluated via laboratory measurements and numerical simulations. These simulations show that the deformations in the gage section are more uniform than in the original design, and the distribution of the average shear stress and strain are highly homogenous. The simulation results yield good agreement with those of closed-form analytic results and the experiments demonstrate that the new SCS geometry and its analytic models are as reliable as other commonly employed specimens. It can also generate higher strain rates under usual loading conditions because of its smaller specimen gage length. The need for care in specimen preparation is also discussed in detail as illuminated by the experimental and simulation results.

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.     

Analytical evaluation of stress-strain test data

An analytical model for deducing the actual stress-strain properties from laboratory test results is discussed. As an illustration, an elastic bilinear material is used for unconfined cylindrical compression test conditions, as simulated with a finite element analysis. The results obtained are applicable for assisting in evaluating measured strength and stiffness properties of some clay soils, concrete test cylinders, concrete cores, and rock cores.The quantitative results of this study can be used for interpreting measured stress-strain data for unconfined compression test conditions. The error in measured results is shown to be influenced by Poisson's ratio, length-to-diameter ratio of the specimen, end condition, and ratio of inelastic modulus to initial elastic modulus. Curves for adjusting the measured results to the theoretical results are presented.Nomenclature D specimen diameter - E _i initial elastic stiffness modulus - E _y elastic stiffness modulus beyond the yield stress, plastic or inelastic modulus - L specimen length - axial strain - _av average strain - _g gage length strain - _y yield strain - Poisson's ratio - compressive stress - _av average stress - _t theoretical compressive stress - _y yield stress - _ym measured stress at the yield strain     

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.     

Application of Rayleigh-Ritz method to forced convection turbulent heat transfer

An analytical model to predict heat transfer rates to an incompressible fluid in turbulent flow, with fully developed velocity profile, between a heated plate and a parallel, insulated plate is developed. The model employs van Driest's mixing length expression near the wall, a constant eddy diffusivitiy in the core and a constant turbulent Prandtl number. An approximate solution obtained by employing Rayleigh-Ritz method is shown to compare well with the exact solution obtained by numerical integration of the differential equations. The results are compared with the available experimental data and analytical solutions.

---

Anwendung der Rayleigh-Ritz-Methode auf die Wärmeübertragung bei erzwungener turbulenter Strömung

Zusammenfassung Es wird ein analytisches Modell zur Berechnung der Wärmeübertragung an ein inkompressibles Fluid in turbulenter Strömung mit voll ausgebildetem Geschwindigkeitsprofil zwischen einer beheizten Platte und einer dazu parallelen isolierten Platte angegeben. Das Modell verwendet van Driest's Ausdruck für die wandnahe Mischungslänge, eine konstante Wirbeldiffusivität im Kern und eine konstante turbulente PrandtlZahl. Eine Näherungslösung nach der Rayleigh-Ritz-Methode läßt sich gut mit der exakten Lösung vergleichen, die durch numerische Integration der Differentialgleichungen erhalten wurde. Die Ergebnisse werden mit verfügbaren Versuchswerten und analytischen Lösungen verglichen.

---

Nomenclature A⁺ dimensionless constant in van Driest formula - a⁺ dimensionless distance from the wall after which the eddy diffusivity of momentum is constant - b half-gap of passage - b⁺ dimensionless half-gap=bu^*/ - C_f skin friction coefficient - C_p constant pressure specific heat - d hydraulic mean diameter defined as 4xarea/perimeter=4b - h convective heat transfer coefficient - K⁺ dimensionless constant in van Driest formula - k fluid thermal conductivity - m mass flow rate of fluid - Nu Nusselt number hd/k - P pressure - Pr Prandtl number=/ - Pr_t turbulent Prandtl number=_m/ - q_w heat flux at wall - Re Reynolds number=v_md/ - T Temperature - u⁺ dimensionless velocity=V_x/u^* - u^* friction velocity= - V_x axial velocity - x axial distance from the entrance - x⁺ dimensionless distance=x/d - y distance from the heated wall - y⁺ dimensionless distance=yu^*/ Greek Symbols thermal molecular diffusivity - function equal to (_H+)/ - boundary layer thickness - _H eddy diffusivity of heat - _m eddy diffusivity of momentum - _m0 uniform eddy diffusivity of momentum in the core - dimensionless temperature - T-T_i/q_wd/k uniform heat flux - T-T_w/T_i-T_w uniform temperature - fluid kinematic viscosity - fluid density - fluid shearing stress - bulk mean temperature—fully developed region - fully developed transverse temperature profile Suffixes 1 fully developed - 2 in the entrance region - i at the inlet - m bulk mean value - w at the heated wall     

Near-wall turbulence closure modeling without “damping functions”

An elliptic relaxation model is proposed for the strongly inhomogeneous region near the wall in wall-bounded turbulent shear flow. This model enables the correct kinematic boundary condition to be imposed on the normal component of turbulent intensity. Hence, wall blocking is represented. Means for enforcing the correct boundary conditions on the other components of intensity and on the k — equations are discussed. The present model agrees quite well with direct numerical simulation (DNS) data. The virtue of the present approach is that arbitrary damping functions are not required.     

Turbulent heat and mass transfer for Pr ≫ 1 and law of turbulent diffusion damping at a solid wall

The question of determining the law of damping for the turbulent diffusion coefficient at a smooth wall according to data on mass and heat transfer for Pr 1 is discussed. It is proved that the hypothesis that this law is determined by the first member of the Taylor series expansion of , namely, / = yⁿ ₊ is valid in the Pr range from 10³ to 10⁵ only under the assumption that the subsequent terms in the expansion have smaller coefficients. A statistical analysis of electrochemical and other experiments devoted to this problem shows that apparently n = 3, but singularities in the experimental results do not permit making a final conclusion. Requirements on a conclusive experiment are formulated on the basis of the analysis made.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 172–175, March–April, 1977.     

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.     

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.     

On the shape of a slender axisymmetric nonstationary cavity

Very few studies have been made of three-dimensional nonstationary cavitation flows. In [1, 2], differential equations were obtained for the shape of a nonstationary cavity by means of a method of sources and sinks distributed along the axis of thin axisymmetric body and the cavity. In the integro-differential equation obtained in the present paper, allowance is made for a number of additional terms, and this makes it possible to dispense with the requirement ¦ In ¦ 1 adopted in [1, 2]. The obtained equation is valid under the weaker restriction 1. In [3], the problem of determining the cavity shape is reduced to a system of integral equations. Examples of calculation of the cavity shape in accordance with the non-stationary equations of [1–3] are unknown. In [4], an equation is obtained for the shape of a thin axisymmetric nonstationary cavity on the basis of a semiempirical approach. In the present paper, an integro-differential equation for the shape of a thin axisymmetric nonstationary cavity is obtained to order ² ( is a small constant parameter which has the order of the transverse-to-longitudinal dimension ratio of the system consisting of the cavity-forming body, the cavity, and the closing body). A boundary-value problem is formulated and an analytic solution to the corresponding differential equation is obtained in the first approximation (to terms of order ² In ), A number of concrete examples is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 38–47, July–August, 1980.I thank V. P. Karlikov and Yu. L. Yakimov for interesting discussions of the work.     

Approximate solution of the boundary problem for the equations of a stationary electrostatic charged-particle beam

The solution is given of the equations of a three-dimensional stationary electrostatic beam of charged particles of like sign filling the region between two nearby curvilinear surfaces. We assume that the flow is nonrotational and nonrelativistic and that the velocity vector is a single-valued function. The solution is constructed in the form of an asymptotic series in powers of the small parameter , which is the ratio of the characteristic transverse (a) and longitudinal (l) dimensions of the problem. The first dimension is taken to be the distance between the electrodes, andl defines the scale at which the geometric and physical parameters (emitter curvature, electric field E on the emitter, and the emission current density J) change noticeably. The emission regimes limited by the space charge (-regime), temperature (T-regime), and the case of nonzero initial velocity (U-regime) are studied. The asymptotic behavior is given by the formulas for the corresponding one-dimensional flow between parallel surface.The solution of the boundary problem for emission in the-regime reduces to determination of the emission current density J for fixed electrode geometry and given accelerating voltage. The corresponding formulas are presented, retaining terms of order ³.Two approximations with respect to are performed for the T- and U-regimes. Here the unknown quantity for given properties of the emitting surface (J) will be the electric field E.The results provided by the constructed expansions are compared with the exact solution for flow from a planar emitter along circular trajectories [1]. As an example we examine the two-dimensional problem of flow between two nearby circular cylindrical electrodes with disruption of the coaxiality.The conventional tensor notations are used.