Thermal reliability testing of functionally gradient materials using thermal shock method

We found a solution of an unsteady two-dimensional heat conduction equation in a functionally gradient material (FGM) which is subjected to a double thermal shock, namely, a local heating of a specimen by a power laser beam and cooling of a heated surface by a water-air spray. We developed an analytical method whereby a coating is described as a laminated plate composed of n layers with the constant material properties within a layer. Temperature distribution in a nonhomogeneous laminated plate is obtained in a form of series using the Laplace–Hankel integral transforms. In order to extend the model of a laminated plate to describe FGM where thermal physical characteristics are continuous functions of spatial coordinate, we considered the limiting case of the obtained temperature distribution when the thickness of the layer iΔ _i → 0, and the number of layers n→∞. This allowed us to obtain the temperature distribution in an easy-to-use analytical form which can be used for determining thermal stresses in FGM. The dependence of the temperature distribution in FGM on the operating parameters of a double thermal shock method, e.g., a duration of heating, laser beam radius, the rate of a spray cooling, is discussed. Received on 3 May 1999     

Drag of a strongly heated sphere at small Reynolds numbers

Using an asymptotic small-perturbation method, the flow around a strongly heated sphere at small Reynolds numbers Re ≪ 1 is considered with account for thermal stresses in the gas in the higher-order approximations, beyond the Stokes one. It is assumed that the value of the Prandtl number Pr is arbitrary and the temperature dependence of the viscosity is described by a power law with an arbitrary exponent. In the O(Re²) and O(Re³ ln(Re)) approximations, the drag force of a heated sphere is found over a wide range of the ratios of sphere’s temperature to the gas free-stream temperature T _W /T _∞. The limits of applicability of the first (in Re) approximation are investigated, including the negative-drag effect, attributable to the action of the thermal stresses. The results are compared with numerical calculations of the flow around a hot sphere. The limits of applicability of the approximations found are examined. Similar results are obtained for the standard Navier-Stokes equations in which the thermal stresses are neglected.     

Analytical Solution for Radial Deformations of Functionally Graded Isotropic and Incompressible Second-Order Elastic Hollow Spheres

We analytically analyze radial expansion/contraction of a hollow sphere composed of a second-order elastic, isotropic, incompressible and inhomogeneous material to delineate differences and similarities between solutions of the first- and the second-order problems. The two elastic moduli are assumed to be either affine or power-law functions of the radial coordinate R in the undeformed reference configuration. For the affine variation of the shear modulus μ, the hoop stress for the linear elastic (or the first-order) problem at the point R=(R _ou R _in (R _ou +R _in )/2)^1/3 is independent of the slope of the μ vs. R line. Here R _in and R _ou equal, respectively, the inner and the outer radius of the sphere in the reference configuration. For μ(R)∝R ⁿ , for the linear problem, the hoop stress is constant in the sphere for n=1. However, no such results are found for the second-order (i.e., materially nonlinear) problem. Whereas for the first-order problem the shear modulus influences only the radial displacement and not the stresses, for the second-order problem the two elastic constants affect both the radial displacement and the stresses. In a very thick homogeneous hollow sphere subjected only to pressure on the outer surface, the hoop stress at a point on the inner surface depends upon values of the two elastic moduli. Thus conclusions drawn from the analysis of the first-order problem do not hold for the second-order problem. Closed form solutions for the displacement and stresses for the first-order and the second-order problems provided herein can be used to verify solutions of the problem obtained by using numerical methods.     

Further Comments on “Combined Forced and Free Convective Flow in a Vertical Porous Channel: The Effects of Viscous Dissipation and Pressure Work”

This note is concerned with the assertion of Barletta and Nield (2009a) that “a fluid with a thermal expansion coefficient greater than that of a perfect gas (β > β _perfect gas) is of marginal or no interest in the framework of convection in porous media”, and that for a remark of Magyari (Transp. Porous Media, 2009) about the forced convection eigenflow solutions, the circumstance β > β _perfect gas does not represent “a sound physical basis”. Here, it is shown, however, that these assertions are in contradiction with the experimentally measured values of β for important technical fluids as e.g., air, nitrogen, carbon dioxide, and ammonia where, in the temperature range between −20 and +100°C, just the inequality β > β _perfect gas holds.     

Flow of power-law fluids over a moving wedge surface with wall mass injection

The boundary layer problem of a power-law fluid flow with fluid injection on a wedge whose surface is moving with a constant velocity in the opposite direction to that of the uniform mainstream is analyzed. The free stream velocity, the injection velocity at the surface, moving velocity of the wedge surface, the wedge angle and the power law index of non-Newtonian fluid are assumed variables. The fourth order Runge–Kutta method modified by Gill is used to solve the non-dimensional boundary layer equations for non-Newtonian flow field. Without fluid injection, for every angle of wedge β, a limiting value for velocity ratio λ _cr (velocity of the wedge surface/velocity of the uniform flow) is found for each power-law index n. The value of λ _cr increases with the increasing wedge angle β. The value of wedge angle also restricts the physical characteristics of the fluid to be used. The effects of the different parameters on velocity profile and on skin friction are studied and the drag reduction is discussed. In case of C = 2.5 and velocity ratio λ = 0.2 for wedge angle β = 0.5 with the fluid with power law-index n = 0.5, 48.8% drag reduction is obtained.     

Equilibrium states in homogeneous turbulence with baroclinic instability

The evolution of energies and fluxes in homogeneous turbulence with baroclinic instability is analyzed using the linear theory. The mean flow corresponds to a vertical shear having a uniform mean velocity gradient, ?U _i /?x _j  = S δ _i1 δ _j3, a system rotation about the vertical axis with rate Ω, Ω _i  = Ωδ _i3, and uniform buoyancy gradients in the spanwise ${(\partial B{/}\partial x_2\,{=}\, N_h^2\,{=}\,-2\Omega S)}The evolution of energies and fluxes in homogeneous turbulence with baroclinic instability is analyzed using the linear theory. The mean flow corresponds to a vertical shear having a uniform mean velocity gradient, ∂U _i /∂x _j  = S δ _i1 δ _j3, a system rotation about the vertical axis with rate Ω, Ω _i  = Ωδ _i3, and uniform buoyancy gradients in the spanwise (?B/?x₂ = N_h² = -2WS){(\partial B{/}\partial x_2\,{=}\, N_h^2\,{=}\,-2\Omega S)} and vertical (?B/?x₃ = N_v²){(\partial B{/}\partial x_3\,{=}\,N_v^2)} directions. Computations based on the rapid distortion theory (RDT) are performed for several values of the rotation number R = 2Ω/S and the Richardson number R_i = N_v²/S² < 1{R_i\,{=}\,N_v^2/S^2 <1 }. It is shown that, during an initial phase, the energies and the buoyancy fluxes are sensitive to the effects of pressure and viscosity. At large time, the ratios of energies, as well as the normalized fluxes, evolve to an asymptotically constant value, while the pressure–strain correlation scaled with the product of the turbulent kinetic energy by the shear rate approaches zero. Accordingly, an analytical parametric study based on the “pressure-less” approach (PLA) is also presented. The analytical study indicates that, when R _i  < 1, there is an exponential instability and equilibrium states of turbulence, in agreement with RDT. The energies and the buoyancy fluxes grow exponentially for large times with the same rate (γ in St units). The asymptotic value of the ratios of energies yielded by RDT is well described by its PLA counterpart derived analytically. At R _i  = 0, the asymptotic value of γ increases with increasing R approaching 2 for high rotation rates. At low rotation rates, an important contribution to the kinetic energy comes from the streamwise kinetic energy, whereas, at high rotation rates, the contribution of the vertical kinetic energy is dominant. When 0 < R _i  < 1 and R 1 0{R\ne 0}, the asymptotic value of γ decreases as R _i increases so as it becomes zero at R _i  = 1.     

Non-isobaric Marangoni boundary layer flow for Cu,Al<Subscript>2</Subscript>O<Subscript>3</Subscript> and TiO<Subscript>2</Subscript> nanoparticles in a water based fluid

In this paper, a non-isobaric Marangoni boundary layer flow that can be formed along the interface of immiscible nanofluids in surface driven flows due to an imposed temperature gradient, is considered. The solution is determined using a similarity solution for both the momentum and energy equations and assuming developing boundary layer flow along the interface of the immiscible nanofluids. The resulting system of nonlinear ordinary differential equations is solved numerically using the shooting method along with the Runge-Kutta-Fehlberg method. Numerical results are obtained for the interface velocity, the surface temperature gradient as well as the velocity and temperature profiles for some values of the governing parameters, namely the nanoparticle volume fraction φ (0≤φ≤0.2) and the constant exponent β. Three different types of nanoparticles, namely Cu, Al₂O₃ and TiO₂ are considered by using water-based fluid with Prandtl number Pr =6.2. It was found that nanoparticles with low thermal conductivity, TiO₂, have better enhancement on heat transfer compared to Al₂O₃ and Cu. The results also indicate that dual solutions exist when β<0.5. The paper complements also the work by Golia and Viviani (Meccanica 21:200–204, 1986) concerning the dual solutions in the case of adverse pressure gradient.     

The onset of buoyancy-driven convection in a ferromagnetic fluid saturated porous medium

The onset of buoyancy-driven convection in an initially quiescent ferrofluid saturated horizontal porous layer in the presence of a uniform vertical magnetic field is investigated. The Brinkman-Lapwood extended Darcy equation with fluid viscosity different from effective viscosity is used to describe the flow in the porous medium. The lower boundary of the porous layer is assumed to be rigid-paramagnetic, while the upper paramagnetic boundary is considered to be either rigid or stress-free. The thermal conditions include fixed heat flux at the lower boundary, and a general convective–radiative exchange at the upper boundary, which encompasses fixed temperature and fixed heat flux as particular cases. The resulting eigenvalue problem is solved numerically using the Galerkin technique. It is found that increase in the Biot number Bi, porous parameter σ, viscosity ratio Λ, magnetic susceptibility χ, and decrease in the magnetic number M ₁ and non-linearity of magnetization M ₃ is to delay the onset of ferroconvection in a porous medium. Further, increase in M ₁, M ₃, and decrease in χ, Λ, σ and Bi is to decrease the size of convection cells.     

Miscible Thermo-Viscous Fingering Instability in Porous Media. Part 2: Numerical Simulations

The development of the thermo-viscous fingering instability of miscible displacements in homogeneous porous media is examined. In this first part of the study dealing with stability analysis, the basic equations and the parameters governing the problem in a rectilinear geometry are developed. An exponential dependence of viscosity on temperature and concentration is represented by two parameters, thermal mobility ratio β _T and a solutal mobility ratio β _C , respectively. Other parameters involved are the Lewis number Le and a thermal-lag coefficient λ. The governing equations are linearized and solved to obtain instability characteristics using either a quasi-steady-state approximation (QSSA) or initial value calculations (IVC). Exact analytical solutions are also obtained for very weakly diffusing systems. Using the QSSA approach, it was found that an increase in thermal mobility ratio β _T is seen to enhance the instability for fixed β _C , Le and λ. For fixed β _C and β _T , a decrease in the thermal-lag coefficient and/or an increase in the Lewis number always decrease the instability. Moreover, strong thermal diffusion at large Le as well as enhanced redistribution of heat between the solid and fluid phases at small λ is seen to alleviate the destabilizing effects of positive β _T . Consequently, the instability gets strictly dominated by the solutal front. The linear stability analysis using IVC approach leads to conclusions similar to the QSSA approach except for the case of large Le and unity λ flow where the instability is seen to get even less pronounced than in the case of a reference isothermal flow of the same β _C , but β _T  = 0. At practically, small value of λ, however, the instability ultimately approaches that due to β _C only.     

Motion of a sphere in an electrically conducting rotating fluid

The combined effect of rotation and magnetic field is investigated for the axisymmetric flow due to the motion of a sphere in an inviscid, incompressible electrically conducting fluid having uniform rotation far upstream. The steady-state linearized equations contain a single parameter α=1/2βR _m, β being the magnetic pressure number and R _m the magnetic Reynolds number. The complete solution for the flow field and magnetic field is obtained and the distribution of vorticity and current density is found. The induced vorticity is O(α⁴) and the current density is O(R _m) on the sphere.     

Slip-flow heat transfer in a microchannel with viscous dissipation

Forced convection flow in a microchannel with constant wall temperature is studied, including viscous dissipation effect. The slip-flow regime is considered by incorporating both the velocity-slip and the temperature-jump conditions at the surface. The energy equation is solved for the developing temperature field using finite integral transform. To increase β_v Kn is to increase the slip velocity at the wall surface, and hence to decrease the friction factor. Effects of the parameters β_v Kn, β, and Br on the heat transfer results are illustrated and discussed in detail. For a fixed Br, the Nusselt number may be either higher or lower than those of the continuum regime, depending on the competition between the effects of β_v Kn and β. At a given β_v Kn the variation of local Nusselt number becomes more even when β becomes larger, accompanied by a shorter thermal entrance length. The fully developed Nusselt number decreases with increasing β irrelevant to β_v Kn. The increase in Nusselt number due to viscous heating is found to be more pronounced at small β_v Kn.     

Numerical study of transient laminar natural convection heat transfer over a sphere subjected to a constant heat flux

Transient laminar natural convection over a sphere which is subjected to a constant heat flux has been studied numerically for high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and a wide range of Prandtl numbers (Pr = 0.02, 0.7, 7, and 100). A plume with a mushroom-shaped cap forms above the sphere and drifts upward continuously with time. The size and the level of temperature of the transient cap and plume stem decrease with increasing Gr and Pr. Flow separation and an associated vortex may appear in the wake of the sphere depending on the magnitude of Gr and Pr. A recirculation vortex which appears and grows until “steady state” is attained was found only for the very high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and the lowest Prandtl number considered (Pr = 0.02). The appearance and subsequent disappearance of a vortex was observed for Gr = 10⁹ and Pr = 0.7. Over the lower hemisphere, the thickness of both the hydrodynamic (δ_H) and the thermal (δ_T) boundary layers remain nearly constant and the sphere surface is nearly isothermal. The surface temperature presents a local maximum in the wake of the sphere whenever a vortex is established in the wake of the sphere. The surface pressure recovery in the wake of the sphere increases with decreasing Pr and with increasing Gr. For very small Pr, unlike forced convection, the ratio δ_T/δ_H remains close to unity. The results are in good agreement with experimental data and in excellent agreement with numerical results available in the literature. A correlation has also been presented for the overall Nusselt number as a function of Gr and Pr.     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

Time-fractional radial heat conduction in a cylinder and associated thermal stresses

The theory of thermal stresses based on the heat conduction equation with the Caputo time-fractional derivative of order 0 < α ≤ 2 is used to investigate axisymmetic thermal stresses in a cylinder. The solution is obtained applying the Laplace and finite Hankel integral transforms. The Dirichlet and two types of Neumann problems with the prescribed boundary value of the temperature, the normal derivative of the temperature, and the heat flux are considered. Numerical results are illustrated graphically.     

Further Study on Pure Shear

It is known that the Cauchy stress tensor T is a pure shear when trT = 0. An elementary derivation is given for a coordinate system such that, when referred to this coordinate system, the diagonal elements of T vanish while the off-diagonal elements τ ₁, τ ₂, τ ₃, are the pure shears. The structure of τ _i (i = 1, 2, 3) depends on one non-dimensional parameter q = 54(detT)² / [tr(T ²)]³, 0 ≤ q ≤ 1. When q = 0, one of the three τ _i vanishes. A coordinate system can be chosen such that the remaining two have the same magnitude or one of the remaining two also vanishes. When q = 1, all three τ _i have the same magnitude. However, there is a one-parameter family of coordinate systems that gives the same three τ _i . For q ≠ 0 or 1, none of the three τ _i vanishes and the three τ _i in general have different magnitudes. Nevertheless, a coordinate system can be chosen such that two of the three τ _i have the same magnitude. ^★Professor Emeritus of University of Illinois at Chicago and Consulting Professor of Stanford University.     

Fatigue gage utilizing slip-initiation phenomenon in electroplated copper foil

Fundamental experiments are carried out to examine the parameter that dominates the slip-band initiation in electroplated copper foil under the condition where the mean stress as well as the stress amplitude varies. In the case of constant-amplitude stressing, the relation between the critical stress for the slip-band initiation σ _p and the number of cyclesN is well represented by σ _p ^α N=C. In other words, the slip bands appear when the total hysteresis energy applied to the copper foil attains a critical value. In the case of variable stresses, the range-pair mainly dominates the occurrence of the slip bands, and Miner's linear cumulative damage rule holds for the accumulation of the fatigue damage for the slip-band initiation. Accordingly, the parameter (Σσ _i ^α n _i/Σn _i)^1/α is equivalent to the critical stress σ _p under constant amplitude stressing, where σ _i andn _i are the stress amplitude and the number of cycles counted by the range-pair method, respectively, and α is the exponent of the σ _p -N relation. Based on these results, the applicability of the copper foil to the fatigue gage that accumulates and indicates a load experience is discussed.     

Localized Thickening of a Compressed Elastic Band

An infinite elastic band is compressed along its unbounded direction, giving rise to a continuous family of homogeneous configurations that is parameterized by the compression rate β < 1 (β = 1 when there is no compression). It is assumed that, for some critical value β ₀, the compression force as a function of β has a strict local extremum and that the linearized equation around the corresponding homogeneous configuration is strongly elliptic. Under these conditions, there are nearby localized deformations that are asymptotically homogeneous. When the compression force reaches a strict local maximum at β ₀, they describe localized thickening and they occur for values of β slightly smaller than β ₀. Since the material is supposed to be hyperelastic, homogeneous and isotropic, the localized deformations are not due to localized imperfections. The method follows the one developed by A. Mielke for an elastic band under traction: interpretation of the nonlinear elliptic system as an infinite dimensional dynamical system in which the unbounded direction plays the role of time, its reduction to a center manifold and the existence of a homoclinic solution to the reduced finite dimensional problem in [A. Mielke, Hamiltonian and Lagrangian fiows on center manifolds, Lecture Notes in Mathematics 1489. Springer, Berlin Heidelberg New York, 1991]. The main difference lies in the fact that Agmon's condition does not hold anymore and therefore the linearized problem cannot be analyzed as in Mielke's work.     

A natural neighbour method based on Fraeijs de Veubeke variational principle for materially non-linear problems

Abstract The natural neighbour method can be considered as one of many variants of the meshless methods. In the present paper, a new approach based on the Fraeijs de Veubeke （FdV） functional, which is initially developed for linear elasticity, is extended to the case of geometrically linear but materially non-linear solids. The new approach provides an original treatment to two classical problems： the numerical evaluation of the integrals over the domain A and the enforcement of boundary conditions of the type ui = hi on Su. In the absence of body forces （Fi = 0）, it will be shown that the calculation of integrals of the type fA .dA can be avoided and that boundary conditions of the type ui = hi on Su can be imposed in the average sense in general and exactly if hi is linear between two contour nodes, which is obviously the case for tTi = O.