Simulation of turbulent flows around a circular cylinder using nonlinear eddy-viscosity modelling: steady and oscillatory ambient flows

Steady incident flow past a circular cylinder for sub- to supercritical Reynolds number has been simulated as an unsteady Reynolds-averaged Navier–Stokes (RANS) equation problem using nonlinear eddy-viscosity modelling assuming two-dimensional flow. The model of Craft et al. (Int. J. Heat Fluid Flow 17 (1996) 108), with adjustment of the coefficients of the ‘cubic’ terms, predicts the drag crisis at a Reynolds number of about 2×10⁵ due to the onset of turbulence upstream of separation and associated changes in Strouhal number and separation positions. Slightly above this value, at critical Reynolds numbers, drag is overestimated because attached separation bubbles are not simulated. These do not occur at supercritical Reynolds numbers and drag coefficient, Strouhal number and separation positions are in approximate agreement with experimental measurements (which show considerable scatter). Fluctuating lift predictions are similar to sectional values measured experimentally for subcritical Reynolds numbers but corresponding measurements have not been made at supercritical Reynolds numbers. For oscillatory ambient flow, in-line forces, as defined by drag and inertia coefficients, have been compared with the experimental values of Sarpkaya (J. Fluid Mech. 165 (1986) 61) for values of the frequency parameter, β=D²/νT, equal to 1035 and 11240 and Keulegan–Carpenter numbers, KC=U₀T/D, between 0.2 and 15 (D is cylinder diameter, ν is kinematic viscosity, T is oscillation period, and U₀ is the amplitude of oscillating velocity). Variations with KC are qualitatively reproduced and magnitudes show best agreement when there is separation with a large-scale wake, for which the turbulence model is intended. Lift coefficients, frequency and transverse vortex shedding patterns for β=1035 are consistent with available experimental information for β≈250−500. For β=11240, it is predicted that separation is delayed due to more prominent turbulence effects, reducing drag and lift coefficients and causing the wake to be more in line with the flow direction than transverse to it. While these oscillatory flows are highly complex, attached separation bubbles are unlikely and the flows probably two dimensional.     

Vibro-impact dynamics near a strong resonance point

Two typical vibratory systems with impact are considered, one of which is a two-degree-of-freedom vibratory system impacting an unconstrained rigid body, the other impacting a rigid amplitude stop. Such models play an important role in the studies of dynamics of mechanical systems with repeated impacts. Two-parameter bifurcations of fixed points in the vibro-impact systems, associated with 1:4 strong resonance, are analyzed by using the center manifold and normal form method for maps. The single-impact periodic motion and Poincaré map of the vibro-impact systems are derived analytically. Stability and local bifurcations of a single-impact periodic motion are analyzed by using the Poincaré map. A center manifold theorem technique is applied to reduce the Poincaré map to a two-dimensional one, and the normal form map for 1:4 resonance is obtained. Local behavior of two vibro-impact systems, near the bifurcation points for 1:4 resonance, are studied. Near the bifurcation point for 1:4 strong resonance there exist a Neimark–Sacker bifurcation of period one single-impact motion and a tangent (fold) bifurcation of period 4 four-impact motion, etc. The results from simulation show some interesting features of dynamics of the vibro-impact systems: namely, the “heteroclinic” circle formed by coinciding stable and unstable separatrices of saddles, T _in, T _on and T _out type tangent (fold) bifurcations, quasi-periodic impact orbits associated with period four four-impact and period eight eight-impact motions, etc. Different routes of period 4 four-impact motion to chaos are obtained by numerical simulation, in which the vibro-impact systems exhibit very complicated quasi-periodic impact motions. The project supported by National Natural Science Foundation of China (50475109, 10572055), Natural Science Foundation of Gansu Province Government of China (3ZS061-A25-043(key item)). The English text was polished by Keren Wang.     

Order Reduction of Parametrically Excited Nonlinear Systems: Techniques and Applications

The basic problem of order reduction of nonlinear systems with time periodic coefficients is considered in state space and in direct second order (structural) form. In state space order reduction methods, the equations of motion are expressed as a set of first order equations and transformed using the Lyapunov–Floquet (L–F) transformation such that the linear parts of new set of equations are time invariant. At this stage, four order reduction methodologies, namely linear, nonlinear projection via singular perturbation, post-processing approach and invariant manifold technique, are suggested. The invariant manifold technique yields a unique ‘reducibility condition’ that provides the conditions under which an accurate nonlinear order reduction is possible. Unlike perturbation or averaging type approaches, the parametric excitation term is not assumed to be small. An alternate approach of deriving reduced order models in direct second order form is also presented. Here the system is converted into an equivalent second order nonlinear system with time invariant linear system matrices and periodically modulated nonlinearities via the L–F and other canonical transformations. Then a master-slave separation of degrees of freedom is used and a nonlinear relation between the slave coordinates and the master coordinates is constructed. This method yields the same ‘reducibility conditions’ obtained by invariant manifold approach in state space. Some examples are given to show potential applications to real problems using above mentioned methodologies. Order reduction possibilities and results for various cases including ‘parametric’, ‘internal’, ‘true internal’ and ‘true combination resonances’ are discussed. A generalization of these ideas to periodic-quasiperiodic systems is included and demonstrated by means of an example.     

Non-linear oscillation of circular cylindrical shells

The method of multiple scales is used to analyze the non-linear forced response of circular cylindrical shells in the presence of a two-to-one internal (autoparametric) resonance to a harmonic excitation having the frequency Ω. If ω_r and a_r denote the frequency and amplitude of a flexural mode and ω_b and a_b denote the frequency and amplitude of the breathing mode, the steady-state response exhibits a saturation phenomenon when ω_b ≈ 2ω_r, if the excitation frequency Ω is near ω_b. As the amplitude ƒ of the excitation increases from zero, a_b increases linearly whereas a_r remains zero until a threshold is reached. This threshold is a function of the damping coefficients and ω_b−2ω_r. Beyond this threshold a_b remains constant (i.e. the breathing mode saturates) and the extra energy spills over into the flexural mode. In other words, although the breathing mode is directly excited by the load, it absorbs a small amount of the input energy (responds with a small amplitude) and passes the rest of the input energy into the flexural mode (responds with a large amplitude). For small damping coefficients and depending on the detunings of the internal resonance and the excitation, the response exhibits a Hopf bifurcation and consequently there are no steadystate periodic responses. Instead, the responses are amplitude- and phase-modulated motions. When Ω ≈ ω_r, there is no saturation phenomenon and at close to perfect resonance, the response exhibits a Hopf bifurcation, leading again to amplitude- and phase-modulated or chaotic motions.     

Linear Stability Analysis of Resonant Periodic Motions in the Restricted Three-Body Problem

The equations of the restricted three-body problem describe the motion of a massless particle under the influence of two primaries of masses 1 −μ and μ, 0≤μ≤ 1/2, that circle each other with period equal to 2π. When μ=0, the problem admits orbits for the massless particle that are ellipses of eccentricity e with the primary of mass 1 located at one of the focii. If the period is a rational multiple of 2π, denoted 2π p/q, some of these orbits perturb to periodic motions for μ > 0. For typical values of e and p/q, two resonant periodic motions are obtained for μ > 0. We show that the characteristic multipliers of both these motions are given by expressions of the form in the limit μ→ 0. The coefficient C(e,p,q) is analytic in e at e=0 and C(e,p,q)=O(e^|p-q|). The coefficients in front of e^|p-q|, obtained when C(e,p,q) is expanded in powers of e for the two resonant periodic motions, sum to zero. Typically, if one of the two resonant periodic motions is of elliptic type the other is of hyperbolic type. We give similar results for retrograde periodic motions and discuss periodic motions that nearly collide with the primary of mass 1 −μ.     

Intermittency and Regularity Issues in 3<Emphasis Type="Italic">D</Emphasis> Navier-Stokes Turbulence

Two related open problems in the theory of 3D Navier-Stokes turbulence are discussed in this paper. The first is the phenomenon of intermittency in the dissipation field. Dissipation-range intermittency was first discovered experimentally by Batchelor and Townsend over fifty years ago. It is characterized by spatio-temporal binary behaviour in which long, quiescent periods in the velocity signal are interrupted by short, active ‘events’ during which there are violent fluctuations away from the average. The second and related problem is whether solutions of the 3D Navier-Stokes equations develop finite time singularities during these events. This paper shows that Leray’s weak solutions of the three-dimensional incompressible Navier-Stokes equations can have a binary character in time. The time-axis is split into ‘good’ and ‘bad’ intervals: on the ‘good’ intervals solutions are bounded and regular, whereas singularities are still possible within the ‘bad’ intervals. An estimate for the width of the latter is very small and decreases with increasing Reynolds number. It also decreases relative to the lengths of the good intervals as the Reynolds number increases. Within these ‘bad’ intervals, lower bounds on the local energy dissipation rate and other quantities, such as ||u(·, t)||_∞ and ||∇u(·, t)||_∞, are very large, resulting in strong dynamics at sub-Kolmogorov scales. Intersections of bad intervals for n≧1 are related to the potentially singular set in time. It is also proved that the Navier-Stokes equations are conditionally regular provided, in a given ‘bad’ interval, the energy has a lower bound that is decaying exponentially in time.Final version 17 March 05. Original version November 03.     

Codimension two bifurcation of a vibro-bounce system

A three-degree-of-freedom vibro-bounce system is considered. The disturbed map of period one single-impact motion is derived analytically. A center manifold theorem technique is applied to reduce the Poincaré map to a three-dimensional one, and the normal form map associated with Hopf-flip bifurcation is obtained. Dynamical behavior of the system, near the point of codimension two bifurcation, is investigated by using qualitative analysis and numerical simulation. It is found that near the point of Hopf-flip bifurcation there exists not only Hopf bifurcation of period one single-impact motion, but also Hopf bifurcation of period two double-impact motion. The results from simulation show that there exists an interesting torus doubling bifurcation near the codimension two bifurcation. The torus doubling bifurcation makes the quasi-periodic attractor associated with period one single-impact motion transform to the other quasi-periodic attractor represented by two attracting closed circles. The torus bifurcation is qualitatively different from the typical torus doubling bifurcation occurring in the vibro-impact systems. Different routes from period one single-impact motion to chaos are observed by numerical simulation.The project supported by the National Natural Science Foundation of China (10172042, 50475109) and the Natural Science Foundation of Gansu Province Government of China (ZS-031-A25-007-Z (key item))     

The mechanisms of ‘neck-like’ deformation in high-speed melt spinning. 2. Effects of polymer crystallization

As is shown in the first paper of the series, the main factor responsible for concentrated (‘neck-like’) deformation in high-speed melt spinning is the gradient of elongational viscosity along the spinline. In the present paper, stress-induced polymer crystallization is analyzed as a potential source of the rapid viscosity increase.A model of crystallization-controlled solidification is proposed, in which viscosity of the polymer increases with the degree of crystallinity, Θ, as

Full-size image (1K)

, reaching infinity (complete solidification) at Θ = Θ_cr. The critical crystallinity level has been interpreted as one required for ‘crosslinking’ of polymer chains present in the melt.In addition to viscosity increase, crystallization modifies the local temperature in the spinline and reduces viscosity.The analysis of stress effects shows that critical crystallization temperature, T_m, and crystallization rate, K, increase with the square of normal stress difference in the spinline, Δp = p_xx − p_rr. The onset of crystallization can be shifted by 20–40 K towards higher temperatures, and crystallization rate can increase by orders of magnitude when high take-up speeds increase the stress level.A simple model illustrating velocity profiles in crystallizing Newtonian jets is discussed.The analysis strongly supports the hypothesis that the high viscosity gradient resulting from rapid stress-induced crystallization provides the major mechanism of ‘neck-like’ deformation.     

Period doubling bifurcation problems in the softening Duffing oscillator with nonstationary excitation

The Duffing oscillators are widely used to mathematically model a variety of engineering and physical systems. A computational analysis has been initiated to explore the effects of nonstationary excitations on the response of the softening Duffing oscillator in the region of the parameter space where the period doubling sequences occur. Significant differences between the stationary and nonstationary responses have been uncovered: (i) the stationary transitions from T to 2T, from 2T to 4T ... etc. branches at the stationary period doubling bifurcations are smooth, in nonstationary cases they exhibit jumps to the near stationary branches at the values of the control parameters greater than those for the stationary; this phenomenon is called penetration (delay or memory). The lengths of the penetrations is being compressed to zero with the increasing number of the iterations. (ii) The stationary and nonstationary responses eventually settle on different limit motions, the nonstationary has modulated components. (iii) The jumps appearing in the stationary bifurcation diagram at 2T from the upper to the lower branches of the (x, f) and (x, ), i.e., (displacement-forcing amplitude) and (displacement-forcing frequency), diagrams have been replaced by continuous transition in the nonstationary diagram climinating thus the discontinuity. Apart from these differences, some specific characteristic nonstationary responses have been observed not encountered in the stationary cases: (iv) the appearance of the window in the nonstationary limit bifurcation diagrams. (v) The nonstationary limit motions located on the upper (lower) branches of the (x, f) or (x, ) diagrams expanded rapidly to the lower (upper) branches. (vi) The stationary and nonstationary bifurcation diagrams are extremely sensitive to the initial conditions, manifested by the mirror reflections, called the flipflop phenomenon. (vii) The nonstationary limit motion has been characterized by a complex phase portrait, the appearance of the Cantor-like set of the limit motion bifurcation plot, and continuous spectral density. For the purpose of comparison, a stationary period doubling sequence T, 2T,..., 2 ⁿ T,... stationary limit motion, _ST which is known to be chaotic has been determined. A far reaching observation has been made in the process of this study: the determination of the nonstationary bifurcations, their branches and limit motions, has been independent of the calculations of the stationary ones, indicating, thus, the existence of an independent class of nonstationary (time varying) dynamics.     

Translational diffusion of rubrene and tetracene in polyisobutylene

The translational diffusion coefficients of rubrene and tetracene in amorphous polyisobutylene (PIB) were measured using the holographic fluorescence recovery after photobleaching technique. Over the temperature range from 400 to 235 K (T _g =205 K), tracer diffusion coefficients from 10^–7 to 10^–14 cm²/s were observed. These diffusion coefficients have essentially the same temperature dependence as the rotational correlation times for these two probes in PIB. Both of these observables have a slightly stronger temperature dependence than does the viscosity. These results contrast strongly with the results of similar experiments on polystyrene and polysulfone. These results are consistent with the hypothesis that local segmental dynamics are more spatially homogeneous in PIB than in polystyrene and polysulfone.Dedicated to Prof. John D. Ferry on the occasion of his 85th birthday.     

Symmetry and secondary bifurcation of elastic structures

We consider non-linear bifurcation problems for elastic structures modeled by the operator equation F[w;α]=0 where F:X×R^k→Y,X,Y are Banach spaces and XY. We focus attention on problems whose bifurcation equations are of the form

f_i(α₁,α₂;λ,μ)=(a_iμ+b_iλ)α_i+p_iα_i³+q_iα_i∑_j=1,j≠i^kα_j+1²+α_ih_i(λ,μ;α₁,α₂,…α_k) i=1,2,…k

which emanates from bifurcation problems for which the linearization of F is Fredholm operators of index 0. Under the assumption of F being odd we prove an important theorem of existence of secondary bifurcation. Under this same assumption we prove a symmetry condition for the reduced equations and consequently we got an existence result for secondary bifurcation. We also include a stability analysis of the bifurcating solutions.     

Fatigue characteristics of ‘bead-on-plate’ of aluminum and Al---Zn---Mg alloy

Fatigue life of ‘bead-on-plate’ on commercially pure aluminum and Al---Zn---Mg alloy have been investigated. In both cases, the fatigue life of ‘Bead-on-Plate’ specimen was drastically reduced. But when the same weld-bead was flushed off, a marginal improvement in fatigue life was observed. For Al---Zn---Mg alloy, the post weld heat treatment was less effective than flushing of weld-bead for improving the fatigue life of welded specimens. Results are analyzed with the help of S-N diagrams and microphotographs.     

Scaling law in saddle-node bifurcations for one-dimensional maps: a complex variable approach

The study of transient dynamical phenomena near bifurcation thresholds has attracted the interest of many researchers due to the relevance of bifurcations in different physical or biological systems. In the context of saddle-node bifurcations, where two or more fixed points collide annihilating each other, it is known that the dynamics can suffer the so-called delayed transition. This phenomenon emerges when the system spends a lot of time before reaching the remaining stable equilibrium, found after the bifurcation, because of the presence of a saddle-remnant in phase space. Some works have analytically tackled this phenomenon, especially in time-continuous dynamical systems, showing that the time delay, τ, scales according to an inverse square-root power law, τ∼(μ−μ _c )^−1/2, as the bifurcation parameter μ, is driven further away from its critical value, μ _c . In this work, we first characterize analytically this scaling law using complex variable techniques for a family of one-dimensional maps, called the normal form for the saddle-node bifurcation. We then apply our general analytic results to a single-species ecological model with harvesting given by a unimodal map, characterizing the delayed transition and the scaling law arising due to the constant of harvesting. For both analyzed systems, we show that the numerical results are in perfect agreement with the analytical solutions we are providing. The procedure presented in this work can be used to characterize the scaling laws of one-dimensional discrete dynamical systems with saddle-node bifurcations.     

A new transient method for analysis of solidification in the continuous casting of metals

Solidification processes involve complex heat and mass transfer phenomena, the modelling of which requires state-of-the art numerical techniques. An efficient and accurate transient numerical method is proposed for the analysis of phase change problems. This method combines both the enthalpy and the enhanced specific heat approaches in incorporating the effects of latent heat released due to phase change. The sensitivity and accuracy of the proposed method to both temporal and spatial discretization is shown together with closed-form solutions and the results from the enhanced specific heat approach. In order to explore the proposed method fully, a non-linear heat release, as is the case for binary alloys, is also examined. The number of operations required for the new transient approach is less than or equal to the enhanced heat capacity method depending on the averaging method adopted. To demonstrate the potential of this new finite-element technique, measurements obtained on operating machines for the casting of zinc, aluminum and steel are compared with the model predictions. The death/birth technique, together with the proper heat-transfer coefficients, were employed in order to model the casting process with minimal error due to the modelling itself.Nomenclature [A] conductance matrix - [B] matrix containing the derivative of the element shape functions - c, C _p specific heat (J kg^–1°C^–1) - effective specific heat (J kg^–1°C^–1) - f(T) local liquid fraction - f thermal load vector - H enthalpy (J kg^–1) - [H] capacitance matrix - h, h _r,h _c heat transfer coefficient (W m^–2°C^–1) - K thermal conductivity (W m^–1°C^–1) - L latent heat of solidification (J kg^–1) - l overall length (m) - N _i shape functions - Q rate of heat generation per unit volume (J m^–3) - q heat flux (W m^–2) - R residual temperature (°C) - T temperature (°C) - T _s solidus temperature (°C) - T _l liquidus temperature (°C) - T _pouring pouring temperature (°C) - T _top temperature at the top of the mould (°C) - T _w temperature of the water spray (°C) - approximated temperature (°C) - T surrounding temperature (°C) - cooling rate (°C/s) - t time (seconds) - x _i,x, y, z spatial variables (m) - t time step (s) - x element size (m) - diffusivity (m²s^–1) - density (kg m^–3) - time marching parameter - _n direction cosines of the unit outward normal to the boundary     

Theorem on perturbation of coisotropic invariant tori of locally Hamiltonian systems and its applications

We study the problem of perturbations of quasiperiodic motions on coisotropic invariant tori in a class of locally Hamiltonian systems. We prove a general KAM-theorem on the perturbation of coisotropic invariant tori for locally Hamiltonian systems. As applications of this theorem, we consider the motion of an electron on a two-dimensional torus under the action of an electromagnetic field and extend results concerning the bifurcation of a Cantor set of coisotropic invariant tori to the case of locally Hamiltonian systems. __________ Translated from Neliniini Kolyvannya, Vol. 8, No. 4, pp. 490–515, October–December, 2005.     

Upper limit of enhancement of boiling heat transfer in a narrow vertical rectangular channel due to passing bubbles

An experimental study has been made of saturated boiling heat transfer for water and R113 in a narrow vertical rectangular channel (2 mm space, 20 mm wide, and 200 mm long) at atmospheric pressure, in which the vertical heated surface (10 mm long and 20 mm wide) is located on one side at a position of 150 mm from its entrance and bubbles are forcibly passed through it at a designated period from 0.33 to 1.0 sec. The experiment shows that the heat transfer coefficients are increased by the bubble passing through the heated surface for the value of thermal diffusivity,a, times period, T₀, of the passing bubbles above about 6×10^–9 m² (a T ₀>6×10^–9 m²) while fora T ₀< 6×10^–9 m², the heat transfer coefficients become independent of the period and the effectiveness of the enhancement of the heat transfer owing to the passing bubble disappears.

---

Die obere Grenze der Verbesserung des Wärmeübergangs beim Sieden in einem vertikalen, rechteckigen Kanal infolge von aufsteigenden Blasen

Zusammenfassung Es wurden Experimente über den Wärmeübergang beim Sättigungssieden mit Wasser und R113 in einem engen, vertikalen, rechteckigen Kanal (2 mm Abstand, 20 mm Breite und 200 mm Länge) bei Umgebungsdruck durchgeführt, wobei die vertikale, beheizte Oberfläche (10 mm lang und 20 mm breit) auf der einen Seite in einem Abstand von 150 mm vom Eintritt angeordnet ist und die Blasen zwangsweise durch den Kanal sich mit einem Periodenabstand von 0,033 bis 1,0 s bewegen. Das Experiment zeigt, daß die Wärmeübergangskoeffizienten durch das Vorbeistreichen der Blasen an der beheizten Oberfläche verbessert werden, wenn das Produkt aus Temperaturleitfähigkeit,a, mal der Periode, T₀, der vorbeistreichenden Blasen größer als 6×10^–9 m² liegt, während unterhalb dieses Wertes der Wärmeübergangskoeffizient unabhängig von der Blasenperiode ist und die Effektivität der Wärmeübergangsverbesserung infolge der Blasenströmung verschwindet.

---

Nomenclature a thermal diffusivity of liquid - ¯h time-averaged heat transfer coefficient - q _w heat flux at wall - T ₀ period of passing bubble - T _w(t) temperature of heated surface - T _w amplitude of heated surface temperature Greek symbols thermal conductivity - thickness of liquid film     

On the laminar forced convection with axial conduction in a circular tube with exponential wall heat flux

The values of the fully developed Nusselt number for laminar forced convection in a circular tube with axial conduction in the fluid and exponential wall heat flux are determined analytically. Moreover, the distinction between the concepts of bulk temperature and mixing-cup temperature, at low values of the Peclet number, is pointed out. Finally it is shown that, if the Nusselt number is defined with respect to the mixing-cup temperature, then the boundary condition of exponentially varying wall heat flux includes as particular cases the boundary conditions of uniform wall temperature and of convection with an external fluid.

---

Über laminare Zwangskonvektion mit Längswärmeleitung in einem Kreisrohr mit exponentiell veränderlichem Wandwärmefluß

Zusammenfassung Es werden die Endwerte der Nusselt-Zahlen für vollausgebildete laminare Zwangskonvektion in einem Kreisrohr mit Längswärmeleitung und exponentiell veränderlichem Wandwärmefluß analytisch ermittelt. Besondere Betonung liegt auf dem Unterschied zwischen den Konzepten für die Mittel- und die Mischtemperatur bei niedrigen Peclet-Zahlen. Schließlich wird gezeigt, daß bei Definition der Nusselt-Zahl bezüglich der Mischtemperatur die Randbedingung exponentiell veränderlichen Randwärmeflusses die Spezialfälle konstanter Wandtemperatur und konvektiven Wärmeaustausches mit einem umgebenden Fluid einschließt.

---

Nomenclature A _n dimensionless coefficients employed in the Appendix - Bi Biot numberBi=h _e r ₀/ - c _n dimensionless coefficients defined in Eq. (17) - c _p specific heat at constant pressure of the fluid within the tube, [J kg^–1 K^–1] - f solution of Eq. (15) - h ₁,h ₂ specific enthalpies employed in Eqs. (2) and (4), [J kg^–1] - h _e convection coefficient with a fluid outside the tube, [W m^–2 K^–1] - rate of mass flow, [kg s^–1] - Nu bulk Nusselt number,2r ₀ q _w /[(T _w –T _b )] - Nu _H fully developed value of the bulk Nusselt number for the boundary condition of uniform wall heat flux - Nu _T fully developed value of the bulk Nusselt number for the boundary condition of uniform wall temperature - Nu ^* mixing Nusselt number,2r ₀ q _w /[(T _w –T _m )] - Nu _C ^* fully developed value of the mixing Nusselt number for the boundary condition of convection with an external fluid - Nu _H ^* fully developed value of the mixing Nusselt number for the boundary condition of uniform wall heat flux - Nu _T ^* fully developed value of the mixing Nusselt number for the boundary condition of uniform wall temperature - Pe Peclet number, 2r ₀/ - q ₀ wall heat flux atx=0, [W m^–2] - q _w wall heat flux, [W m^–2] - r radial coordinate, [m] - r ₀ radius of the tube, [m] - s dimensionless radius,s=r/r ₀ - T temperature, [K] - T ₀ temperature constant employed in Eq. (14), [K] - T reference temperature of the fluid external to the tube, [K] - T _b bulk temperature, [K] - T _m mixing or mixing-cup temperature, [K] - T _w wall temperature, [K] - u velocity component in the axial direction, [m s^–1] - mean value ofu, [m s^–1] - x axial coordinate, [m] Greek symbols thermal diffusivity of the fluid within the tube, [m² s^–1] - exponent in wall heat flux variation, [m^–1] - dimensionless parameter - dimensionless temperature =(T _w –T)/(T _w –T _b ) - ^* dimensionless temperature ^*=(T _w –T)/(T _w –T _m ) - thermal conductivity of the fluid within the tube, [W m^–1 K^–1] - density of the fluid within the tube, [kg m^–3]     

On the existence of low-period orbits in <Emphasis Type="Italic">n</Emphasis>-dimensional piecewise linear discontinuous maps

Discontinuous maps occur in many practical systems, and yet bifurcation phenomena in such maps is quite poorly understood. In this paper, we report some important results that help in analyzing the border collision bifurcations that occur in n-dimensional discontinuous maps. For this purpose, we use the piecewise linear approximation in the neighborhood of the plane of discontinuity. Earlier, Feigin had made a similar analysis for general n-dimensional piecewise smooth continuous maps. In this paper, we extend that line of work for maps with discontinuity to obtain the general conditions of existence of period-1 and period-2 fixed points before and after a border collision bifurcation. The application of the method is then illustrated using a specific example of a two-dimensional discontinuous map. This work was supported in part by the BRNS, Department of Atomic Energy (DAE), Government of India under project no. 2003/37/11/BRNS.     

Rayleigh waves obtained by the indeterminate couple-stress theory

Rayleigh waves in a linear elastic couple-stress medium are investigated; the constitutive equations involve a length parameter l that characterizes the microstructure of the material. With , c_T=conventional transversal speed and q=wave number, an explicit expression is derived for the relation between , lq and Poisson's ratio ν. The Rayleigh speed turns out to be dispersive and always larger than the conventional Rayleigh speed. It is of interest that when lq=1 and ν≥0, it always holds that . The displacement field is investigated and it is shown that no Rayleigh wave motions exist when lq→∞ and when lq=1, ν≥0. Moreover, a principal change of the displacement field occurs when lq passes unity. The peculiarity that no Rayleigh wave motions exist when lq=1, ν≥0 may support the criticism by Eringen (1968) against the couple-stress theory adopted here as well as in much recent literature.