Computation of rheological properties of suspensions of rigid rods: stress growth after inception of steady shear flow

The orientation distribution and stress growth for a suspension of rigid rods (or dumbbells) in a Newtonian solvent are calculated for inception of steady shear flow. Galerkin's method, with spherical harmonics as trial functions, is used in the spatial coordinates to obtain a system of ordinary differential equations in time which is solved by the spectral method. The method is applicable over a wide range of dimensionless shear rates (Peclet numbers) and has been coded with standard system-solvent and eigensystem packages. For sufficiently large Peclet numbers, the results give the well known rigid-dumbbell prediction of an overshoot in the shear viscosity and normal stress differences. This overshoot is then followed by an undershoot. An explicit analytical approximation for the fluid stresses is presented which is reasonably accurate for Peclet numbers less than unity.     

The effect of deformation on the effective conductivity of a dilute suspension of drops in the limit of low particle peclet number

The effective thermal conductivity of a dilute suspension of slightly deformed droplets is calculated in the limit of small particle Peclet number for the undisturbed bulk shear, u = γy, and the linear bulk temperature distribution, T = αy. Two distinct cases of small deformation are considered; deformation dominated by interfacial tension forces, and deformation dominated by viscous forces in the drop. The results show that the presence of deformation can cause a fundamental change in form of the dominant, flow-induced contribution to the effective conductivity.     

Investigation of viscosity and thermal conductivity of alumina nanofluids with addition of SDBS

The present study focused on thermal conductivity and viscosity of alumina nanoparticles, at low volume concentrations of 0.01–1.0 % dispersed in the mixture of ethylene glycol and water (mass ratio, 60:40). Sodium dodeobcylbenzene sulfonate (SDBS) was applied for better dispersion and stability of alumina nanoparticles and study of its influence on both thermal conductivity and viscosity. The thermal conductivity established polynomial enhancement pattern with increase of volume concentration up to 0.1 % while linear enhancement was obtained at higher concentrations. In addition, thermal conductivity was enhanced with the rise of temperature. However, the augmentation was negligible compared to that obtained with increase of volume concentration. In contrast, viscosity data showed remarkable reduction with increase of temperature. Meanwhile, viscosity of nanofluids enhanced with loading of alumina nanoparticles. Thermal conductivity and viscosity measurements showed higher values over theoretical predictions. Results showed SDBS at different concentrations has distinct influence on thermal conductivity and viscosity of nanofluid.     

Endochronic analysis for finite elasto-plastic deformation and application to metal tube under torsion and metal rectangular block under biaxial compression

The ordinary differential constitutive equations of endochronic theory are extended to simulate elasto-plastic deformation in the range of finite strain using the concept of corotational rate. Different corotational stress rates (Jaumann, Cotter-Rivlin, Truesdell, Dienes and Mandel) are incorporated into the theory. In addition, a new formulation of the plastic spin, which can be used in the Mandel stress rate, is derived. Theoretical simulations of the axial effects for various materials subjected to simple and pure torsional loading cases are discussed in this study. It is shown that the endochronic theory incorporated with the Mandel stress rate yields the most satisfactory result, as indicated from comparison with the experimental data found in literature.

Finally, theoretical investigation of the deformation subjected to finite proportional and non-proportional biaxial compression is presented. The true relationship between stress and strain can be converted to a nominal stress-strain relationship for biaxial loading through the explicit transformation equations derived in this paper. Experimental data tested by Khan and Wang [1990] (“An Experimental study of Large finite Plastic Deformation in Annealed 1100 Aluminum During Proportional and Non-proportional Biaxial Compression” Int. J. Plasticity, 6, 485) are suitably described by the theory demonstrated from a comparison with the theoretical prediction according to rigid-plastic and elastic-plastic models employed by Huang and Khan [1991]. “An Analysis of Finite Elastic-Plastic Deformation under Biaxial Compression”, Int. J. Plasticity, 7, 219).     

Variation of fiber orientation in turbulent flow inside a planar contraction with different shapes

In this study, we have investigated the influence of shape of planar contractions on the orientation distribution of stiff fibers suspended in turbulent flow. To do this, we have employed a model for the orientational diffusion coefficient based on the data obtained by high-speed imaging of suspension flow at the centerline of a contraction with flat walls. This orientational diffusion coefficient depends only on the contraction ratio and turbulence intensity. Our measurements show that the turbulence intensity decays exponentially independent of the contraction angle. This implies that the turbulence variation in the contraction is independent of the shape, consistent with the results by the rapid distortion theory and the experimental results of axisymmetric contractions. In order to determine the orientation anisotropy, we have solved a Fokker–Planck type equation governing the orientation distribution of fibers in turbulent flow. Although the turbulence variation and the orientational diffusion are independent of the contraction shape, the results show that the variation of the orientation anisotropy is dependent on shape. This can be explained by the variation of the rotational Péclet number, Pe_r, inside the contractions. This quantity is a measure of the importance of the mean rate of the strain relative to the orientational diffusion. We have shown that when Pe_r < 10 turbulence can significantly influence the evolution of the orientation anisotropy. Since in contractions with identical inlet conditions the streamwise position where Pe_r = 10 depends on the shape, the orientation anisotropy is dependent on the variation of rate of strain in a given contraction. We demonstrate the shape effect by considering contraction with flat walls as well as three contractions with different mean rate of strain variation.     

Measurement of heat transfer at the phase interface of condensing bubbles

Using holographic interferometry and high-speed cinematography the heat transfer at the phase interface of vapor bubbles condensing in a subcooled liquid of the same substance was measured with ethanol, propanol, refrigerant R113 and water. To evaluate the axisymmetric temperature field around the bubble from the interference fringe field, the Abel integral method is not sufficient. A correction procedure considering the light deflection caused by the local temperature gradient has been developed and applied to calculate the heat transfer coefficient. The measurements were performed in the range 2 < Pr < 15 and I < Ja < 120. Measured data could be well-correlated as functions of the Prandtl and Reynolds numbers for the heat transfer coefficient and as functions of the Reynolds, Prandtl and Jacob numbers for the bubble collapse time.     

Modal interactions and energy transfers in isotropic turbulence as predicted by local energy transfer theory

Energy transfer processes in decaying, three-dimensional, isotropic turbulence are investigated using numerical results from local energy transfer (LET) theory. The study covers a wide range of evolved, microscale Reynolds numbers (5 < R_λ < 250). It is found that the energy transfer is mainly local (between scales of similar size), but there are also some signs of nonlocal transfer at higher Reynolds number. The nature of the underlying triad-wavenumber interactions, on the other hand, seems to depend on both the Reynolds number and the wavenumber range of interest. In the energy containing and dissipation ranges, both local (all three scales of the triad interaction are of comparable size) and nonlocal (one scale is much larger than the remaining two) interactions are important, with the latter becoming more dominant as the Reynolds number increases. But our nonlocal interactions tend to be less severe than those observed by Domaradzki and Rogallo. More significantly, in the inertial range of high Reynolds number flows, the LET theory predicts dominance of local and near-local interactions. While this is contrary to the result from eddy damped quasi-normal Markovian theory that the important triad interactions are mainly nonlocal, it is closer to the Kolmogorov picture of turbulence. Another interesting result is that, despite their inherent differences, the LET theory and the full simulation of Ohkitani and Kida predict inertial-range values for the energy transfer locality function in fairly good agreement, not only with each other, but also with the analytical closure theory result for infinite Reynolds number, stationary turbulence by Kraichnan. The calculated values reveal that the contributions to the net energy transfer are predominantly from near-local interactions (scale ratios ≈ 4), which is indicative of cancellation of large numbers of highly nonlocal interactions.     

Laminar convective heat transfer of non-Newtonian nanofluids with constant wall temperature

Nanofluids are obtained by dispersing homogeneously nanoparticles into a base fluid. Nanofluids often exhibit higher heat transfer rate in comparison with the base fluid. In the present study, forced convection heat transfer under laminar flow conditions was investigated experimentally for three types of non-Newtonian nanofluids in a circular tube with constant wall temperature. CMC solution was used as the base fluid and γ-Al₂O₃, TiO₂ and CuO nanoparticles were homogeneously dispersed to create nanodispersions of different concentrations. Nanofluids as well as the base fluid show shear thinning (pseudoplastic) rheological behavior. Results show that the presence of nanoparticles increases the convective heat transfer of the nanodispersions in comparison with the base fluid. The convective heat transfer enhancement is more significant when both the Peclet number and the nanoparticle concentration are increased. The increase in convective heat transfer is higher than the increase caused by the augmentation of the effective thermal conductivity.     

Turbulent forced convection heat transfer of non-Newtonian nanofluids

Forced convection heat transfer of non-Newtonian nanofluids in a circular tube with constant wall temperature under turbulent flow conditions was investigated experimentally. Three types of nanofluids were prepared by dispersing homogeneously γ-Al₂O₃, TiO₂ and CuO nanoparticles into the base fluid. An aqueous solution of carboxymethyl cellulose (CMC) was used as the base fluid. Nanofluids as well as the base fluid show shear-thinning (pseudoplastic) rheological behavior. Results indicate that the convective heat transfer coefficient of nanofluids is higher than that of the base fluid. The enhancement of the convective heat transfer coefficient increases with an increase in the Peclet number and the nanoparticle concentration. The increase in the convective heat transfer coefficient of nanofluids is greater than the increase that would be observed considering strictly the increase in the effective thermal conductivity of nanofluids. Experimental data were compared to heat transfer coefficients predicted using available correlations for purely viscous non-Newtonian fluids. Results show poor agreement between experimental and predicted values. New correlation was proposed to predict successfully Nusselt numbers of non-Newtonian nanofluids as a function of Reynolds and Prandtl numbers.     

A prediction model for the effective thermal conductivity of nanofluids considering agglomeration and the radial distribution function of nanoparticles

Conventional heat transfer fluids usually have low thermal conductivity, limiting their efficiency in many applications. Many experiments have shown that adding nanosize solid particles to conventional fluids can greatly enhance their thermal conductivity. To explain this anomalous phenomenon, many theoretical investigations have been conducted in recent years. Some of this research has indicated that the particle agglomeration effect that commonly occurs in nanofluids should play an important role in such enhancement of the thermal conductivity, while some have shown that the enhancement of the effective thermal conductivity might be accounted for by the structure of nanofluids, which can be described using the radial distribution function of particles. However, theoretical predictions from these studies are not in very good agreement with experimental results. This paper proposes a prediction model for the effective thermal conductivity of nanofluids, considering both the agglomeration effect and the radial distribution function of nanoparticles. The resulting theoretical predictions for several sets of nanofluids are highly consistent with experimental data.     

Critical heat flux and turbulent mixing in hexagonal tight rod bundles

Experimental and theoretical investigations have been performed on critical heat flux (CHF) and turbulent mixing in tight, hexagonal, 7-rod bundles. Freon-12 was used as working fluid due to its low latent heat, low critical pressure and well known properties. It has been found that the two-phase mixing coefficient depends mainly on mass flux. It increases with decreasing mass flux and ranges from 0.01 to 0.04 for the test conditions considered. More than 900 CHF data points have been obtained in a large range of parameters: pressure 1.0–3.0 MPa and mass flux 1.0–6.0 Mg/m²s. The effect of different parameters on CHF has been analysed. It has been found that the effect of pressure, mass flux and vapour quality on CHF is similar to that observed in circular tubes. Nevertheless, the CHF in the tight rod bundle is much lower than that in a circular tube of the same equivalent hydraulic diameters. The effect of wire wraps on CHF is mainly dependent on local vapour qualities and subsequently on flow regimes. Based on subchannel flow conditions, the effect of radial power distribution on CHF is small. Comparison of the test results with CHF prediction methods underlines the need for further work.     

Experimental investigations and theoretical determination of thermal conductivity and viscosity of Al2O3/water nanofluid

Experimental investigations and theoretical determination of effective thermal conductivity and viscosity of Al₂O₃/H₂O nanofluid are reported in this paper. The nanofluid was prepared by synthesizing Al₂O₃ nanoparticles using microwave assisted chemical precipitation method, and then dispersing them in distilled water using a sonicator. Al₂O₃/water nanofluid with a nominal diameter of 43 nm at different volume concentrations (0.33–5%) at room temperature were used for the investigation. The thermal conductivity and viscosity of nanofluids are measured and it is found that the viscosity increase is substantially higher than the increase in thermal conductivity. Both the thermal conductivity and viscosity of nanofluids increase with the nanoparticle volume concentration. Theoretical models are developed to predict thermal conductivity and viscosity of nanofluids without resorting to the well established Maxwell and Einstein models, respectively. The proposed models show reasonably good agreement with our experimental results.     

Impingement heat transfer at a circular cylinder due to a submerged slot jet of water

An experimental investigation was carried out on the heat transfer due to a submerged slot jet of water impinging on a circular cylinder in crossflow. The cylinder diameter and the slot width are of the same order of magnitude, specifically D_s = 2.0 and 3.0 mm and D_c = 2.5 and 3.0 mm. The experimental apparatus allowed variation of the slot width, the cylinder diameter, and the distance from nozxle exit to heater. Conditions of impingement from the bottom (ascending flow) were taken into consideration as well as impingement from above (descending flow). The Nusselt number was determined as a function of Reynolds and Prandtl numbers in the range 1.5 × 10³ < Re < 2.0 × 10⁴, 2.7 < Pr < 7.0, and 1.5 ≤ z/D_s ≤ 10. The experimental data were correlated with a simple equation that fits 90% of the data with a precision of 20%.     

Effect of variable thermal conductivity and viscosity on single phase convective heat transfer in slip flow

For a variety of fields in which micro-mechanical systems and electronic components are used, fluid flow and heat transfer at the microscale needs to be understood and modeled with an acceptable reliability. In general, models are prepared by making some extensions to the conventional theories by including the scaling effects that become important for microscale. Some of these effects are; axial conduction, viscous dissipation, and rarefaction. In addition to these effects, temperature variable thermal conductivity and viscosity may become important in microscale gas flows due to the high temperature gradients that may exist in the fluid. For this purpose, simultaneously developing, single phase, laminar and incompressible air flow in a microtube and in the micro gap between parallel plates is numerically analyzed. Navier–Stokes and energy equations are solved and the variation of Nusselt number along the channel is presented in tabular and graphical forms as a function of Knudsen, Peclet, and Brinkman numbers, including temperature variable thermal conductivity and viscosity.     

Effect of an Inserted Porous Layer Located at a Wall of a Parallel Plate Channel on Forced Convection Heat Transfer

A theoretical study is performed on heat and fluid flow in partially porous medium filled parallel plate channel. A uniform symmetrical heat flux is imposed onto the boundaries of the channel partially filled with porous medium. The dimensional forms of the governing equations are solved numerically for different permeability and effective thermal conductivity ratios. Then, the governing equations are made dimensionless and solved analytically. The results of two approaches are compared and an excellent agreement is observed, indicating correctness of the both solutions. An overall Nusselt number is defined based on overall thermal conductivity and difference between the average temperature of walls and mean temperature to compare heat transfer in different channels with different porous layer thickness, Darcy number, and thermal conductivity ratio. Moreover, individual Nusselt numbers for upper and lower walls are also defined and obtained. The obtained results show that the maximum overall Nusselt number is achieved for thermal conductivity ratio of 1. At specific values of Darcy number and thermal conductivity ratio, individual Nusselt numbers approach to infinity since the value of wall temperatures approaches to mean temperature.     

Experimental investigation of heat transfer and pressure drop characteristics of flow through spirally fluted tubes

Spirally fluted tubes are used extensively in the design of tubular heat exchangers. In previous investigations, results for tubes with flute depths e/D_vi < 0.2 were reported, with most correlations applicable for Re ≥ 5000. This paper presents the results of an experimental investigation of the heat transfer and pressure drop characteristics of spirally fluted tubes with the following tube and flow parameter ranges: flute depth e/D_vi = 0.1−0.4, flute pitch p/D_vi = 0.4−7.3, helix angle θ/90° = 0.3−0.65, Re = 500−80,000, and Pr = 2−7. The heat transfer coefficients inside the fluted tube were obtained from measured values of the overall heat transfer coefficient using a nonlinear regression scheme. The friction factor data obtained consisted of 507 data points. The proposed correlation for the friction factor predicts 96% of the database within ±20%. The heat transfer correlation for the range 500 ≤ Re ≤ 5000 predicts 76% of the database (178 data points) within ±20%, and the correlation for the higher Re range predicts 97% of the 342 data points within ±20%. Comparison of heat transfer and friction data show that these tubes are most effective in the laminar and transition flow regimes. The present results show that the increase of flute depth in the range considered does not improve heat transfer.     

Measurement of the effective thermal conductivity of a Mg–MgH2 packed bed with oscillating heating

The magnesium–magnesium hydride–hydrogen-system (Mg–MgH₂–H₂) offers, because of its combined hydrogen and heat storage capacity, the possibility to design hydride heat pumps and heat stores. For such industrial application systems based on cylindrically formed reactors filled with an active magnesium powder, the effective thermal conductivity limits the time in which the metal hydride alloy is charged and discharged with hydrogen. Determination of this transport coefficient is of fundamental importance for the optimum design of magnesium hydride reactors. The complex interrelation of the different transport mechanisms in a metal hydride packed bed and the hitherto undefined rule that the solid effective thermal conductivity behaves as a function of the hydrogen concentration, requires a reliable and simple-to-realize measuring method so as to determine the effective thermal conductivity of a magnesium hydride bed. In the present study, a report is given for the first time on the initiation of a measuring technique with oscillating change of temperature in a non-permeated packed bed of fine-grained material. The measurement of the effective thermal conductivity can ensue by tailoring the problem-specific mathematical result to the experimentally recorded temperature-time function. The effective thermal conductivity of the magnesium hydride bed varies between 2 and 8 W/(m K) in a temperature range of 523–653 K.     

Influence of medium viscosity and adsorbed polymer on the reversible shear thickening transition in concentrated colloidal dispersions

The influence of medium viscosity on the onset of shear thickening of silica dispersions is investigated with two different methods. In the first method, the sample temperature is varied over a narrow range for two different suspensions. For the first suspension, the stress at the onset of shear thickening, or the critical stress, was found to be independent of sample viscosity, and the shear viscosity scaled with Peclet number, as expected. The critical stress for the second suspension was not independent of sample viscosity, and the Peclet number scaling was only moderately successful. The differences were attributed to changes in particle interactions with temperature. In the second method, the molecular weight of an oligomeric silicone oil medium is varied. In principle, this method should maintain constant chemical interactions as medium viscosity varies; however the polymer is found to adsorb onto the silica surface and delay shear thickening to higher stresses with increasing molecular weight. The critical stress for the highest molecular weight systems, which is highly dependent on particle loading, overlays with an effective volume fraction based on the hydrodynamic diameter of the polymer-stabilized colloids. The results of both methods suggest that if all other properties of the dispersion are held constant, critical stress is independent of medium viscosity.     

Numerical analysis of transient laminar forced convection of nanofluids in circular ducts

In this study, forced convection heat transfer characteristics of nanofluids are investigated by numerical analysis of incompressible transient laminar flow in a circular duct under step change in wall temperature and wall heat flux. The thermal responses of the system are obtained by solving energy equation under both transient and steady-state conditions for hydro-dynamically fully-developed flow. In the analyses, temperature dependent thermo-physical properties are also considered. In the numerical analysis, Al₂O₃/water nanofluid is assumed as a homogenous single-phase fluid. For the effective thermal conductivity of nanofluids, Hamilton–Crosser model is used together with a model for Brownian motion in the analysis which takes the effects of temperature and the particle diameter into account. Temperature distributions across the tube for a step jump of wall temperature and also wall heat flux are obtained for various times during the transient calculations at a given location for a constant value of Peclet number and a particle diameter. Variations of thermal conductivity in turn, heat transfer enhancement is obtained at various times as a function of nanoparticle volume fractions, at a given nanoparticle diameter and Peclet number. The results are given under transient and steady-state conditions; steady-state conditions are obtained at larger times and enhancements are found by comparison to the base fluid heat transfer coefficient under the same conditions.