DETERMINATION OF THE EFFECTIVE RADIAL THERMAL DIFFUSIVITY FOR EVALUATING ENHANCED HEAT TRANSFER IN TUBES UNDER NON-NEWTONIAN LAMINAR FLOW

Abstract

Enhanced heat transfer in tubes under laminar flow conditions can be found in coils or corrugated tubes or in the presence of high wall relative roughness, curves, pipe fittings or mechanical vibration. Modeling these cases can be complex because of the induced secondary flow. A modification of the Graetz problem for non-Newtonian power-law flow is proposed to take into account the augmented heat transfer by the introduction of an effective radial thermal diffusivity. The induced mixing was modeled as an increased radial heat transfer in a straight tube. Three experiments using a coiled tube and a tubular heat exchanger with high relative wall roughness are presented in order to show how this parameter can be obtained. Results were successfully correlated with Reynolds number. This approach can be useful for modeling laminar flow reactors (LFR) and tubular heat exchangers available in the chemical and food industries.

Keywords
Heat transfer; Laminar flow; Diffusivity; Laminar flow reactor; Graetz problem

INTRODUCTION

Heat transfer to non-Newtonian viscous fluids under laminar flow in tubes is a recurrent problem in many chemical and food engineering applications, such as laminar flow reactors (LFR), tube sterilizers or tubular heat exchangers. Heat is transferred through the liquid by conduction in the radial direction and by advection and conduction in the axial direction. Due to the high consistency of the fluid, free convection is seldom significant (Skelland, 1967; Chhabra and Richardson, 2008).

The well known Graetz problem of laminar flow heat exchange in a straight tube has been extensively studied and various correlations for the Nusselt number, either theoretical or empirical, are available in the literature for Newtonian and non-Newtonian fluids (Lyche and Bird, 1956; Jones, 1971; Mahalingam et al., 1975; Richardson, 1979; Weigand and Gassner, 2007). In the case of coiled tubes or corrugated tubes under laminar flow, the development of secondary flows improves fluid mixing and, consequently, heat transfer. However, the mathematical modeling of these cases is more complex and studies often rely on empirical correlations (Barba et al., 2002; Vicente et al., 2004; Naphon, 2007; Pimenta and Campos, 2013; Pawar and Sunnapwar, 2013).

Enhanced heat transfer in laminar flow can also be associated with high relative wall roughness of the tube and the presence of curves, pipe fittings or mechanical vibration (Bergles and Joshi , 1983). In order to model these systems, the purpose of this work was to adapt the Graetz problem to take into account the heat transfer augmentation in the form of an effective heat transfer diffusivity. In this case, a simpler model could be used to simulate heat transfer in real tubular systems. Three experiments are presented to show how this effective diffusivity can be determined through this simple approach and what the extent of the heat transfer enhancement can be.

An analogous approach was taken by Dantas et al. (2014) to model mass dispersion in laminar flow tubular reactors. Dispersion inducing elements were assumed to increase the radial mass diffusivity in the reactor. By comparing experimental residence time distribution (RTD) data with the model prediction of tracer dispersion in the reactor, it was possible to determine the effective radial mass diffusivity and the mean residence time. With these parameters, the LFR model would provide a more reliable prediction of the RTD of the system, without the need of modeling complex flow distributions in the reactor.

MATHEMATICAL MODELING

The objective is to derive a dimensionless model for the problem of heat transfer to a fluid flowing in a tube under developed laminar flow. A tube with circular cross section, length L and internal radius R is considered. The temperature of the internal surface of the tube is assumed to be uniform (T_w). Firstly, the Graetz problem is derived for non-Newtonian power-law flow in a straight smooth tube with radial diffusion of heat and then the effective radial thermal diffusivity is introduced as a parameter to take into account mechanisms that promote mixing and enhance heat transfer in laminar flow, such as coiled or corrugated tubes.

The following equation represents the thermal energy conservation for a homogenous liquid with a velocity field Bird et al., 2002; Incropera et al., 2006): and a temperature gradient (

where ρ is the density, C_P is the specific heat capacity at constant pressure, k is the thermal conductivity, T is the temperature, t is the time and Eq. (1), from left to right, represent: 1) the accumulation of thermal energy, 2) the transport of thermal energy due to advection, 3) the transport of thermal energy due to diffusion and 4) the generation or consumption of thermal energy. is the uniform volumetric thermal energy generation. The terms in

Considering steady-state conditions, uniform thermo-physical properties and negligible thermal energy generation (viscous dissipation), Eq. (1) can be expressed as:

where α = k/(ρ C_P) is the thermal diffusivity of the fluid.

Assuming incompressible inelastic flow with developed velocity profile, the fluid velocity would have only one axial component, v_z(r). Moreover, the temperature field can be considered to be axisymmetric around the z axis, i.e., T = T(z,r). Regarding the axial conduction of heat, this transport mechanism can be neglected for most engineering applications. In this case, Eq. (2) is reduced to (Chhabra and Richardson, 2008; Bird et al., 2002):

which is known as the Graetz problem for heat transfer in laminar flow.

In order to solve this differential equation, three boundary conditions are required; which are: 1) specification of the tube inlet temperature (T_in), 2) the symmetry condition in the center of the tube, and 3) specification of the tube wall temperature (T_w). These conditions are summarized as:

The laminar velocity profile depends on the fluid viscosity. The rheological parameters of a power-law fluid are the consistency coefficient K and the flow behavior index n. Since average thermo-physical properties are considered, the temperature effects on the rheological parameters are not taken into account. Consequently, the following power-law velocity profile can be assumed as uniform along the tube (Chhabra and Richardson, 2008):

where v_m is the bulk velocity, which is related to the tube pressure drop and the consistency coefficient.

The Graetz problem for power-law laminar flow can be expressed in this dimensionless form:

using the following dimensionless variables:

and the parameter Pe:

which is the reactor Peclet number. This parameter represents the ratio between the axial transport of thermal energy by advection and the radial transport of thermal energy by diffusion (Fogler, 2006).

By introducing the modified Peclet number for radial diffusion of heat defined as:

the problem can be expressed in the following form:

This dimensionless problem has two parameters: the flow behavior index of the fluid (n) and the modified Peclet number (Pe'). Alternatively, the modified Peclet number can be expressed as the Graetz number (Dantas et al., 2014):

where Re is the Reynolds number, Pr is the Prandtl number, D = 2 R is the tube diameter and μis the fluid viscosity. However, the Peclet form seems more suitable for diffusivity studies.

The average flow behavior index of the fluid (n) can be obtained from a rheological study correlating shear rate and shear stress for proper temperature and shear conditions. In the case of a Newtonian fluid, n = 1.0. The modified Peclet number for radial diffusion of heat (Pe') can be calculated through Eq. (9) using an average value for the heat diffusivity of the fluid.

The solution of the model provides the temperature distribution in the tube, θ(Z,Y). The mixing cup average temperature in the tube (T_m or θ_m) can be evaluated numerically with:

The average Nusselt number based on the tube diameter can be calculated with:

where h is the average convective coefficient. In order to obtain this dimensionless expression for Nu, the following definition of the convective coefficient was employed:

where q is the heat transfer rate, A is the heat transfer area, ΔT is the inlet/outlet change in T_m, ΔT_lm is the log-mean of the temperature difference between T_m and T_w (LMTD) and Q is the volumetric flow rate (Bird et al., 2002; Incropera et al., 2006).

Heat Transfer Enhancement

As discussed in the introduction, the heat transfer to the fluid can be enhanced when using coiled or corrugated tubes. Mixing in laminar flow can also be induced by high relative wall roughness, mechanical vibration or the presence of fittings or curves. Since the mathematical modeling of these systems can be rather complex (Bergles and Joshi, 1983; Pimenta and Campos, 2013), the use of the straight tube model derived in this section with an "effective radial thermal diffusivity" (α_eff ≥ α) is suggested. The augmentation of the heat transfer due to mixing in laminar flow could be modeled as a higher thermal diffusivity acting in the radial direction. The physical phenomena are different, but this approach could lead to a good approximation of the heat transfer enhancement using a simpler mathematical model. This would leave the model with one empirical parameter to be determined through heat transfer experiments in the equipment. In this case, the theoretical Pe' should be replaced in the model by:

which is an experimental parameter and depends on the effective radial thermal diffusivity.

When the outlet temperature of the fluid is measured (θ_m(Z=1)), the model in Eqs. (10) and (11) can be solved to determine the corresponding modified Peclet number. Equation (16) would then provide the effective radial thermal diffusivity (α_eff). Since α_eff = k_eff /(ρ C_P), the improvement in heat transfer could also be expressed as an effective radial thermal conductivity (k_eff). Moreover, a heat transfer enhancement factor can be defined as:

This factor would be useful for the evaluation of corrugated or coiled tubes for improved heat transfer in laminar flow.

MATERIALS AND METHODS

In order to test and validate the proposed model, three heat transfer experiments (A, B1 and B2) were conducted and the enhancement factor in Eq. (17) was determined and analyzed. In Experiment A, a stainless-steel coiled tube with 9.3 mm internal diameter, external 12.7 mm diameter, 200 mL volume and 2.94 m calculated linear length was used. The diameter of the coil was 10.5 cm (eight turns). Thermocouples with exposed junctions (IOPE, Brazil) were placed in the inlet and outlet of the tube using union tees (John Guest, UK) and were connected to a cDAQ-9172 data acquisition system with three NI9211 modules (National Instruments, USA) and a computer running a data-logger software. The outlet thermocouple was placed in a 90° bend to provide mixing to the fluid leaving the coil. Thermocouples were previously calibrated for the temperature range of 0 to 90 °C.

The coil was connected to the hot water circulating bath of a FT-43A laboratory plate pasteurizer (Armfield, UK). The equipment has a centrifugal pump, a rotameter and a heating element (1.5 kW). Tested volumetric flow rates were Q = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 L/min of distilled water. The coil was immersed in an ice-water bath (27 L) with two propeller mixers and ice was added constantly in order to maintain the bath temperature constant and uniform. The temperature of the external surface of the tube (T_ext) was assessed with an exposed junction thermocouple. The average temperature of the internal surface of the tube (T_w) was calculated taking into account the thermal resistance of a cylindrical wall (Incropera et al., 2006):

where R_ext is the external radius of the tube and k_tube is the thermal conductivity of the metal.

The heat transfer rate across the tube wall (q) was calculated with:

where T_in and T_out are the measured inlet and outlet temperatures of the fluid. The experimental Nusselt number was calculated based on the definition in Eq. (14) with the fluid thermal diffusivity.

Once steady-state operation was verified, the temperatures were recorded for 30 s and average values were calculated. Two tests were performed for each flow rate, with inlet temperatures between 50 °C and 65 °C.

For Experiments B1 and B2, a stainless steel double-pipe heat exchanger was used. This exchanger is the cooling section of a liquid food thermal processing unit (laboratory scale). More details are presented by Dantas et al. (2014). The inner tube has an internal diameter of 4.0 mm and an external diameter of 6.0 mm. The annular space has an internal diameter of 21.0 mm. Flow was provided by an eccentric screw pump model 3NE10A (Netzsch, Brazil) and tested flow rates were Q = 20, 30, 40 and 50 L/h. The fluids tested were a 80/20 (w/w) Glycerin/Water mixture in Experiment B1, which has Newtonian behavior, and a 1.0% (w/w) solution of carboxymethylcellulose (CMC) in Experiment B2, which is a pseudoplastic fluid.

The exchanger of the heating section worked with pressurized hot water and steam injection in a closed loop circuit and a 32-16A centrifugal pump (KSB, Brazil). The cooling section, which was monitored in this study, was connected to a MAS-5-RI-220/C chiller (Mecalor, Brazil) with a cooling power of 5.8 kW and a flow rate between 5 and 6 m³/h. The flow rate was kept at its maximum to provide a uniform wall temperature around the internal tube (v_m ≈ 9 m/s in the annulus). The temperature of the external surface of the tube was assumed to be equal to the average temperature of the cooling water because of the high velocity in the annular space. The average temperature of the internal surface was calculated using Eqs. (18) and (19).

The heat exchanger of the cooling section consists of four straight tubes with a total effective length for heat transfer of 6.7 m. The tubes were connected by 180° bends with a diameter of 12.5 cm. The fluid inlet and outlet temperatures were obtained using thermocouples (IOPE, Brazil) and the same data acquisition equipment of Experiment A. Inlet and outlet temperatures of the cooling water were obtained using two thermo-resistances and a temperature display (H.Cidade, Brazil). Once steady-state operation was observed, the temperatures were measured for 30 s and average values were calculated.

The dimensionless model presented in Eqs. (10) and (11) was solved using the finite difference method to discretize the variables. For n = 1.0, a first-order backward scheme with 150 equally spaced intervals was used for variable Z (axial) while a second order centered scheme with 50 intervals was used for variable Y (radial). For the simulations with n < 1.0 or log(Pe') ≤ 0.0, the number of intervals had to be increased. Different numbers of intervals were tested in order to achieve reliable results with less computational time. The model was implemented in the software gPROMS 3.2 (Process Systems Enterprise, UK). For each simulation, values of the parameters n and Pe' were specified.

For determining the experimental values of the modified Peclet number (Pe'_exp), this parameter was treated as a variable in the simulation and the measured outlet temperature was specified as θ_m(Z=1) to yield zero degrees of freedom. The enhancement factor was subsequently calculated using Eqs. (16) and (17). Results were correlated with the Reynolds number to analyze the behavior of the heat transfer enhancement. The Reynolds number for tube flow was evaluated with:

For the power-law non-Newtonian flow, Re was calculated using the concept of generalized viscosity μ_g instead of the fluid viscosity μ in Eq. (20). The generalized viscosity is defined as:

where K is the consistency coefficient of the power-law model, D_eq is the equivalent diameter of the duct and ξ and υ are geometrical parameters of the duct (D_eq = 2 R, ξ = 8 and υ = 3 for a circular tube) (Delplace and Leuliet, 1995; Carezzato et al., 2007).

Thermo-physical properties of water were evaluated using the correlations compiled by Gut and Pinto (2003). Density, heat capacity and thermal conductivity of pure glycerin were obtained from Yaws (2003) and the properties of the glycerin/water mixture were weighted based on the volume fraction (ρ and k) or on the mass fraction (C_P). Viscosity of the glycerin/water mixture was obtained from Cheng (2008). Properties of the CMC solution were determined by this research group and are available elsewhere (Carezzato et al., 2007). All properties were evaluated at inlet and outlet conditions and the average was used in the calculations.

RESULTS AND DISCUSSION

Model Simulation Results

The proposed model was simulated for different values of n and Pe'. Figure 1 presents simulation results showing the average dimensionless temperature along the tube length as affected by the modified Peclet number for a Newtonian fluid (n = 1.0) and a pseudoplastic fluid with n = 0.3 for comparison.

Figure 1
Simulation results for the average temperature (θ_m) as a function of the tube length (Z) and modified Peclet number (Pe') for n = 1.0 (black lines) and n = 0.3 (grey lines).

Lower values of Peclet are related to high thermal diffusivities and, thus, to a more efficient heat transfer. Figure 2(a) shows the effect of Pe' on the outlet temperature of the tube. Pseudoplasticity is verified to slightly improve the heat exchange, which can be explained by the flatter velocity profile.

Figure 2
Simulation results showing the effect of the modified Peclet number (Pe') on (a) the tube outlet temperature and (b) on the Nusselt number for n = 1.0 (black lines) and n = 0.3 (grey lines).

Figure 2(b) presents the effect of Pe' on the Nusselt number and, for decreasing Pe', the curve is observed to converge asymptotically to the constant Nusselt value of fully developed flow (velocity and thermal profiles). For n = 1.0, limiting Nusselt in Nu = 3.66 (Incropera et al., 2006). Low values of Pe' can be associated with long tubes, in which the entrance effect can be neglected.

Model Validation

Table 1 presents the range of the variables in the experiments. For Experiment A, the Reynolds number varied between 885 at 0.3 L/min to 3856 at 1.0 L/min, thus including the laminar and transition flow regimes. For all conditions, the experimental value of the modified Peclet number was obtained through model simulation with the specification of the outlet dimensionless temperature of the tube as described in the previous section. The enhancement factor was calculated through Eq. (17).

Table 1
Range of the variables in the experimental tests.

Experiment	A: Coiled tube with water	B1: Double pipe exchanger with glycerin/water mixture	B2: Double pipe exchanger with 1.0% CMC solution
Tin (°C)	52 to 64	60 to 84	58 to 85
Tout (°C)	17 to 32	14 to 25	9.3 to 18
Tw (°C)	16 to 37	6.2 to 8.1	5.7 to 10
q (kW)	0.81 to 1.4	1.0 to 2.4	1.3 to 3.1
vm (m/s)	0.074 to 0.25	0.44 to 1.1	0.44 to 1.1
n (-)	1.0	1.0	0.43 to 0.47
Re (-)	885 to 3856	42 to 200	37 to 170
Pe' exp (-)	0.77 to 2.4	1.7 to 3.3	1.2 to 2.2
Fheat (-)	3.5 to 6.6	1.3 to 2.0	1.1 to 2.1

The results of Experiment A are presented in Figure 3 as a function of the Reynolds number. As expected, there was a significant improvement in the heat transfer because of the coil curves that induce secondary flows and mixture. This problem could be modeled as a straight tube under laminar flow with an effective thermal diffusivity up to six times higher than the actual thermal diffusivity of water. The increase in the flow rate caused an increase in the enhancement factor because of the larger degree of mixing. Because of the small data range and experimental error it was not possible to observe in Figure 3 the end of the laminar flow regime and the beginning of the transition flow regime.

Four tests were conducted for each flow rate using the double pipe heat exchanger for Experiments B1 and B2. The temperature increase of the cold water stream was less than 0.9 °C because of the high flow rate. The low Reynolds numbers in Table 1 suggest laminar flow regime.

Figure 3
Heat transfer enhancement factor (F_heat) as a function of the Reynolds number (Re) for Experiment A (water cooling in the coiled tube).

Figure 4 presents the enhancement factor obtained in these experiments as a function of the Reynolds number. Once again, an improved heat transfer was observed (F_heat > 1.0) and the enhancement factor also increased with higher flow rates. The effective radial thermal diffusivity was up to two times higher than the fluid thermal diffusivity. This higher heat transfer efficiency in Experiments B1 and B2 can be attributed to the large relative roughness of the tube: 0.0018 in / 4.0 mm = 0.011 (Perry and Chilton, 1973). The wall turbulence induced by the tube roughness makes an important contribution to the heat transfer. Moreover, the presence of the 180° bends and the tip of thermocouples along the fluid path also contributed to the fluid mixing.

Figure 4
Heat transfer enhancement factor (F_heat) as a function of the Reynolds number (Re) for Experiments (a) B1 (glycerin/water mixture) and (b) B2 (CMC solution) in the double pipe heat exchanger.

Dantas et al. (2014) determined the effective mass diffusivity in the same equipment and with the same fluids used in Experiments B1 and B2. RTD data for the isothermal flow of a 80/20 (w/w) glycerin/water mixture and 1.0% (w/w) solution of carboxymethylcellulose (CMC) was obtained for flow rates between 10 and 50 L/h in one hairpin of the exchanger with 3.64 m length. Methylene blue dissolved in the test fluid was used as tracer and its concentration was assessed by a spectrophotometer. Process modeling was used to simulate the flow dispersion of the tracer in a LFR with radial mass diffusivity and power-law velocity profile. Comparison between experimental data and simulation results provided the effective mass diffusivity, also expressed as a modified Peclet number. Results showed that the radial dispersion was higher than expected in this equipment, as was verified in this work on thermal diffusion. This behavior was also attributed to the large relative roughness of the tube.

CONCLUSIONS

The enhanced heat transfer in tubes under laminar flow with induced mixing by the presence of coils, corrugated walls, curves, pipe fittings, high relative wall toughness or mechanical vibration was modeled as a straight tube problem with an effective radial thermal diffusivity. This variable was expressed as a modified Peclet number, a parameter to be determined from heat transfer experiments.

In the three experiments reported herein, it was possible to determine the modified Peclet number and correlate it with the Reynolds number of the flow. For the flow in the coiled tube, the effective thermal diffusivity was up to six times higher than the fluid property. In the case of the double-pipe heat exchanger with high relative wall toughness, the thermal diffusivity almost doubled.

For these experiments, the proposed model could predict the outlet temperature of the tube for a given effective radial thermal diffusivity or enhancement factor (adjusted correlations presented in Figures 3 and 4). However, it was not verified if this simplified model could predict the temperature distribution along the tube with confidence. If experimental data of the average temperature are available at more axial positions, the sum of squared errors on the predicted temperature could be minimized in order to adjust Pe'_exp.

The advantage of using the effective thermal diffusivity is that a laminar flow reactor that has mixing inducing elements can be modeled as a straight tube reactor without the need of describing a complex geometry or conditions. The work by Kechichian et al. (2012) can be used as an example, where the thermal processing equipment for liquid foods used in Experiments B1 and B2 was modeled and simulated as a non-isothermal laminar flow reactor to determine the temperature and lethality distribution along the food product path. The model relies on the thermal and mass diffusivities in the food product and the experiments in this work can contribute with important parameters for the process model simulation and validation. Preliminary data have already shown that the thermal effectiveness of the two heat exchangers is higher than expected because of the enhanced heat transfer to the food product (Dantas et al., 2012).

NOMENCLATURE

Greek Letters

ACKNOWLEDGMENTS

The authors would like to acknowledge the financial support from the São Paulo Research Foundation (FAPESP) under grant 2009/07934-3.

REFERENCES

• Barba, A., Rainieri, S. and Spiga, M., Heat transfer enhancement in a corrugated tube. International Communications in Heat and Mass Transfer, 29 (3), p. 313-322 (2002).
• Bergles, A. E. and Joshi, S. D., Augmentation Techniques for Low Reynolds Number in-Tube Flow. In: Kakaç, S., Shah, R. K. and Bergles, A. E. (Eds.), Low Reynolds Number Flow Heat Exchangers. Hemisphere Pub. Co., Washington, p. 695-720 (1983).
• Bird, R. B., Stewart, W. E. and Lightfoot, E. N., Transport Phenomena. 2nd Ed., John Wiley & Sons, New York, p. 290-332 (2002).
• Carezzato A., Alcantara, M. R., Telis-Romero, J., Tadini, C. C. and Gut, J. A. W., Non-Newtonian heat transfer on a plate heat exchanger with generalized configurations. Chemical Engineering and Technology, 30(1) p. 21-26 (2007).
• Cheng, N. S., Formula for the viscosity of a glycerol-water mixture. Industrial & Engineering Chemistry Research, 47(9) p. 3285-3288 (2008).
• Chhabra, R. P. and Richardson, J. F., Non-Newtonian Flow and Applied Rheology: Engineering Applications. 2nd Ed., Butterworth-Heinemann, p. 316-342 (2008).
• Dantas, J. A. T. A., Ferrão, E. S. P. and Gut, J. A. W., Model validation of continuous thermal processing of a pseudoplastic food under non-ideal laminar conditions. EFFoST Annual Meeting Abstracts. Elsevier, Berlin, p. 4.41 (2011).
• Dantas, J. A. T. A., Gut, J. A. W. and Pegoraro, P. R., Determination of the effective radial mass diffusivity in tubular reactors under non-Newtonian laminar flow using residence time distribution data. International Journal of Heat and Mass Transfer, 71, p. 18-25 (2014).
• Delplace, F. and Leuliet, J. C., Generalized Reynolds number for the flow of power law fluids in cylindrical ducts of arbitrary cross-section. The Chemical Engineering Journal, 56(2) p. 33-37 (1995).
• Fogler, H. S., Elements of Chemical Reaction Engineering. 4th Ed., Prentice Hall, Upper Saddle River,p. 945-1005 (2006).
• Gut, J. A. W. and Pinto, J. M., Modeling of plate heat exchangers with generalized configurations. International Journal of Heat and Mass Transfer, 46(14), p. 2571-2585 (2003).
• Incropera, F. P., DeWitt, D. P., Bergman, T. L. and Lavine, A. S., Fundamentals of Heat and Mass Transfer. 6th Ed., John Wiley & Sons, New York, p. 485-558 (2006).
• Jones, A. S., Extensions to the solution of the Graetz problem. International Journal of Heat and Mass Transfer, 14(4), p. 619.623 (1971).
• Kechichian, V., Crivellari, G. P., Gut J. A. W. and Tadini C. C., Modeling of continuous thermal processing of a non-Newtonian liquid food under diffusive laminar flow in a tubular system. International Journal of Heat and Mass Transfer, 55(21-22), p. 5783-5792 (2012).
• Lyche, B. C. and Bird, R. B., Graetz-Nusselt problem for a power-law non-Newtonian fluid. Chemical Engineering Science, 6, p. 35-41 (1956).
• Mahalingam, R., Tilton, L. O. and Coulson, J. M., Heat transfer in laminar flow of non-Newtonian fluids. Chemical Engineering Science, 30(8), p. 921-929 (1975).
• Naphon, P., Laminar convective heat transfer and pressure drop in the corrugated channels. International Communications in Heat and Mass Transfer, 34(1), p. 62-71 (2007).
• Pawar, S. S. and Sunnapwar, V. K., Experimental studies on heat transfer to Newtonian and non-Newtonian fluids in helical coils with laminar and turbulent flow. Experimental Thermal and Fluid Science, 44, p. 792-804 (2013).
• Perry, R. H. and Chilton, C. H., Chemical Engineers' Handbook. 5th Ed., McGraw-Hill Kogakusha, Tokyo, p. 5.20-5.40 (1973).
• Pimenta, T. A. and Campos, J. B. L. M., Heat transfer coefficients from Newtonian and non-Newtonian fluids flowing in laminar regime in a helical coil. International Journal of Heat and Mass Transfer, 58(1-2), p. 676-690 (2013).
• Richardson, S. M., Extended Leveque solutions for flows of power law fluids in pipes and channels. International Journal of Heat and Mass Transfer, 22(10), p. 1417-1423 (1979).
• Skelland, A. H. P., Non-Newtonian Flow and Heat Transfer. Wiley & Sons, New York, p. 358-433 (1967).
• Vicente, P. G., García, A. and Viedma, A., Mixed convection heat transfer and isothermal pressure drop in corrugated tubes for laminar and transition flow. International Communications in Heat and Mass Transfer, 31(5), p. 651-662 (2004).
• Weigand, B. and Gassner, G., The effect of wall conduction for the extended Graetz problem for laminar and turbulent channel flows. International Journal of Heat and Mass Transfer, 50(5-6), p. 1097-1105 (2007).
• Yaws, C. L., Yaws' Handbook of Thermodynamic and Physical Properties of Chemical Compounds. Knovel (2003).

Publication Dates

• Publication in this collection
  Apr-Jun 2015

History

• Received
  24 Feb 2014
• Reviewed
  15 Aug 2014
• Accepted
  20 Aug 2014