Experimental/analytical characterization of composite tubes under combined loading

The paper outlines the results of an investigation to characterize the response of P75/934 graphite/epoxy tubes with a stacking sequence of [15/0/±10/0/−15] _s under pure torsion and combined axial/torsion loading. The experimentally observed nonlinear response and path-dependent failure are discussed in terms of material nonlinearities at the ply level and first-ply failure loads with the help of an analytical model. Paper was presented at the 1989 SEM Spring Conference on Experimental Mechanics held in Cambridge, MA on May 28–June 1.     

On Extension and Torsion of Strain-Stiffening Rubber-Like Elastic Circular Cylinders

This paper is concerned with investigation of the effects of strain-stiffening on the response of solid circular cylinders in the combined deformation of torsion superimposed on axial extension. The cylinders are composed of incompressible isotropic nonlinearly elastic materials. Our primary focus is on materials that undergo severe strain-stiffening in the stress-stretch response. In particular, we consider two particular phenomenological constitutive models for such materials that reflect limiting chain extensibility at the molecular level. The axial stretch γ and twist that can be sustained in cylinders composed of such materials are shown to be constrained in a coupled fashion. It is shown that, in the absence of an additional axial force, a transition value γ=γ _t of the axial stretch exists such that for γ<γ _t , the stretched cylinder tends to elongate on twisting whereas for γ>γ _t , the stretched cylinder tends to shorten on twisting. These results are in sharp contrast with those for classical models such as the Mooney-Rivlin (and neo-Hookean) models that predict that the stretched circular cylinder always tends to further elongate on twisting. We also obtain results for materials modeled by the well-known exponential strain-energy widely used in biomechanics applications. This model reflects a strain-stiffening that is less abrupt than that for the limiting chain extensibility models. Surprisingly, it turns out that the results in this case are somewhat more complicated. For a fixed stiffening parameter, provided that the stretch is sufficiently small, the stretched bar always tends to elongate on twisting in the absence of an additional axial force. However, for sufficiently large stretch, the cylinder tends to shorten on undergoing sufficiently small twist but then tends to elongate on further twisting. These results are of interest in view of the widespread use of exponential models in the context of the mechanics of soft biological tissues. The special case of pure torsion is also briefly considered. In this case, the resultant axial force required to maintain pure torsion is compressive for all the models discussed here. In the absence of such a force, the bar would elongate on twisting reflecting the celebrated Poynting effect.      

Gaseous nitromethane and nitromethane-oxygen mixtures: A new detonation structure

Detonation in gaseous nitromethane (NM) and mixed with O₂ has been studied. Experiments were performed in a preheated steel tube at an initial temperatureT ₀∼=390 K for different initial pressuresP ₀ (1.7≥P ₀≥5 10⁻² bar). Different selfsustained detonation regimes were obtained, from multiheaded mode to spinning and galloping mode in marginal conditions. These chemical systems were characterized by a specific detonation cellular structure very different from that currently observed with classical gaseous C _n H _m /O₂/N₂ mixtures. All modes of detonation propagation in rich NM/O₂ mixtures exhibit a double scale cellular structure. The pattern of this double scale structure is particularly clear in the case of spinning mode. An abridged version of this paper was presented at the 15th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Boulder, Colorado, from July 30 to August 4, 1995     

Behavior of Pin-loaded Laminated Composites

In this study, an investigation was carried out to determine the effects of joint geometry and fiber orientation on the failure strength and failure mode in a pinned joint laminated composite plate. Behavior of pin-loaded laminated composites with different stacking sequence and different dimensions has been observed experimentally. E/glass–epoxy composites were manufactured to fabricate the specimens. Mechanical properties of the composites were characterized under tension, compression and in-plane shear in static loading conditions. Laminated composites were loaded through pins. Single-hole pin-loaded specimens were tested for their tensile response and width-to-hole diameter (W/D) and edge distance-to-hole diameter (E/D) ratios evaluated. A series of experiments was performed with six different material configurations ([0/±45]_s–[90/±45]_s, [0/90/0]_s–[90/0/90]_s and [90/0]_2s–[±45]_2s), in all, over 120 specimens. E/D ratios and W/D ratios of plates were changed from 1 to 5 and 2 to 5, respectively. Failure propagation and failure type were observed on the specimens. The influence of the joint geometry on the strength of the pin-loaded composites was assessed. When laminated composite plates were loaded to final failure, three basic failure modes consisting of net-tension, shear out and bearing failure were observed for the different geometric dimensions. All the connections tested showed that the fiber orientations have a definite influence on the position around hole circumference at which failure initiated. Net-tension failure occurred for specimens that had small width and large end distance. When the width was increased, the specimens which had small end distances failed in the shear-out modes. When the end distance was increased, bearing failure developed in addition to shear-out failure. The experimental results showed that the ultimate load capacities of E/glass–epoxy laminate plates with pin connection were increased by increasing W and E. However, increasing the E/D and W/D ratios beyond a critical value has an insignificant effect on the ultimate load capacity of the connection.     

Simple shear in nonlinear Cosserat elasticity: bifurcation and induced microstructure

We discuss the simple shear problem for a geometrically exact Cosserat model. In contrast to linear Cosserat elasticity, where the unique solution is available in closed form we exhibit a multitude of solutions to the nonlinear problem, even if the two fields of deformations φ and microrotations remain homogeneous. This motivates a search for new conditions on the microrotations which single out a unique, physically reasonable, response. The influence of material parameters, notably the Cosserat couple modulus μ _c and the internal length scale L _c on the response is also studied. For small Cosserat couple modulus μ _c  > 0 we observe a pitchfork bifurcation of the homogeneous response and for vanishing internal length L _c  = 0 and zero Cosserat couple modulus μ _c  = 0 the Cosserat model may show highly oscillating “microstructure” solutions which are energetically better than the homogeneous response. Thus, the large scale nonlinear Cosserat limit is not necessarily a classical limit.      

The Asymptotic Finite-dimensional Character of a Spectrally-hyperviscous Model of 3D Turbulent Flow

We obtain attractor and inertial-manifold results for a class of 3D turbulent flow models on a periodic spatial domain in which hyperviscous terms are added spectrally to the standard incompressible Navier–Stokes equations (NSE). Let P _m be the projection onto the first m eigenspaces of A =−Δ, let μ and α be positive constants with α ≥3/2, and let Q _m =I − P _m , then we add to the NSE operators μ A _φ in a general family such that A _φ≥Q _m A ^α in the sense of quadratic forms. The models are motivated by characteristics of spectral eddy-viscosity (SEV) and spectral vanishing viscosity (SVV) models. A distinguished class of our models adds extra hyperviscosity terms only to high wavenumbers past a cutoff λ _m0 where m ₀ ≤ m, so that for large enough m ₀ the inertial-range wavenumbers see only standard NSE viscosity. We first obtain estimates on the Hausdorff and fractal dimensions of the attractor (respectively and ). For a constant K _α on the order of unity we show if μ ≥ ν that and if μ ≤ ν that where ν is the standard viscosity coefficient, l ₀ = λ₁^−1/2 represents characteristic macroscopic length, and is the Kolmogorov length scale, i.e. where is Kolmogorov’s mean rate of dissipation of energy in turbulent flow. All bracketed constants and K _α are dimensionless and scale-invariant. The estimate grows in m due to the term λ _m /λ₁ but at a rate lower than m ^3/5, and the estimate grows in μ as the relative size of ν to μ. The exponent on is significantly less than the Landau–Lifschitz predicted value of 3. If we impose the condition , the estimates become for μ ≥ ν and for μ ≤ ν. This result holds independently of α, with K _α and c _α independent of m. In an SVV example μ ≥ ν, and for μ ≤ ν aspects of SEV theory and observation suggest setting for 1/c within α orders of magnitude of unity, giving the estimate where c _α is within an order of magnitude of unity. These choices give straight-up or nearly straight-up agreement with the Landau–Lifschitz predictions for the number of degrees of freedom in 3D turbulent flow with m so large that (e.g. in the distinguished-class case for m ₀ large enough) we would expect our solutions to be very good if not virtually indistinguishable approximants to standard NSE solutions. We would expect lower choices of λ _m (e.g. with a > 1) to still give good NSE approximation with lower powers on l ₀/l _ε, showing the potential of the model to reduce the number of degrees of freedom needed in practical simulations. For the choice , motivated by the Chapman–Enskog expansion in the case m = 0, the condition becomes , giving agreement with Landau–Lifschitz for smaller values of λ _m then as above but still large enough to suggest good NSE approximation. Our final results establish the existence of a inertial manifold for reasonably wide classes of the above models using the Foias/Sell/Temam theory. The first of these results obtains such an of dimension N > m for the general class of operators A _φ if α > 5/2. The special class of A _φ such that P _m A _φ = 0 and Q _m A _φ ≥ Q _m A ^α has a unique spectral-gap property which we can use whenever α ≥ 3/2 to show that we have an inertial manifold of dimension m if m is large enough. As a corollary, for most of the cases of the operators A _φ in the distinguished-class case that we expect will be typically used in practice we also obtain an , now of dimension m ₀ for m ₀ large enough, though under conditions requiring generally larger m ₀ than the m in the special class. In both cases, for large enough m (respectively m ₀), we have an inertial manifold for a system in which the inertial range essentially behaves according to standard NSE physics, and in particular trajectories on are controlled by essentially NSE dynamics.      

Global stability of a virus dynamics model with Beddington–DeAngelis incidence rate and CTL immune response

In this paper, the global stability of virus dynamics model with Beddington–DeAngelis infection rate and CTL immune response is studied by constructing Lyapunov functions. We derive the basic reproduction number R ₀ and the immune response reproduction number R ⁰ for the virus infection model, and establish that the global dynamics are completely determined by the values of R ₀. We obtain the global stabilities of the disease-free equilibrium E ₀, immune-free equilibrium E ₁ and endemic equilibrium E ^∗ when R ₀≤1, R ₀>1, R ⁰>1, respectively.     

Detonability of THDCPD-exo–air mixtures

In the frame of industrial risk and propulsive application, the detonability study of JP10–air mixtures was performed. The simulation and measurements of detonation parameters were performed for THDCPD-exo/air mixtures at various initial pressure (1 bar < P ₀ < 3 bar) and equivalence ratio (0.8 < Φ < 1.6) in a heated tube (T ₀ ~ 375 K). Numerical simulations of the detonation were performed with the STANJAN code and a detailed kinetic scheme of the combustion of THDCPD. The experimental study deals with the measurements of detonation velocity and cell size λ. The measured velocity is in a good agreement with the calculated theoretical values. The cell size measurements show a minimum value for Φ ~ 1.2 at every level of initial pressure studied and the calculated induction length L _i corresponds to cell size value with a coefficient k = λ/L _i = 24 at P ₀ = 1 bar. Based on the comparison between the results obtained during this study and those available in the literature on the critical initiation energy E _c, critical tube diameter d _c and deflagration to detonation transition length L _DDT, we can conclude that the detonability of THDCPD–air mixtures corresponds to that of hydrocarbon–air mixtures.

---

This paper is based on the work presented at the 33rd International Pyrotechnics Seminar, IPS 2006, Fort Collins, July 16–21, 2006.     

Effect of different sized transverse square grooves on a turbulent boundary layer

The effect of three different sized transverse square grooves (5, 10, and 20 mm) on a turbulent boundary layer was investigated at two values of momentum thickness Reynolds numbers (R _θ =1,000 and 3,000) using hot-wire anemometry. The ratios of the groove depth to the boundary layer thickness (d/δ ₀) are approximately 0.07, 0.13, and 0.27. Wall shear stress (τ _w), mean velocity (U), and turbulence intensity downstream of the grooves are compared to those on a corresponding smooth wall The effects of the grooves are more significant at the higher R _θ , with the most pronounced effects caused by the largest size groove. There is an increase in mean velocity (U), streamwise (u′/U ₀), and wall-normal (ν′/U ₀) turbulence intensities in the near-wall region immediately downstream of the grooves. The increase propagates outwards in the layer as the streamwise distance increases downstream of the grooves. The increase in ν′/U ₀ is much more significant than that of u′/U ₀, which is also evident in the spectra of u′ and ν′. There is an increase in τ _w over the smooth wall value immediately downstream of the grooves at R _θ =1,000, with the increase being more pronounced as the groove size increases. The growth of the internal layer downstream of the grooves is found to scale with the groove size, and is more rapid at R _θ =3,000. Published online: 23 November 2002     

Effects of Large Scales Motion in Unstable Stratified Shear Flows

Homogeneous turbulence under unstable uniform stratification (N ² < 0) and vertical shear is investigated by using the linear theory (or the so-called rapid distortion theory, RDT) for an initial isotropic turbulence over a range −∞ ≤ R _i =N ²/S ² ≤ 0. The initial potential energy is zero and P _r =1 (i.e. the molecular Prandtl number).One-dimensional (streamwise) k ₁−spectra, especially Θ₃₃(k ₁) (i.e., that of the vertical kinetic energy, are investigated. In agreement with previous experiments, it is found that the unstable stratification affects the turbulence quantities at all scales. A significant increase of the vertical kinetic energy is observed at low wavenumbers k ₁ (i.e. large scales) due to an increase of the stratification . The effect of the shear (S) is appreciable only at high wavenumbers k ₁ (i.e. small scales).Based on the importance of the spectral components with k ₁ = 0, the asymptotic forms of Θ _ij (k ₁ = 0) or equivalently the so-called “two-dimensional” energy components (2DEC) are analyzed in detail. The asymptotic form for the ratio of 2DEC is compared to the long-time limit of the ratio of real energies. In the unstable stratified shearless case (S=0,N ² ≠ 0) the variances and the covariances of the velocity and the density are derived analytically in terms of the Weber functions, while when S ≠ 0 and N ² ≠ 0 they are obtained numerically (−100 ≤ R _i <0 and . The results are discussed in connection to previous experimental results in unstable stratified open channel flows cooled from above by Komori et al. Phy Fluids 25, 1539–1546 (1982).It is shown that the Richardson number dependence of the long-time limit of the ratios of real energies is well described by this “simple” model (i.e. the dependence of the long-time limit of 2DEC on R _i ). For example, the ratio of the potential energy to the kinetic energy (q ²/2), approaches −R _i /(1−R _i ), the ratio of turbulent energy production by buoyancy forces to production by shearing forces (i.e. the flux Richardson number, R _if ), approaches R _i . Also, the Richardson number dependence of the principal angle (β) of the Reynolds stress tensor and the angle (β_ρ) of the scalar flux vector is fairly predicted by this model .On the other hand, it is found that the above ratios are insensitive to viscosity, while the ratios ɛ /q ² and , depend on the viscosity and they evolve asymptotically like t ⁻¹. The turbulent Froude number, F _rt =(L _Oz /L _E )^2/3, where L _Oz and L _E are the Ozmidov length scale and the Ellison length scale, respectively, evolves asymptotically like t ^−1/3.     

Exact Connections of Warping Fields and Torsional Rigidities of Symmetric Composite Bars

Saint-Venant's torsion of symmetric cylindrical bars consisting of two or four homogeneous phases is studied. A symmetric section is meant that the cross section of the cylindrical bar possesses reflectional symmetry with respect to one or more axes. Each constituent region may have different shear modulus. The idea of the analysis is to superimpose suitably reflected potentials to obtain the torsion solution of the same composite section but with different moduli. For two-phase sections, we show that, if the warping fields for a given symmetric section with phase shear moduli μ₁ and μ₂ are known a priori, then the warping fields for the same configuration but with a different set of constituent moduli μ₁ ^′ and μ₂ ^′ are readily found through simple linear superpositions. Further, suppose that the torsional rigidities T(μ₁,μ₂) and T(μ₁ ^′,μ₂ ^′) for any two sets of phase moduli can be measured by some experimental tests or evaluated through numerical procedures, then the torsional rigidity for any other combinations of constituent moduli T(μ₁ ^′′,μ₂ ^′′) can be exactly determined without any recourse to the field solutions of governing differential equations. Similar procedures can be applied to a 4-phase symmetric section. But the coefficients of superposition are only found for a few branches. Specifically, we find that depending on the conditions of μ and μ^′, admissible solutions can be divided into three categories. When the correspondence between the warping field is known to exist, a link between the torsional rigidities can be established as well. This revised version was published online in June 2006 with corrections to the Cover Date.     

Entire Solutions with Merging Fronts to Reaction–Diffusion Equations

We deal with a reaction–diffusion equation u _t = u _xx + f(u) which has two stable constant equilibria, u = 0, 1 and a monotone increasing traveling front solution u = φ(x + ct) (c > 0) connecting those equilibria. Suppose that u = a (0 < a < 1) is an unstable equilibrium and that the equation allows monotone increasing traveling front solutions u = ψ₁(x + c ₁ t) (c ₁ < 0) and ψ₂(x + c ₂ t) (c ₂ > 0) connecting u = 0 with u = a and u = a with u = 1, respectively. We call by an entire solution a classical solution which is defined for all . We prove that there exists an entire solution such that for t≈ − ∞ it behaves as two fronts ψ₁(x + c ₁ t) and ψ₂(x + c ₂ t) on the left and right x-axes, respectively, while it converges to φ(x + ct) as t→∞. In addition, if c > − c ₁, we show the existence of an entire solution which behaves as ψ₁( − x + c ₁ t) in and φ(x + ct) in for t≈ − ∞.     

Similarities Between Rotation and Stratification Effects on Homogeneous Shear Flow

Linear theory is applied to examine rotation and buoyancy effects on homogeneous turbulent shear flows with given vertical velocity shear, S=d/dx ₃. In the rotating shear case (where the rotation vector is perpendicular to the plane of the mean flow, Ω _i =Ωδ _{i 2}), general solutions for the Fourier components of the fluctuating velocity are proposed. These solutions are compared with those proposed in the literature for the Fourier components of the fluctuating velocity and density in the case of a homogeneous stratified shear flow with vertical density gradient, S _ρ=d/dx ₃. It is shown that from the normal mode stability stand point the Bradshaw parameter B=2Ω/S(1+2Ω/S) (in the rotating shear case) and the Richardson number R _i (in the statified shear case) play similar roles in identifying the stability for all the wave components except in the case where Ω·k=0, for which rotation has no effects on the flow. Analysis of the long-time behavior of the non-dimensional spectral density of energy, e ^g , is carried out. In the stable case, e ^g has decaying oscillations or undergoes a power law decay in time. Analytical solutions for the streamwise two-dimensional energy ℰ _ii ¹/2 (i.e. the limit at k ₁=0 of the one-dimensional energy spectra) are proposed. At large time, ℰ _ii ¹(t)/ℰ _ii ¹(0) oscillates around the value (3R _i +1)/(4R _i ) except at R _i =1 it stays constant in time. Similar behavior for ℰ _ii ¹(t)/ℰ _ii ¹(0) is also observed in the rotating shear case (ℰ _ii ¹(t)/ℰ _ii ¹(0) oscillates around the value (1+4B)/(4B)). Due to the behavior of the dimensionless spectral density of energy in both flow cases, the turbulent kinetic energy, /2, the production rate, ?, and the rate due to the buoyancy forces, ℬ, are split into two parts, , ?=?₁+?₂, ℬ=ℬ₁+ℬ₂ (in the stratified shear case, both ?₁ and ℬ₁ vanish when R _i >?, while in the rotating shear case one has ℬ=0). It is shown that when rotation is “cyclonic” (i.e. Ω/S>0), part reaches maximum magnitudes at St ≈2, independent of the B value, and the first time to a zero crossing of ?₂ occurs at this particular value. When rotation is “anticyclonic” (i.e. Ω/S<0) one finds St ≈1.6 instead of St ≈2. In the stratified shear case, both ?₂ and ℬ₂ cross zero at Nt=St ≈2, and part reaches maximum magnitudes at this particular value. These results and in particular those for the turbulent kinetic energy are compared with previous direct numerical simulation (DNS) results in homogeneous stratified shear flows. Received 30 July 2001 and accepted 19 February 2002     

Dirichlet problem with Φ-Laplacian and mixed singularities

Assume that A₁, A₂ ⊂ ℝ are closed intervals containing 0, ϕ is an increasing odd homeomorphism with ϕ (ℝ) = ℝ, and T ∈ (0, ∞). We study a singular Dirichlet problem of the form {fx080-01} and prove the existence of its smooth solution satisfying the conditions {fx080-02}, where ƒ satisfies the Carathéodory conditions on the set (0, T) × D and can have time singularities at t = 0 and t = T and space singularities at x = 0, y = 0. Published in Neliniini Kolyvannya, Vol. 11, No. 1, pp. 81–95, January–March, 2007.     

Singular periodic impulse problems

We establish an existence principle for the impulsive periodic boundary-value problem {fx029-01}, where g ∈ C(0, ∞) can have a strong singularity at the origin. Furthermore, we assume that 0 < t ₁ < … < t _m < T, e ∈ L ₁[0, T], c ∈ ℝ, J _i and M _i , i = 1, 2, …, m, are continuous mappings of G[0, T] × G[0, T] into ℝ, and G[0, T] denotes the space of functions regulated on [0, T]. The presented principle is based on an averaging procedure similar to that introduced by Manásevich and Mawhin for singular periodic problems with p-Laplacian. Published in Neliniini Kolyvannya, Vol. 11, No. 1, pp. 32–44, January–March, 2007.     

On the use of crack-tip-opening displacement to predict the fracture strength of notched graphite/epoxy laminates

The fracture strength and crack-opening displacement of notched graphite/epoxy laminates were measured experimentally using the center-cracked tension-specimen geometry. Four replicate tests were conducted for a variety of laminate stacking sequences, thicknesses, and notch lengths. Most laminates exhibited extensive notch-tip damage prior to fracture. Values of crack-tip-opening displacement (CTOD) at fracture were estimated from values of crack-opening displacement measured at the crack center line. CTOD was independent of specimen crack length for the [0/±45/90] _s , [0/±45/90]_15s , [0/±45] _s , [0/±45/]_15s , and [0/90]_24s laminates. In addition, notched laminate strength was accurately predicted using a Dugdale-type model along with the estimated CTOD.Paper was presented at V International Congress on Experimental Mechanics held in Montreal, Quebec, Canada on June 10–15, 1984.     

Influence of molecular parameters on the stress dependence of viscous and elastic properties of polypropylene melts in shear

The stress dependencies of the steady-state viscosity η and, particularly, that of the steady-state elastic compliance J _e of various linear isotactic polypropylenes (PP) and one long-chain branched PP are investigated using creep-recovery tests. The creep stresses applied range from 2 to 10,000 Pa. In order to discuss the stress-dependent viscosity η and elastic compliance J _e with respect to the influence of the weight average molar mass M _w and the polydispersity factor M _w/M _n the PP are characterized by SEC–MALLS. For the linear PP, linear steady-state elastic compliances J_e⁰J_{\rm e}^0 in the range of 10^− 5–10^− 3 Pa^− 1 are obtained depending on the molar mass distribution. J_e⁰J_{\rm e}^0 of the LCB-PP is distinctly higher and comes to lie at around 10^− 2 Pa^− 1. J_e⁰J_{\rm e}^0 is found to be independent of M _w but strongly dependent on polydispersity. η and J _e decrease with increasing stress. For the linear PP, J _e as a function of the stress τ is temperature independent. The higher M _w/M _n the stronger is the shear thinning of η and the more pronounced is the stress dependence of J _e. For the LCB-PP, the strongest stress dependence of η and J _e is observed. Furthermore, for all PP J _e reacts more sensitively to an increasing stress than η. A qualitative explanation for the stronger stress dependence of J _e compared to η is given by analyzing the contribution of long relaxation times to the viscosity and elasticity.     

Streamwise evolution of an inclined cylinder wake

The streamwise evolution of an inclined circular cylinder wake was investigated by measuring all three velocity and vorticity components using an eight-hotwire vorticity probe in a wind tunnel at a Reynolds number Re_d of 7,200 based on free stream velocity (U _∞) and cylinder diameter (d). The measurements were conducted at four different inclination angles (α), namely 0°, 15°, 30°, and 45° and at three downstream locations, i.e., x/d = 10, 20, and 40 from the cylinder. At x/d = 10, the effects of α on the three coherent vorticity components are negligibly small for α ≤ 15°. When α increases further to 45°, the maximum of coherent spanwise vorticity reduces by about 50%, while that of the streamwise vorticity increases by about 70%. Similar results are found at x/d = 20, indicating the impaired spanwise vortices and the enhancement of the three-dimensionality of the wake with increasing α. The streamwise decay rate of the coherent spanwise vorticity is smaller for a larger α. This is because the streamwise spacing between the spanwise vortices is bigger for a larger α, resulting in a weak interaction between the vortices and hence slower decaying rate in the streamwise direction. For all tested α, the coherent contribution to [`(v<sup>2</sup>)] \overline{{v^{2}}} is remarkable at x/d = 10 and 20 and significantly larger than that to [`(u²)] \overline{{u^{2}}} and [`(w<sup>2</sup>)]. \overline{{w^{2}}}. This contribution to all three Reynolds normal stresses becomes negligibly small at x/d = 40. The coherent contribution to [`(u²)] \overline{{u^{2}}} and [`(v²)] \overline{{v^{2}}} decays slower as moving downstream for a larger α, consistent with the slow decay of the coherent spanwise vorticity for a larger α.     

A unified approach to closed-form solutions of moving heat-source problems

A closed-form model for the computation of temperature distribution in an infinitely extended isotropic body with a time-dependent moving-heat sources is discussed. The temperature solutions are presented for the sources of the forms: (i) ₀₁(t)=₀ exp(−λt), (ii) ₀₂(t) =₀(t/t ^*)exp(−λt), and ₀₃(t)=₀[1+a cos(ωt)], where λ and ω are real parameters and t ^* characterizes the limiting time. The reduced (or dimensionless) temperature solutions are presented in terms of the generalized representation of an incomplete gamma function Γ(α,x;b) and its decomposition C _Γ and S _Γ. The solutions are presented for moving, -point, -line, and -plane heat sources. It is also demonstrated that the present analysis covers the classical temperature solutions of a constant strength source under quasi-steady state situations. Received on 13 June 1997