Fluid flow and heat transfer in a two-dimensional wavy channel

A numerical study is presented for the laminar fully developed flow and heat transfer in a two-dimensional wavy channel. The effects of the geometry, Reynolds and Prandtl number on the flow field and heat transfer are investigated. The channel is characterized by a wavy wall, heated at uniform heat flux, and an opposite wall, being plane and adiabatic. The extent of the wall waviness and the distance between the channel walls are found to significantly affect the streamlines contours as well as the heat transfer coefficients. Comparisons with the straight channel, in the same flow rate and heat transfer conditions, have been performed. Pressure drop of the wavy channel is found to be always larger than the value characteristic of a straight channel, while heat transfer performance decreases or increases depending on the values of the parameters (geometry, Reynolds and Prandtl numbers).     

Turbulent flow and convective heat transfer in a wavy wall channel

This paper reports the numerical modeling of turbulent flow and convective heat transfer over a wavy wall using a two equations eddy viscosity turbulence model. The wall boundary conditions were applied by using a new zonal modeling strategy based on DNS data and combining the standard k– turbulence model in the outer core flow with a one equation model to resolve the near-wall region.It was found that the two-layer model is successful in capturing most of the important physical features of a turbulent flow over a wavy wall with reasonable amount of memory storage and computer time. The predicted results show the shortcomings of the standard law of the wall for predicting such type of flows and consequently suggest that direct integrations to the wall must be used instead. Moreover, Comparison of the predicted results of a wavy wall with that of a straight channel, indicates that the averaged Nusselt number increases until a critical value is reached where the amplitude wave is increased. However, this heat transfer enhancement is accompanied by an increase in the pressure drop.     

Mass transfer characteristics in an axisymmetric wavy-walled tube for pulsatile flow with backward flow

Under the pulsatile flow with backward flow (PFBF) conditions, flow mixing and mass transfer characteristics are experimentally investigated in an axisymmetric wavy-walled tube at a net flow Reynolds number from 50 to 1,000. An electrochemical technique is employed to measure the mass transfer rate. An optimal Strouhal number corresponding to the peak value of the mass transfer enhancement factor is observed, which is independent of the oscillatory fraction of the flow rate, but decreases with the increasing net flow Reynolds number. It was found that the mass transfer enhancement under PFBF has the similar characters of resonant enhancement in two-dimensional (2-D) channels, but there also exists an essential difference since no self-sustained oscillation occurs in the wavy-walled tube.     

Numerical investigation of heat transfer and fluid flow characteristics inside a wavy channel

The periodically fully developed laminar heat transfer and fluid flow characteristics inside a two-dimensional wavy channel in a compact heat exchanger have been numerically investigated. Calculations were performed for Prandtl number 0.7, and Reynolds number ranging from 100 to 1,100 on non-orthogonal non-staggered grid systems, based on SIMPLER algorithm in the curvilinear body-fitted coordinates. Effects of wavy heights, lengths, wavy pitches and channel widths on fluid flow and heat transfer were studied. The results show that overall Nusselt numbers and friction factors increase with the increase of Reynolds numbers. According to the local Nusselt number distribution along channel wall, the heat transfer may be greatly enhanced due to the wavy characteristics. In the geometries parameters considered, friction factors and overall Nusselt number always increase with the increase of wavy heights or channel widths, and with the decrease of wavy lengths or wavy pitches. Especially the overall Nusselt number significantly increase with the increase of wavy heights or channel widths, where the flow may become into transition regime with a penalty of strongly increasing in pressure drop. An erratum to this article can be found at     

Natural convection heat transfer from sinusoidal wavy surfaces

This paper presents the results of an experimental study of the natural convection heat transfer characteristics of sinusoidal wavy surfaces on vertical plates maintained at a constant temperature. Local heat transfer coefficients were obtained with a Mach-Zehnder interferometer. The heat transfer from the wavy surfaces, compared to a plane plate of equal projected area, increased with increasing amplitude-to-wavelength ratio. The heat transfer was increased by about 15 percent at an amplitude-to-wavelength ratio of 0.3; for this case a flow instability was detected. A quantitative comparison with a previously published numerical investigation is also presented. In general, there is agreement between the two studies.     

Oscillatory flow and mass transfer within asymmetric and symmetric channels with sinusoidal wavy walls

Mass transfer for oscillatory flow was studied experimentally in channels with two different geometries, i.e., a periodically converging-diverging channel and a serpentine channel, both having sinusoidal wavy walls. The experiments were carried out under the following conditions: 10<Re<500 and 0.008<St<0.05. The channel geometries were found to have an important effect on the flow patterns and the mass transfer rates. At low Strouhal numbers of less than 0.023, the mass transfer rates for both channels were almost identical, regardless of different flow patterns and wall shear stresses. At high Strouhal numbers, however, the serpentine channel had a smaller mass transfer rate than the converging-diverging channel. The mass transfer characteristics were explained in terms of the vortex dynamics, wall shear stresses and fluid mixing based on numerical analysis and flow visualizations. The serpentine channel yields a better mass transfer and pumping power performance than the converging and diverging channel at low Strouhal numbers.

---

Schwingungsströmung und Stofftransport innerhalb asymmetrischer und symmetrischer Kanäle mit sinusförmig gewellten Wandungen

Zusammenfassung Der Stofftransport bei Schwingungsströmungen in Kanälen mit zwei verschiedenen Geometrien experimentell untersucht, d.h. in einem periodisch konvergierenden und divergierenden Kanal und einem schlangenförmig gewundenen Serpentinenkanal, wobei die Kanäle jeweils sinus förmig gewellte Wandungen aufwiesen. Die Versuche wurden unter folgenden Bedingungen, ausgeführt; 10<Re<500 und 0.008<St<0.05. Es zeigte sich, daß die Kanalgeometrien einen erheblichen Einfluß auf die Strömungsmuster und Stofftransportraten haben. Bei niedrigen Strouhal-Zahlen unter 0.023 waren die Stofftransportraten beider Kanäle praktisch identisch, und zwar unabhängig von unterschiedlichen Strömungsmustern und Wand-Scherspannungen. Bei hohen Strouhal-Zahlen dagegen zeigte der Serpentinenkanal eine geringere Stofftransportrate als der Konvergenz-Divergenz-Kanal. Diese Stofftransport-Charakteristik wurde, basierend auf numerischer Analyse und Sichtbarmachung der Strömung, erklärt in Form von Vortexdynamic, Wand-Scherspannungen und Flüssigkeitsmischung. Bei niedrigen Strouhal-Zahlen erbringt der Serpentinenkanal ein besseres Stofftransport- und Pumpleistungsverhalten als der Konvergenz-Divergenz-Kanal.

---

Nomenclature A area of mass transfer surface - a wave amplitude of wavy wall - c _b concentration of the ferricyanide ion - D diameter of piston - D molecular diffusivity of the ferricyanide ion - F Faraday constant - f frequency of oscillation - H _min minimum spacing between wavy walls - Re Reynolds number, Eq. (4) - s length of stroke of piston - Sc Schmidt number - Sh Sherwood number, Eq. (5) - St Strouhal number, Eq. (3) - T period of oscillation - U _p peak velocity based onH _min - W width of wavy wall Greek symbols wavelength of wavy wall - kinematic viscosity - _w wall vorticity - _max vortex strength - angle of misalignment of the two channel walls This work was supported in part by a Grant-in-Aid for Science Research (No. 63750889 and No. 03302031) from the ministry of Education, Science and Culture of Japan. The author acknowledge with thanks the assistance of graduate student Shigeki Matsune in the computations.     

Vortices and wall shear stresses in asymmetric and symmetric channels with sinusoidal wavy walls for pulsatile flow at low Reynolds numbers

Numerical simulation and flow visualization were performed to study the dynamical behavior of vortices generated in channels with two different geometries, i.e., a periodically converging–diverging channel and serpentine channel, both having sinusoidal wavy walls. This system for pulsatile flow is used to enhance heat and mass transfer in very viscous liquids. The numerical results predict well the dynamical behavior of vortices and agree with the flow visualizations. For both channels, the vortex expands in each furrow of the channel walls during the deceleration phase and shrinks during the acceleration phase, which leads to fluid exchange between the vortex and the mainstream. The time-averaged vortex strength and wall shear stresses increase, as the frequency of fluid oscillation increases under a fixed oscillatory fraction of the flow rate. However, above a certain value of the frequency, they reversely decrease due to viscous effects. This frequency for the serpentine channel is smaller than that for the converging–diverging channel. The channel geometries are found to have an important effect of the flow characteristics.     

Slow two-phase flow through a sinusoidal channel

The two-phase flow through a symmetric sinusoidal channel is studied by means of a regular perturbation analysis, where the small parameter is defined as the ratio between the amplitude of variation of the channel wall and the average thickness of the non-wetting phase. Results are valid for Reynolds numbers of the same order of magnitude as that of the expansion parameter. It is thus found that the fluid-fluid interface presents a wavy shape characterized by an amplitude and a phase-shift with respect to the fixed solid-fluid interface. Instabilities of the two-phase flow can arise for large values of the viscosity, flow rate and phase thickness ratios. Results are expected to be a first step towards the understanding of the hydrodynamics of trickle bed reactors, where several flow regimes are possible.     

On Stokes slip flow through a transversely wavy channel

The Stokes flow through a wavy or corrugated channel with surface slip is studied. The correct Navier's partial slip condition is applied and perturbation solutions about the small amplitude to channel width ratio are obtained. As in Stokes slip flow over a sphere, the resistance is not zero even when slip is infinite. The resistance (due to the interaction of waviness and slip) is larger when the corrugations of the two plates are out of phase than that when they are in phase.     

Mass transfer during catalytic reaction in electroosmotically driven flow in a channel microreactor

Analytical solution for concentration profile in a microreactor is obtained during heterogeneous catalytic reaction. Reaction occurs in rectangular microchannel with catalyst-coated walls. Flow is induced electroosmotically in the microchannel. A general solution is obtained for first order reaction using a power series solution. Profiles of conversion, cup-mixing concentration of reactant, etc. and variation of Sherwood number is analyzed as function of operating variables. Analytical solution is compared with numerical results.     

Experiments on the linear instability of flow in a wavy channel

The effects of wall corrugation on the stability of wall-bounded shear flows have been examined experimentally in plane channel flows. One of the channel walls has been modified by introduction of the wavy wall model with the amplitude of 4% of the channel half height and the wave number of 1.02. The experiment is focused on the two-dimensional travelling wave instability and the results are compared with the theory [J.M. Floryan, Two-dimensional instability of flow in a rough channel, Phys. Fluids 17 (2005) 044101 (also: Rept. ESFD-1/2003, Dept. of Mechanical and Materials Engineering, The University of Western Ontario, London, Ontario, Canada, 2003)]. It is shown that the flow is destabilized by the wall corrugation at subcritical Reynolds numbers below 5772, as predicted by the theory. For the present corrugation geometry, the critical Reynolds number is decreased down to about 4000. The spatial growth rates, the disturbance wave numbers and the distribution of disturbance amplitude measured over such wavy wall also agree well with the theoretical results.     

Localised dynamics of laminar pulsatile flow in a rectangular channel

The exploitation of flow pulsation in low-Reynolds number micro/minichannel flows is a potentially useful technique for enhancing cooling of high power photonics and electronics devices. Although the mechanical and thermal problems are inextricably linked, decoupling of the local instantaneous parameters provides insight into underlying mechanisms. The current study performs complementary experimental and analytical analyses to verify novel representations of the pulsating channel flow solutions, which conveniently decompose hydrodynamic parameters into amplitude and phase values relative to a prescribed flow rate, for sinusoidally-pulsating flows of Womersley numbers 1.4 ≤ Wo ≤ 7.0 and a fixed ratio of oscillating flow rate amplitude to steady flow rate equal to 0.9. To the best of the authors’ knowledge, the velocity measurements – taken using particle image velocimetry – constitute the first experimental verification of theory over two dimensions of a rectangular channel. Furthermore, the wall shear stress measurements add to the very limited number of studies that exist for any vessel geometry. The amplification of the modulation component of wall shear stress relative to a steady flow (with flow rate equal to the amplitude of the oscillating flow rate) is an important thermal indicator that may be coupled with future heat transfer measurements. The positive half-cycle time- and space-averaged value is found to increase with frequency owing to growing phase delays and higher amplitudes in the near-wall region of the velocity profiles. Furthermore, the local time-dependent amplification varies depending on the regime of unsteadiness: (i) For quasi-steady flows, the local values are similar during acceleration and deceleration though amplification is greater near the corners over the interval 0 – 0.5π. (ii) At intermediate frequencies, local behaviour begins to differ during accelerating and decelerating periods and the interval of greater wall shear stress near the corners lengthens. (iii) Plug-like flows experience universally high amplifications, with wall shear stress greater near the corners for the majority of the positive half-cycle. The overall fluid mechanical performance of pulsating flow, measured by the ratio of bulk mean wall shear stress and pressure gradient amplifications, is found to reduce from an initial value of 0.97 at Wo = 1.4 to 0.28 at Wo = 7.0, demonstrating the increasing work input required to overcome inertia.     

Elongational effects in the flow of viscoelastic fluid through a wavy channel

The numerical simulation of the viscoelastic flow through a wavy channel was carried out using the modified Giesekus model. It was found that the excess pressure loss relates to the stretch-thickening properties of elongational viscosity and the geometry of the wavy channel through a large elongational component of the flow at the winding part of the channel. The profiles of the axial component of the velocity become significantly asymmetric when the excess pressure loss occurs. Furthermore, the velocity profiles of a 0.1 wt% solution of polyacrylamide were measured using laser Doppler velocimetry. The results of these measurements are compared to the numerical results. Received: 30 June 1998 Accepted: 20 May 1999     

Laminar heat transfer characteristics of internally finned tube with sinusoidal wavy fin

Comparative numerical study of laminar heat transfer characteristics of annular tubes with sinusoidal wavy fins has been conducted both experimentally and numerically with Re = 299–1,475. The uniform heat flux is imposed on the tube outside wall surface. Two tube materials (copper and stainless steel) are considered. It is found that the fluid temperature profile is not linear but convex along the flow direction due to the axial heat conduction in tube wall, and the effects of axial heat conduction on the heat transfer decreases with an increase in Reynolds number or decrease in tube wall thermal conductivity. The axial distributions of local Nusselt number could reach periodically fully developed after 3–5 cycles. The convectional data reduction method based on the traditional method should be improved for tube with high thermal conductivity or low Reynolds numbers, Otherwise, the heat transfer performance of internally finned tube may be underestimated.     

Study of steady viscous flow through a wavy channel: non-orthogonal coordinates

In this paper we consider the steady flow of a viscous fluid through a channel bounded by two sinusoidally varying plates differing in phase by π and separated by a mean distance 2h. For the non-varying channel, the classical parabolic velocity profile for the fully developed flow is well known. An attempt here is made to analyze the flow in a generalized non-orthogonal coordinate system that renders the wavy channels as plane walls. Continuity equation and Navier-Stokes equations are presented in the generalized coordinate system and simplified through use of small perturbation under small Reynolds number approximation. Flow characteristics such as centerline velocity and drag force have been evaluated and discussed.     

On solute dispersion in pulsatile flow through a channel with absorbing walls

The paper presents the longitudinal dispersion of passive tracer molecules released in an incompressible viscous fluid flowing through a channel with reactive walls under the action of a periodic pressure gradient. A finite-difference implicit scheme is adopted to solve the unsteady advection-diffusion equation based on the Aris-Barton method of moments for all time period. Here it is shown how the spreading of tracers is influenced by the shear flow, lateral diffusion about its mean position due to the action of absorption at both the walls. The analysis has been performed for three different cases: steady, periodic and the combined effect of steady and periodic currents, separately. The results show that for all cases the dispersion coefficient asymptotically reaches a stationary state after a certain critical time and it achieves a stationary state at earlier instant of time, when absorption at the walls increases. The axial distributions of mean concentration are determined from the first four central moments by using Hermite polynomial representation for all three different flow velocities.     

Heat transfer enhancement in a channel with a built-in rectangular cylinder

A two dimensional numerical investigation of the unsteady laminar flow pattern and forced convective heat transfer in a channel with a built-in rectangular cylinder is presented. The channel in the entrance region has a length to plate spacing of ten. The computations were made for several Reynolds number and two rectangular cylinder aspect ratios. Hydrodynamic behavior and heat transfer results are obtained by solution of the complete Navier-Stokes and energy equation. The results show that these flow exhibits laminar self-sustained oscillations for Reynolds numbers above the critical one. This study show that oscillatory separated flows result in a significant heat transfer enhancement but also in a significant pressure drop increase.

---

Erhöhung des Wärmeübergangs in einem Spaltkanal mit quer eingebautem Rechteckprisma

Zusammenfassung Es wird eine zweidimensionale numerische Untersuchung des instationären Wärmeübergangs und Druckverlustes im laminar durchströmten Spaltkanal mit quer eingebautem Rechteckprisma dargelegt und zwar für verschiedene Reynoldszahlen und zwei Prismenabmessungen. Als Lösung der Navier-Stokes- und der Energiegleichung resultieren selbsterregt oszillieren de Strömungs- und Temperaturfelder, verbunden mit starker Erhöhung des Wärmeübergangs und des Druckverlustes.

---

List of symbols C _f skin friction coefficient, Eq. (11) - C _D drag coefficient, Eq. (11) - D drag [N/m] - f _app apparent friction factor, Eq. (10) - h cylinder height [m] - H channel height [m] - k thermal conductivity of cylinder [W/mK] - k ₀ thermal conductivity of air [W/mK] - l cylinder length [m] - L channel length [m] - Nu Nusselt number, Eq. (7) - P dimensionless pressure - Pr Prandtl number of air - Re Reynolds number, Eq. (6) - t time [s] - T temperature [K] - T _b bulk temperature [K], Eq. (8) - U, V dimensionless velocity components - X, Y dimensionless coordinates Greek symbols thermal diffusivity [m²/s] - velocity factor, Eq. (11) - dimensionless temperature, Eq. (5) - fluid density [kg/m³] - kinematic viscosity [m²/s] - dimensionless time, Eq. (5) - difference     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.