Mixed convection adjacent to non-isothermal vertical surfaces

Study of laminar boundary layer in mixed convection from vertical plates is carried out. The surface temperature along the vertical plate is assumed to vary arbitrarily with vertical distance. Perturbation technique is used to solve the governing boundary layer equations. The differentials of the wall temperature are used as perturbation elements, which are functions of vertical distance, to obtain universal functions. The universal functions are valid for any type of vertical wall temperature variation. Heat transfer rates and fluid velocity inside the boundary layer can be expressed and calculated using these universal functions. Heat transfer rates are obtained for the special cases of power-law variation of the wall temperature. The effect of the governing parameter (Gr(y)/Re-y(2)) and the power index of the power-law wall temperature variation on heat transfer rates is studied. For the purpose of validation, the mixed convection results obtained by the present technique pertaining to the special cases of isothermal vertical wall are compared with those obtained by similarity analysis reported in literature, and the agreement is found to be good. (C) 2012 Elsevier Ltd. All rights reserved.

Effect of isothermal aging and electromigration on the microstructural evolution of solder interconnections during thermomechanical loading

The effect of different pre-aging treatments on the microstructural evolution of lead-free solder and growth of interfacial intermetallic compound layers under thermal cycling has been investigated in this work. The results show that there are distinct differences in the microstructural changes between samples with no pretreatment, samples that have experienced thermal annealing at 125A degrees C for 750 h before thermal cycling, and those that have had direct current (DC) stressing for 750 h as pretreatment. The microstructural evolution of the solder matrix is rationalized by utilizing the science of microstructures and analysis of the influence of electron flow on the precipitation phenomena. The finite-element method is utilized to understand the loading conditions imposed on the solder interconnections during cyclic stressing. The growth of intermetallic reaction layers is further analyzed by utilizing quantitative thermodynamic calculations coupled with kinetic analysis. The latter is based on the changes in the intrinsic diffusion fluxes of elements induced by current flow and alloying elements present in the system. With this concurrent approach the differences seen in thermal cycling behavior between the different pre-aging treatments can be explained.

Isothermal transformation of beta-phase in Cu-rich Cu-Al-Sn alloys

Microstructural changes resulting from isothermal decomposition of the beta-phase have been studied in Cu-rich binary Cu-Al and ternary Cu-Al-Sn alloys containing up to 3 at.% Sn at temperatures from 873 to 673 K. Results are presented as TTT diagrams. The decomposition occurs in several stages, each of which involves the establishment of metastable equilibrium between beta and one or more of the product phases alpha, beta(1) and gamma(2). Addition of Sn has been shown to increase the stability of the ordered beta(1)-phase in relation to beta. In alloys containing more than 2 at.% Sn, the beta(1) emerges as a stable phase. At low Sn concentrations beta(1) is metastable. An important new finding is the existence of three-phase equilibrium microstructure containing alpha, beta(1) and gamma(2). Increasing addition of Sn alters the morphology of beta(1) from rosettes to dendrites and finally to Widmanstatten needles.

Transition and acoustic response of recirculation structures in an unconfined co-axial isothermal swirling flow

This paper reports the first observations of transition from a pre-vortex breakdown (Pre-VB) flowreversal to a fully developed central toroidal recirculation zone in a non-reacting, double-concentric swirling jet configuration and its response to longitudinal acoustic excitation. This transition proceeds with the formation of two intermediate, critical flow regimes. First, a partially penetrated vortex breakdown bubble (VBB) is formed that indicates the first occurrence of an enclosed structure as the centre jet penetration is suppressed by the growing outer roll-up eddy; resulting in an opposed flow stagnation region. Second, a metastable transition structure is formed that marks the collapse of inner mixing vortices. In this study, the time-averaged topological changes in the coherent recirculation structures are discussed based on the non-dimensional modified Rossby number (Ro(m)) which appears to describe the spreading of the zone of swirl influence in different flow regimes. Further, the time-mean global acoustic response of pre-VB and VBB is measured as a function of pulsing frequency using the relative aerodynamic blockage factor (i.e., maximum radial width of the inner recirculation zone). It is observed that all flow modes except VBB are structurally unstable as they exhibit severe transverse radial shrinkage (similar to 20%) at the burner Helmholtz resonant modes (100-110 Hz). In contrast, all flow regimes show positional instability as seen by the large-scale, asymmetric spatial shifting of the vortex core centres. Finally, the mixing transfer function M (f) and magnitude squared coherence lambda(2)(f) analysis is presented to determine the natural couplingmodes of the system dynamic parameters (u', p'), i.e., local acoustic response. It is seen that the pre-VB flow mode exhibits a narrow-band, low pass filter behavior with a linear response window of 100-105 Hz. However, in the VBB structure, presence of critical regions such as the opposed flow stagnation region alters the linearity range with the structure showing a response even at higher pulsing frequencies (100-300 Hz). (C) 2013 AIP Publishing LLC.

Transition in vortex breakdown modes in a coaxial isothermal unconfined swirling jet

This paper reports first observations of transition in recirculation pattern from an open-bubble type axisymmetric vortex breakdown to partially open bubble mode through an intermediate, critical regime of conical sheet formation in an unconfined, co-axial isothermal swirling flow. This time-mean transition is studied for two distinct flow modes which are characterized based on the modified Rossby number (Ro(m)), i.e., Ro(m) <= 1 and Ro(m) > 1. Flow modes with Ro(m) <= 1 are observed to first undergo cone-type breakdown and then to partially open bubble state as the geometric swirl number (S-G) is increased by similar to 20% and similar to 40%, respectively, from the baseline open-bubble state. However, the flow modes with Ro(m) > 1 fail to undergo such sequential transition. This distinct behavior is explained based on the physical significance associated with Ro(m) and the swirl momentum factor (xi). In essence, xi represents the ratio of angular momentum distributed across the flow structure to that distributed from central axis to the edge of the vortex core. It is observed that xi increases by similar to 100% in the critical swirl number band where conical breakdown occurs as compared to its magnitude in the S-G regime where open bubble state is seen. This results from the fact that flow modes with Ro(m) <= 1 are dominated by radial pressure gradient due to swirl/rotational effect when compared to radial pressure deficit arising from entrainment (due to the presence of co-stream). Consequently, the imparted swirl tends to penetrate easily towards the central axis causing it to spread laterally and finally undergo conical sheet breakdown. However, the flow modes with Ro(m) > 1 are dominated by pressure deficit due to entrainment effect. This blocks the radial inward penetration of imparted angular momentum thus preventing the lateral spread of these flow modes. As such these structures fail to undergo cone mode of vortex breakdown which is substantiated by a mere 30%-40% rise in xi in the critical swirl number range. (C) 2014 AIP Publishing LLC.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

Formation of Prestellar Cores via Non-Isothermal Gas Fragmentation

Sheet-like clouds are common in turbulent gas and perhaps form via collisions between turbulent gas flows. Having examined the evolution of an isothermal shocked slab in an earlier contribution, in this work we follow the evolution of a sheet-like cloud confined by (thermal) pressure and gas in it is allowed to cool. The extant purpose of this endeavour is to study the early phases of core-formation. The observed evolution of this cloud supports the conjecture that molecular clouds themselves are three-phase media (comprising viz. a stable cold and warm medium, and a third thermally unstable medium), though it appears, clouds may evolve in this manner irrespective of whether they are gravitationally bound. We report, this sheet fragments initially due to the growth of the thermal instability (TI) and some fragments are elongated, filament-like. Subsequently, relatively large fragments become gravitationally unstable and sub-fragment into smaller cores. The formation of cores appears to be a three stage process: first, growth of the TI leads to rapid fragmentation of the slab; second, relatively small fragments acquire mass via gas-accretion and/or merger and third, sufficiently massive fragments become susceptible to the gravitational instability and sub-fragment to form smaller cores. We investigate typical properties of clumps (and smaller cores) resulting from this fragmentation process. Findings of this work support the suggestion that the weak velocity field usually observed in dense clumps and smaller cores is likely seeded by the growth of dynamic instabilities. Simulations were performed using the smooth particle hydrodynamics algorithm.

Acoustic response of vortex breakdown modes in a coaxial isothermal unconfined swirling jet

The present experimental work is concerned with the study of amplitude dependent acoustic response of an isothermal coaxial swirling jet. The excitation amplitude is increased in five distinct steps at the burner's Helmholtz resonator mode (i.e., 100 Hz). Two flow states are compared, namely, sub-critical and super-critical vortex breakdown (VB) that occur before and after the critical conical sheet breakdown, respectively. The geometric swirl number is varied in the range 2.14-4.03. Under the influence of external pulsing, global response characteristics are studied based on the topological changes observed in time-averaged 2D flow field. These are obtained from high resolution 2D PIV (particle image velocimetry) in the longitudinal-mid plane. PIV results also illustrate the changes in the normalized vortex core coordinates (r(vcc)/(r(vcc))(0) (Hz), y(vcc)/(y(vcc))(0) (Hz)) of internal recirculation zone (IRZ). A strong forced response is observed at 100 Hz (excitation frequency) in the convectively unstable region which get amplified based on the magnitude of external forcing. The radial extent of this forced response region at a given excitation amplitude is represented by the acoustic response region (b). The topological placement of the responsive convectively unstable region is a function of both the intensity of imparted swirl (characterized by swirl number) and forcing amplitude. It is observed that for sub-critical VB mode, an increase in the excitation amplitude till a critical value shifts the vortex core centre (particularly, the vortex core moves downstream and radially outwards) leading to drastic fanning-out/widening of the IRZ. This is accompanied by similar to 30% reduction in the recirculation velocity of the IRZ. It is also observed that b < R (R: radial distance from central axis to outer shear layer-OSL). At super-critical amplitudes, the sub-critical IRZ topology transits back (the vortex core retracts upstream and radially inwards) and finally undergoes a transverse shrinkage ((r(vcc))/(r(vcc))(0 Hz) decreases by similar to 20%) when b >= R. In contrast, the vortex core of super-critical breakdown mode consistently spreads radially outwards and is displaced further downstream. Finally, the IRZ fans-out at the threshold excitation amplitude. However, the acoustic response region b is still less than R. This is explained based on the characteristic geometric swirl number (S-G) of the flow regimes. The super-critical flow mode with higher S-G (hence, higher radial pressure drop due to rotational effect which scales as Delta P similar to rho u theta(2) and acts inwards towards the center line) compared to sub-critical state imposes a greater resistance to the radial outward spread of b. As a result, the acoustic energy supplied to the super-critical flow mode increases the degree of acoustic response at the pulsing frequency and energizes its harmonics (evident from power spectra). As a disturbance amplifier, the stronger convective instability mode within the flow structure of super-critical VB causes the topology to widen/fan-out severely at threshold excitation amplitude. (C) 2015 AIP Publishing LLC.

A novel interferometric technique to estimate thermal diffusivity of `optically transparent solid' using `Isothermal surface velocimetry'

One-dimensional transient heat flow is interpreted as a procession of `macro-scale translatory motion of indexed isothermal surfaces'. A new analytical model is proposed by introducing velocity of isothermal surface in Fourier heat diffusion equation. The velocity dependent function is extracted by revisiting `the concept of thermal layer of heat conduction in solid' and `exact solution' to estimate thermal diffusivity. The experimental approach involves establishment of 1 D unsteady heat flow inside the sample through Step-temperature excitation. A novel self-reference interferometer is utilized to separate a `unique isothermal surface' in time-varying temperature field. The translatory motion of the said isothermal surface is recorded using digital camera to estimate its velocity. From the knowledge of thermo-optic coefficient, temperature of the said isothermal surface is predicted. The performance of proposed method is evaluated for Quartz sample and compared with literature.

Modeling of the non-isothermal liquid droplet impact on a heated solid substrate with heterogeneous wettability

A comprehensive numerical investigation on the impingement and spreading of a non-isothermal liquid droplet on a solid substrate with heterogeneous wettability is presented in this work. The time-dependent incompressible Navier-Stokes equations are used to describe the fluid flow in the liquid droplet, whereas the heat transfer in the moving droplet and in the solid substrate is described by the energy equation. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the time-dependent incompressible Navier-Stokes equation and the energy equation in the time-dependent moving domain. Moreover, the Marangoni convection is included in the variational form of the Navier-Stokes equations without calculating the partial derivatives of the temperature on the free surface. The heterogeneous wettability is incorporated into the numerical model by defining a space-dependent contact angle. An array of simulations for droplet impingement on a heated solid substrate with circular patterned heterogeneous wettability are presented. The numerical study includes the influence of wettability contrast, pattern diameter, Reynolds number and Weber number on the confinement of the spreading droplet within the inner region, which is more wettable than the outer region. Also, the influence of these parameters on the total heat transfer from the solid substrate to the liquid droplet is examined. We observe that the equilibrium position depends on the wettability contrast and the diameter of the inner surface. Consequently. the heat transfer is more when the wettability contrast is small and/or the diameter of inner region is large. The influence of the Weber number on the total heat transfer is more compared to the Reynolds number, and the total heat transfer increases when the Weber number increases. (C) 2015 Elsevier Ltd. All rights reserved.

Stability of ZrTiCuNiBe Bulk Metallic Glass Upon Isothermal Annealing Near the Glass Transition Temperature

The stability of Zr41Ti14Cu12.5Ni10Be22.5 bulk metallic glass (BMG) upon isothermal annealing near the glass transition temperature has been investigated by using x-ray diffraction, differential scanning calorimetry, and the pulse echo overlap method. The density, elastic constants, and thermodynamic parameters as well as their annealing time dependence have been determined. The microstructural and properties changes of the annealed BMG were checked by acoustic measurement. Obvious structural and property changes were observed with prolonged annealing of the BMG near the glass transition temperature.