On the thermodynamics of elastic-plastic materials with temperature-dependent moduli and yield stresses

Summary The subject of this article is the thermodynamics of perfect elastic-plastic materials undergoing unidimensional, but not necessarily isothermal, deformations. The first and second laws of thermodynamics are employed in a form in which only the following quantities appear: the temperature , the elastic strain _e, the plastic strain _p, the elastic modulus (gq), the yield strain (gq), the heat capacity (_e, _p,), the latent elastic heat _e(_e, _p, ), and the latent plastic heat _p(_e, _p, ). Relations among the response functions , , , _e, and _p are derived, and it is shown that a set of these relations gives a necessary and sufficient condition for compliance with the laws of thermodynamics. Some observations are made about the existence and uniqueness of energy and entropy as functions of state.Dedicated to Clifford Truesdell on the occasion of his 60th birthdayThis research was supported by the U.S. National Science Foundation.     

Parameter identification in a soil with constant diffusivity

We consider infiltration into a soil that is assumed to have hydraulic conductivity of the form K = K = K_se^h and water content of the form = K – _r. Here h denotes capillary pressure head while K_s, , and _r represent soil specific parameters. These assumptions linearize the flow equation and permit a closed form solution that displays the roles of all the parameters appearing in the hydraulic function K and . We assume K_s and _r to be known. A measurement of diffusivity fixes the product of and resulting in a parameter identification problem for one parameter. We show that this parameter identification problem, in some cases, has a unique solution. We also show that, in some cases, this parameter identification problem can have multiple solutions, or no solution. In addition it is shown that solutions to the parameter identification problem can be very sensitive to small changes in the problem data.     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

Determination of water retention in stratified porous materials

Predicted and measured water-retention values,(), were compared for repacked, stratified core samples consisting of either a sand with a stone-bearing layer or a sand with a clay loam layer in various spatial orientations. Stratified core samples were packed in submersible pressure outflow cells, then water-retention measurements were performed between matric potentials,, of 0 to -100 kPa. Predictions of() were based on a simple volume-averaging model using estimates of the relative fraction and() values of each textural component within a stratified sample. In general, predicted() curves resembled measured curves well, except at higher saturations in a sample consisting of a clay loam layer over a sand layer. In this case, the model averaged the air-entry of both materials, while the air-entry of the sample was controlled by the clay loam in contact with the cell's air-pressure inlet. In situ, avenues for air-entry generally exist around clay layers, so that the model should adequately predict air-entry for stratified formations regardless of spatial orientation of fine versus coarse layers. Agreement between measured and predicted volumetric water contents,, was variable though encouraging, with mean differences between measured and predicted values in the range of 10%. Differences in of this magnitude are expected due to variability in pore structure between samples, and do not indicate inherent problems with the volume averaging model. This suggets that explicit modeling of stratified formations through detailed characterization of the stratigraphy has the potential of yielding accurate() values. However, hydraulic-equilibration times were distinctly different for each variation in spatial orientation of textural layering, indicating that transient behavior during drainage in stratified formations is highly sensitive to the stratigraphic sequence of textural components, as well as the volume fraction of each textural component in a formation. This indicates that prolonged residence times of water, nutrients, and pollutants are likely within finer-textured layers, when conditions have resulted in drainage of underlying coarser-textured strata.     

Forced oscillations of a rotating shaft with nonlinear spring characteristics and internal damping (1/2 order subharmonic oscillations and entrainment)

Nonlinear forced oscillations of a rotating shaft with nonlinear spring characteristics and internal damping are studied. In particular, entrainment phenomena at the critical speeds of 1/2 order subharmonic oscillations of forward and backward whirling modes are investigated. A self-excited oscillation appears in the wide range above the major critical speed. The amplitude of this oscillation reaches a limit value and then a self-sustained oscillation occurs. In the vicinity of a 1/2 order subharmonic oscillation of a forward whirling mode, a self-excited oscillation is entrained by a subharmonic oscillation. In the vicinity of a 1/2 order subharmonic oscillation of a backward whirling mode, either a self-excited oscillation or a subharmonic oscillation occurs.Experiments were made by an elastic rotating shaft with a disc. Nonlinearity in its restoring force was due to an angular clearance of a bearing and internal damping was due to friction between the shaft and an inner ring of the bearing. A self-excited oscillation was observed in the range above the major critical speed and this self-excited oscillation was entrained by a 1/2 order subharmonic oscillation of a forward whirling mode.Nomenclature O–xyz rectangular coordinate system - , _x, _y inclination angle of a shaft and its projections on the xz- and yz-planes - _x, _y inclination angles in rotating coordinates - , polar coordinates - I _p polar moment of inertia of a rotor - I diametral moment of inertia of a rotor - i _p ratio of I _p to I - dynamic unbalance of a rotor - rotating speed (angular velocity) - F magnitude of a dynamic unbalance force, F = (1 – i _p )² - c external damping coefficient - h internal damping coefficient - t time - D _x , D _y internal damping terms in stationary coordinates - D _x , D _y internal damping terms in rotating coordinates - N _x , N _y nonlinear terms in restoring forces     

Stress analysis of plates with a circular hole reinforced by flange reinforcing member

This paper deals with the problem of stress analysis of plates with a circular hole reinforced by flange reinforcing member. The so called flange reinforcing member here means that the reinforcing member is built up by setting shapes or bars with any section shape on both sides of the plates along the edge of the hole. Two cases of external loads are considered. In one case the external loads are stressesσ_X^(∞),σ_Y^(∞),and τ_XY^(∞) acting at infinite point of the plate, and in the other the external loads are linear distributed normal stresses. The procedure of solving the problems mentioned above consists of three steps. Firstly, the reinforcing member is taken out from the plates and considered to be a circular bar being solved to determine its deformation under the action of radial force q₀(θ) and tangential force t₀(θ) which are forces acting upon each other between reinforcing member and plate. Secondly, the displacements of plate with a circular hole under the action of q₀(θ) and t₀(θ) and external loads are determined. Finally, forces q₀(θ) and t₀(θ) are obtained by the compatibility of deformations between reinforcing member and plate. Then the internal forces and displacements of reinforcing member and plate are deduced from q₀(θ) and t₀(θ) obtained.     

Nonstationary vibration of a rotating shaft with nonlinear spring characteristics during acceleration through a critical speed (A critical speed of a summed-and-differential harmonic oscillation)

Nonstationary vibration of a flexible rotating shaft with nonlinear spring characteristics during acceleration through a critical speed of a summed-and-differential harmonic oscillation was investigated. In numerical simulations, we investigated the influence of the angular acceleration , the initial angular position of the unbalance _n and the initial rotating speed on the maximum amplitude. We also performed experiments with various angular accelerations. The following results were obtained: (1) the maximum amplitude depends not only on but also on _n and : (2) when the initial angular position _n changes. the maximum amplitude varies between two values. The upper and lower bounds of the maximum amplitude do not change monotonously for the angular acceleration: (3) In order to always pass the critical speed with finite amplitude during acceleration. the value of must exceed a certain critical value.Nomenclature O-xyz rectangular coordinate system - , ₁, ₁ inclination angle of rotor and its projections to thexy- andyz-planes - I _r polar moment of inertia of rotor - I diametral moment of inertia of rotor - i _r ratio ofI _r toI - dynamic unbalance of rotor - directional angle of fromx-axis - c damping coefficient - spring constant of shaft - N _nt ,N _nt nonlinear terms in restoring forees in ₁ and ₁ directions - ₄ representative angle - a small quantity - V. V _u .V _N potential energy and its components corresponding to linear and nonlinear terms in the restoring forees - directional angle - _n coefficients of asymmetrical nonlinear terms - _n coefficients of symmetrical nonlinear terms - coefficients of asymmetrical nonlinear terms experessed in polar coordinates - coefficients of symmetrical nonlinear terms expressed in polar coordinates - rotating speed of shaft - t time - _n initial angular position of att=0 - p natural frequency - p ₁.p _t natural frequencies of forward and backward precessions - , ₁, ₁ total phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - , ₁, ₁ phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - P, R _t ,R _b amplitudes of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - difference between phases ( = _f – _u) - acceleration of rotor - initial rotating speed - t _t ,r _b amplitudes of nonstationary oscillation during acceleration - (r _t )_max, (r _b )_max maximum amplitudes of nonstationary oscillation during acceleration - (r ₁ ¹ )_max, (r _b ¹ )_max maximum value of angular acceleration of non-passable case - ₀ critical value over which the rotor can always pass the critical speed - p ₁,p ₂,p ₃,p ₄ natural frequencies of experimental apparatus     

Flow of a simple non-newtonian fluid past a sphere

Summary Creeping flow past a sphere is solved for a limiting case of fluid behaviour: an abrupt change in viscosity.List of Symbols d _ij Component of rate-of-deformation tensor - F _d Drag force exerted on sphere by fluid - G ^(d) Coefficients in expression for _ij in terms of d _ij - G _YOJK ^(d) Coefficients in power series representing G ^(d) - R Radius of sphere - r Spherical coordinate - V Velocity of fluid very far from sphere - v _i Component of the velocity vector - x Dimensionless radial distance, r/R - x _i Rectangular Cartesian coordinate - Dimensionless quantity defined by (26) - ^(d) Potential defined by (7) - Value of x denoting border between Regions 1 and 2 as a function of - ₁, ₂ Lower and upper limiting viscosities defined by (10) - Spherical coordinate - * Value of for which =1 - Value of denoting border between regions 1 and 2 as a function of x - Newtonian viscosity - _ij Component of the stress tensor - Spherical coordinate - ₁, ₂ Stream functions defined by (12) and (14) - Second and third invariants of the stress tensor and of the rate-of-deformation tensor, defined by (3)     

The secondary flow induced around a sphere in an oscillating stream of elastico-viscous liquid

The present paper is devoted to the theoretical study of the secondary flow induced around a sphere in an oscillating stream of an elastico-viscous liquid. The boundary layer equations are derived following Wang's method and solved by the method of successive approximations. The effect of elasticity of the liquid is to produce a reverse flow in the region close to the surface of the sphere and to shift the entire flow pattern towards the main flow. The resistance on the surface of the sphere and the steady secondary inflow increase with the elasticity of the liquid.Nomenclature a radius of the sphere - b ^ik contravariant components of a tensor - e contravariant components of the rate of strain tensor - F() see (47) - G total nondimensional resistance on the surface of the sphere - g _ik covariant components of the metric tensor - f, g, h secondary flow components introduced in (34) - k ₀ measure of relaxation time minus retardation time (elastico-viscous parameter) - K =k ₀ ²/V ₀ ² , nondimensional parameter characterizing the elasticity of the liquid - n measure of the ratio of the boundary layer thickness and the oscillation amplitude - N, T defined in (44) - p arbitrary isotropic pressure - p _ik covariant components of the stress tensor - p ^ik contravariant components of the stress tensor associated with the change of shape of the material - R =V ₀ a/v, the Reynolds number - S =a/V ₀, the Strouhall number - r, , spherical polar coordinates - u, v, w r, , component of velocity - t time - V(, t) potential velocity distribution around the sphere - V ₀ characteristic velocity - u, v, t, y, P nondimensional quantities defined in (15) - reciprocal of s - density - defined in (32) - defined in (42) - ₀ limiting viscosity for very small changes in deformation velocity - complex conjugate of - oscillation frequency - = ₀/, the kinematic coefficient of viscosity - , defined in (52) - (, y) stream function defined in (45) - =(NT/2n)^1/2 y - /t convective time derivative (1) _ik      

Analysis of turbulent boundary layer interaction with an external supersonic flow behind a step

An integral method of analyzing turbulent flow behind plane and axisymmetric steps is proposed, which will permit calculation of the pressure distribution, the displacement thickness, the momentum-loss thickness, and the friction in the zone of boundary layer interaction with an external ideal flow. The characteristics of an incompressible turbulent equilibrium boundary layer are used to analyze the flow behind the step, and the parameters of the compressible boundary layer flow are connected with the parameters of the incompressible boundary layer flow by using the Cowles-Crocco transformation.A large number of theoretical and experimental papers devoted to this topic can be mentioned. Let us consider just two [1, 2], which are similar to the method proposed herein, wherein the parameter distribution of the flow of a plane nearby turbulent wake is analyzed. The flow behind the body in these papers is separated into a zone of isobaric flow and a zone of boundary layer interaction with an external ideal flow. The jet boundary layer in the interaction zone is analyzed by the method of integral relations.The flow behind plane and axisymmetric steps is analyzed on the basis of a scheme of boundary layer interaction with an external ideal supersonic stream. The results of the analysis by the method proposed are compared with known experimental data.Notation x, y longitudinal and transverse coordinates - X, Y transformed longitudinal and transverse coordinates - , ^*, ^** boundary layer thickness, displacement thickness, momentum-loss thickness of a boundary layer - , ^*, ^** layer thickness, displacement thickness, momentum-loss thickness of an incompressible boundary layer - u, velocity and density of a compressible boundary layer - U, velocity and density of the incompressible boundary layer - , stream function of the compressible and incompressible boundary layers - , dynamic coefficient of viscosity of the compressible and incompressible boundary layers - r₁ radius of the base part of an axisymmetric body - r radius - R transformed radius - M Mach number - friction stress - p pressure - a speed of sound - s enthalpy - v Prandtl-Mayer angle - P Prandtl number - P_t turbulent Prandtl number - r₂ radius of the base sting - b step depth - =0 for plane flow - =1 for axisymmetric flow Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 33–40, May–June, 1971.In conclusion, the authors are grateful to M. Ya. Yudelovich and E. N. Bondarev for useful comments and discussions.     

The three dimensional incompressible laminar boundary layer on a spinning cone

Summary The motion of an incompressible viscous fluid induced by a spinning cone is analytically studied and similar solutions of the relevant steady state boundary equations are obtained. Some of the numerical results are shown to be obtainable from the Karman-Cochran solution for the infinite disc.Symbols and Notation p Pressure - p Pressure at infinity - p ₀ Pressure at the wall - Density - Transverse component of velocity - Normal component of velocity - Radial component of velocity - Angular velocity - Semi-vertex angle - Re Reynolds number with respect to _o - _o Transverse component of velocity at the cone surface - Kinematic viscosity This research is sponsored by the Air Force Office of Scientific Research, Fluid Mechanics Division, under Contract Number AF 18(600)-498.     

Perturbation analysis for periodic heat transfer in radiating fins

A perturbation analysis is presented for periodic heat transfer in radiating fins of uniform thickness. The base temperature is assumed to oscillate around a mean value. The perturbation expansion is carried out in terms of dimensionless amplitude of the base temperature oscillation. The zero-order problem which is nonlinear, and corresponds to the steady state fin behaviour, is solved by quasilinearization. A method of complex combination is used to reduce both the first and the second order problems to two, coupled linear boundary value problems which are subsequently solved by a noniterative numerical scheme. The second-order term is composed of an oscillatory component with twice the frequency of base temperature oscillation and a time-independent term which causes a net change in the steady state values of temperature and heat transfer rate. Within the range of parameters used, the net effect is to decrease the mean temperature and increase the mean heat transfer rate. This is in constrast to the linear case of convecting fins where the mean values are unaffected by base temperature oscillations. Detailed numerical results are presented illustrating the effects of fin parameter N and dimensionless frequency B on temperature distribution, heat transfer rate, and time-average fin efficiency. The time-average fin efficiency is found to reduce significantly at low N and high B.

---

Störungsanalyse für periodische Wärmeübertragung an Strahlungsrippen

Zusammenfassung Eine Störungsanalyse wird für periodische Wärmeübertragung in Strahlungsrippen gleicher Dicke vorgelegt. Die Fußtemperatur wird als um einen Mittelwert schwingend angenommen. Die Störungsentwicklung wird in Termen einer dimensionslosen Amplitude e dieser Schwingung angesetzt. Das Problem nullter Ordnung, das nichtlinear ist und dem stationären Verhalten der Rippe entspricht, wird durch Quasilinearisierung gelöst. Eine Methode der komplexen Kombination wird angewandt, um die Probleme erster und zweiter Ordnung auf zwei gekoppelte Grenzwertprobleme zu reduzieren, die nacheinander nach einem nichtiterativen Schema gelöst werden. Der Term zweiter Ordnung besteht aus einer Schwingungskomponente mit der doppelten Frequenz der Schwingung der Fußtemperatur und einem zeitunabhängigen Term, der eine Nettoänderung der stationären Werte der Temperatur und der Wärmeübertragung verursacht. Im verwendeten Bereich der Parameter tritt eine Abnahme der mittleren Temperatur und eine Zunahme der mittleren Wärmeübertragung auf. Das steht im Gegensatz zum linearen Fall der Konvektionsrippe, bei dem die Mittelwerte durch Schwingungen der Fußtemperatur nicht beeinflußt werden. Detaillierte numerische Ergebnisse zeigen die Einflüsse des Rippenparameters N und der dimensionslosen Frequenz B auf Temperatur Verteilung, Wärmeübertragung und zeitliches Mittel des Rippengütegrades. Dieses zeitliche Mittel nimmt merklich ab bei kleinem N und hohem B.

---

Nomenclature b fin thickness - B dimensionless frequency, L²/ - E emissivity - f₀, f₁ functions of X - g₀, g₁, g₂ functions of X - h₀, h₁, h₂ functions of X - k thermal conductivity - L fin Length - N fin parameter, 2EL²T_bm/bk - q heat transfer rate - Q dimensionless heat transfer rate, qL/kbT_bm - t time - T temperature - T_b fin base temperature - T_S effective sink temperature - T_bm mean fin base temperature - x axial distance - X dimensionless axial distance, x/L - dimensionless amplitude of base temperature (s. Eq.2) - thermal diffusivity - instantaneous fin efficiency - time-average fin efficiency - _ss steady state fin efficiency - dimensionless temperature, T/T_bm - ₀ zero-order approximation - ₁ first-order approximation - ₂ second-order approximation - _2s steady component of ₂ - , ₁, ₂ constants - complex function of X - ₁ real part of - ₂ imaginary part of - complex function of X - ₁ real part of Y - ₂ imaginary part of - dimensionless time, t/L² - frequency of base temperature oscillation     

Determination of specific thrust in unchoked regimes and design of a separationless convergent nozzle contour

The variation of the specific thrust R_Y on the angle of inclination of the wall is analyzed within the framework of the ideal gas model using the results of specific impulse and flow rate calculations for conical convergent nozzles. It is shown that in unchoked regimes nozzles with different have almost the same values of R_Y for both subcritical and supercritical pressure ratios _c. On the interval _C < 6 typical of convergent nozzles conical convergent nozzles with =30–90° have almost the same value of the specific thrust, maximal relative to the R_Y of nozzles with < 30°. In the presence of viscosity forces local boundary layer separation may occur in the neighborhood of the entrance section of the convergent nozzle. A method of constructing a separationless convergent nozzle contour with enhanced thrust is developed on the basis of a boundary layer separation criterion. The separationless contour is determined for given values of the flow rate, specific heat ratio, Reynolds number, wall temperature and initial boundary layer displacement thickness.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 158–164, January–February, 1990.     

The value of the critical exponent for reaction-diffusion equations in cones

Let D R ^N be a cone with vertex at the origin i.e., D = (0, )x where S ^N–1 and x D if and only if x = (r, ) with r=¦x¦, . We consider the initial boundary value problem: u _t = u+u ^p in D×(0, T), u=0 on Dx(0, T) with u(x, 0)=u ₀(x) 0. Let ₁ denote the smallest Dirichlet eigenvalue for the Laplace-Beltrami operator on and let ₊ denote the positive root of (+N–2) = ₁. Let p ^* = 1 + 2/(N + ₊). If 1 < p < p ^*, no positive global solution exists. If p>p ^*, positive global solutions do exist. Extensions are given to the same problem for u _t=+¦x¦ u ^p .This research was supported in part by the Air Force Office of Scientific Research under Grant # AFOSR 88-0031 and in part by NSF Grant DMS-8 822 788. The United States Government is authorized to reproduce and distribute reprints for governmental purposes not withstanding any copyright notation therein.     

A rotating hot-wire technique for spatial sampling of disturbed and manipulated duct flows

The documentation and control of flow disturbances downstream of various open inlet contractions was the primary focus with which to evaluate a spatial sampling technique. An X-wire probe was rotated about the center of a cylindrical test section at a radius equal to one-half that of the test section. This provided quasi-instantaneous multi-point measurements of the streamwise and azimuthal components of the velocity to investigate the temporal and spatial characteristics of the flowfield downstream of various contractions. The extent to which a particular contraction is effective in controlling ingested flow disturbances was investigated by artificially introducing disturbances upstream of the contractions. Spatial as well as temporal mappings of various quantities are presented for the streamwise and azimuthal components of the velocity. It was found that the control of upstream disturbances is highly dependent on the inlet contraction; for example, reduction of blade passing frequency noise in the ground testing of jet engines should be achieved with the proper choice of inlet configurations.List of symbols K _uv correlation coefficient= - P percentage of time that an azimuthal fluctuating velocity derivative dv/d is found - U streamwise velocity component U=U (, t) - V azimuthal or tangential velocity component due to flow and probe rotation V=V (, t) - mean value of streamwise velocity component - U _m resultant velocity from and - mean value of azimuthal velocity component induced by rotation - u fluctuating streamwise component of velocity u=u(, t) - v fluctuating azimuthal component of velocity v = v (, t) - u phase-averaged fluctuating streamwise component of velocity u=u(0) - v phase-averaged fluctuating azimuthal component of velocity v=v() - û average of phase-averaged fluctuating streamwise component of velocity (u()) over cases I-1, II-1 and III-1 û = û() - average of phase-averaged fluctuating azimuthal component of velocity (v()) over cases I-1, II-1 and III-1 - u fluctuating streamwise component of velocity corrected for non-uniformity of probe rotation and/or phase-related vibration u = u(0, t) - v fluctuating azimuthal component of velocity corrected for non-uniformity or probe rotation and/or phase-related vibration v=v (, t) - u ² rms value of corrected fluctuating streamwise component of velocity - rms value of corrected fluctuating azimuthal component of velocity - phase or azimuthal position of X-probe     

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.     

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

Scalar Minimizers with Fractal Singular Sets

Lack of regularity of local minimizers for convex functionals with non-standard growth conditions is considered. It is shown that for every >0 there exists a function aC() such that the functional admits a local minimizer uW^1,p() whose set of non-Lebesgue points is a closed set with dim()>N–p–, and where 1<p<N<N+<q<+.     

Displacement type unified equation for bending, buckling and vibration of conical shells and its applications

By using Donnell's simplication and starting from the displacement type equations of conical shells, and introducing a displacement functionU(s,,) (In the limit case, it will be reduced to cylindrical shell displacement function introduced by V. S. Vlasov) and a generalized loadq,(s,,),the equations of conical shells are changed into an eighth—order solvable partial differential equation about the displacement functionU(s,,). As a special case, the general bending problem of conical shells on Winkler foundation has been studied. Detailed numerical results and boundary coefficients for edge unit loads are obtained.The project supported by the National Natural Science Foundation of China.     

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.