Analytical solution for transient laminar fully developed free convection in vertical channels

Using Green's function method, analytical solutions for transient fully developed natural convection in open-ended vertical circular and two-parallel-plate channels are presented. Different fundamental boundary conditions for these two configurations have been investigated and the corresponding fundamental solutions are obtained. These fundamental solutions may be used to obtain solutions satisfying more general thermal boundary conditions. In terms of the obtained unsteady temperature and velocity profiles, the transient volumetric flow rate, mixing cup emperature and local nusselt number are estimated.Zusammenfassung Für oben und unten offene vertikale Kanäle mit Kreisquerschnitt bzw. als Parallelplattenanordnung werden unter Verwendung der Methode der Greenschen Funktionen analytische Lösungen für die nichtstationäre, vollausgebildete, natürliche Konvektion gefunden und zwar unter Zugrundelegung verschiedener Fundamental-Randbedingungen bezüglich beider Konfigurationen. Die so ermittelten Fundamentallösungen können zur Gewinnung von Lösungen für allgemeine Randbedingungen dienen. Der zeitlich veränderliche Volumenstrom, die Mischtemperatur und die Nusselt-Zahl werden mit Bezug auf die erhaltenen nichtstationären Profile für Temperatur und Geschwindigkeit näher analysiert.

---

Analytische Lösung für die nichtstationäre vollausgebildete laminare freie Konvektion invertikalen Kanälen

Nomenclature a local heat transfer coefficient based on the area of the heat transfer surface,q/(T _w –T ₀)=±(T/y)_w/(T_w–T₀), minus and plus signs apply respectively for heating and cooling in case of parallel-plate channel and vice versa in case of a tube - average heat transfer coefficient over the channel - c _p specific heat of fluid at constant pressure - f volumetric flow rate, for circular channels and or two-parallel-plate channels - F dimensionless volumetric flow rate,f/(2lvGr^*) for circular channels forfw/(lvGr ^*) for two-parallel-plate channels - g gravitational body force per unit mass (acceleration) - G Green's function - Gr Grashof number,±g(T _w–T₀)w³/v² in case of an isothermal boundary of±gqw ⁴/2kv² in case of a uniform heat flux (UHF) on the heat transfer boundary, the plus and minus signs apply to upward (heating) and downward (cooling) flows, respectively. ThusGr is a positive number in both cases. - Gr ^* modified Grashof number,wGr/l - h heat gained or lost by fluid from the entrance up to a particular elevation in the channel, ₀ fc _p(T _m–T ₀) for all cases - J ₀ Bessel function of zero order - k thermal conductivity of fluid - l height of channel - L dimensionless height of channel,1/Gr ^* - Nu local Nusselt number,|a| w/k - average Nusselt number, - p pressure of fluid inside the channel at any cross-section - p pressure defect at any point,p–p _s - p ₀ pressure of fluid at the channel entrance - p _s hydrostatic pressure, ₀ gz where the minus and plus signs are for upward (heating) and downward (cooling) flows, respectively - p dimensionless pressure defect at any point(pw ⁴)/(₀ l ²² Gr ²) - Pr Prandtl number,c _p/k - q heat flux at the heat transfer surface,q=±k(T/y)_w where the minus and plus signs are, respectively, for cooling and heating in case of circular pipe and vice versa in case of a parallel-plate channel - Ra Rayleigh number,GrPr - Ra ^* modified Rayleigh number,Gr ^*Pr - t time - T fluid temperature at any point - T _m mixing-cup (mixed-mean) temperature over any cross section, for circular channels, and for two-parallelplate channels - T ₀ initial and channel-inlet fluid temperature - T _w temperature of the heat-transfer wall - u axial velocity component at any point - U dimensionless axial velocity,uw ²/(lvGr^*) - w radius of circular tube or width (between plates) of parallel-plate channel - y radial or transverse coordinate - y dimensionless radial or transverse coordinate,y/w - z axial coordinate - Z dimensional axial coordinate,z/(lGr ^*) Greek symbols constant appears in Eq. (8) - parameter appears in Eq. (9) which equals the integration of with respect to or volumetric coefficient of thermal expansion - _n eigenvalues - parameter appears in Eq. (7) - _n eigenvalues - parameter appears in Eq. (12) - _n eigenvalues - parameter appears in Eq. (9) - dimensionless temperature,(T–T ₀)/(T_w–T₀) in case of an isothermal heat transfer boundary and(T–T ₀)/(qw/2k) for UHF boundary - _m dimensionless mixing cup temperature,(T _m–T₀)/(T_w–T₀) in case of an isothermal heat transfer boundary and(T _m–T₀)/(qw/2k) for UHF boundary - _w dimensionless temperature of the heat-transfer wall, equals unity in case of an isothermal heat transfer boundary and(T _w–T₀)/(qw/2k) for a UHF boundary - _n eigenvalues - dynamic viscosity of fluid - kinematic viscosity of fluid, /₀ - fluid density at temperatureT,₀[1–(T–T ₀)] - ₀ fluid density atT ₀ - demensionless time,tk/(cw²)     

Analytical solutions for laminar fully developed flows in double-sine shaped ducts

Fully developed, constant property, laminar flows in double-sine shaped ducts are considered. This cross section represents a limiting inter-plate channel geometry in plate heat exchangers. Accurate analytical solutions based on the Galerkin integral method are presented. Heat transfer with both T and H1 thermal boundary conditions is analyzed; they simulate the most fundamental practical heating/cooling applications. Velocity and temperature distributions, along withfRe, Nu _T , andNu _H1 results for flows in double-sine ducts of different aspect ratios (1/8 ≤γ ≤ 8) are presented. Effects of the relative cross-sectional geometry and thermal boundary conditions are delineated. A comparison of the thermal-hydraulic performance with that of other compact channel geometries is made. The results suggest an optimum (Nu/fRe) performance in a double-sine duct of aspect ratio near unity.     

A variational principle for fully developed laminar heat transfer in uniform channels I

Summary A variational principle for forced convection in laminar flow in uniform channels is formulated. The concepts of thermal potential, dissipation function and thermal force have been introduced and a Lagrangian formulation of the thermal flow field using generalized coordinates is given. Special consideration is given to heat transfer in plug flow when the thermal field is a linear function of the generalized coordinates. The concept of normal and ignorable coordinates is introduced. An example showing the use of Lagrangian equations is worked out.     

Developing and fully developed turbulent flow in ribbed channels

Wall-mounted roughness features, such as ribs, are often placed along the walls of a channel to increase the convective surface area and to augment heat transfer and mixing by increasing turbulence. Depending on the relative roughness size and orientation, the ribs also have varying degrees of increased pressure losses. Designs that use ribs to promote heat transfer encompass the full range of having only a few streamwise ribs, which do not allow fully developed flow conditions, to multiple streamwise ribs, which do allow the flow to become fully developed. The majority of previous studies have focused on perturbing the geometry of the rib with little attention to the spatially and temporally varying flow characteristics and their dependence on the Reynolds number. A staggered rib-roughened channel study was performed using time-resolved digital particle image velocimetry (TRDPIV). Both the developing (entry region) and a fully developed region were interrogated for three Reynolds numbers of 2,500, 10,000, and 20,000. The results indicate that the flow was more sensitive to Reynolds number at the inlet than within the fully developed region. Despite having a similar mean-averaged flowfield structure over the full Reynolds number range investigated, the population and distribution of coherent structures and turbulent dissipation within the fully developed region were also found to be Reynolds number dependent. Exploring the time-accurate flow characteristics revealed that in addition to vortices shed from the rib shear layer, the region of the rib wake was governed by a periodic process of bursting of the wake vortices resulting in the intermittent ejection of the inter-rib recirculation region into the core flow. This periodic process was the driving mechanism resulting in mixing and heat transfer augmentation. A quadrant-splitting burst analysis was also performed to determine the characteristic frequency and duration of inter-rib bursting as well as the wake shedding frequency, both of which were determined to be Reynolds number dependent.     

Experimental analysis of turbulent flow structure in a fully developed rib-roughened rectangular channel with PIV

An experimental study was carried out to investigate the turbulent water-flow structure over one-side micro-repeated ribs in a narrow two-dimensional rectangular channel by particle image velocimetry (PIV). Two rib pitch-to-height ratios p/k of 10 and 20 were investigated while the rib height was held constant at 4 mm. The rib height-to-channel equivalent diameter ratio k/D_e was 0.1. The streamwise mean velocity and turbulent kinetic energy fields in the fully developed flow region of the channel were calculated at three different positions of x=0, 3.8, and -6.2 mm, which corresponded to center, downstream, and upstream of the rib, respectively, and for two Reynolds numbers Re of 7,000 and 20,000. A large-scale turbulent eddy was generated by the rib promoter and then propagated into the mainstream flow, which led to the deformation of the velocity profile. Downstream of the rib, rotating and counter-rotating eddies were also generated by the rib promoter. The enhancement of the turbulent kinetic energy was not changed when the Reynolds number increased from 7,000 to 20,000 between p/k=20 and 10. The reattachment length L_R was measured from velocity vector fields in the developing, fully developed, and exit regions of the flow over the range Re=1,400-50,000. The results showed that the ratio p/k and the Reynolds number had no significant effect on the reattachment length beyond a critical value of Re=15,000, where L_R was found to be approximately 4 times the rib height.     

Study on transient local entropy generation in pulsating fully developed laminar flow through an externally heated pipe

This study presents the investigation of transient local entropy generation rate in pulsating fully developed laminar flow through an externally heated pipe. The flow inlet to the pipe is considered as pulsating at a constant period and amplitude (only the velocity oscillates). The simulations are extended to include different pulsating flow cases (sinusoidal flow, step flow, and saw-down flow). To determine the effects of the mean velocity, the period and the amplitude of the pulsating flow on the entropy generation rate, the pulsating flow is examined for various cases of these parameters. Two-dimensional flow and temperature fields are computed numerically with the help of the fluent computational fluid dynamics (CFD) code. In addition to this CFD code, a computer program has been developed to calculate numerically the entropy generation and other thermodynamic parameters by using the results of the calculations performed for the flow and temperature fields. In all investigated cases, the irreversibility due to the heat transfer dominates. The step flow constitutes the highest temperature (about 919 K) and generates the highest total entropy rate (about 0.033 W/K) within the pipe. The results of this study indicate that in the considered situations, the inverse of square of temperature (1/T ²) is more dominant on the entropy generation than the temperature gradients, and that the increase of the mean velocity of the pulsating flow has an adverse effect on the ratio of the useful energy transfer rate to irreversibility rate.     

Bifurcation of developed flow in rectangular channels rotating about the transverse axis

The linear problem of the stability of Poiseuille flow between two infinite plates rotating about an axis parallel to the plates and normal to the flow direction was studied in [1, 2]. It was established that the flow is least stable with respect to disturbances in the form of standing waves known as Taylor eddies. The experimental data of [1, 3] and the results [4] of a numerical integration of the Navier-Stokes equations for channels with cross sections highly elongated in the direction of the axis of rotation are in good agreement with the conclusions of the linear theory. In the case of channels with a side ratio of the order of unity, which is of greater practical interest, the primary flow becomes essentially three-dimensional and evolves with variation of the governing criteria: the Reynolds and Rossby numbers. Obviously, this seriously complicates the use of the methods of the linear theory. The effect of the ratio of the sides of the cross section on the stability of the primary flow regime was studied experimentally in [3]. The present article describes the results of an investigation of the problem based on a numerical method of integrating the nonlinear Navier-Stokes equations. Moreover, an asymptotic estimate of the stability limit of the primary regime, based on a local condition of inviscid instability of rotating flows, is presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 27–32, September–October, 1985.     

Flow visualization and LDV measurement of fully developed laminar flow in helically coiled tubes

Secondary flow structure in helically coiled tubes is characterized by laser light-sheet flow visualization photographs and laser Doppler velocimetry measured velocity vector field. The torsion-to-curvature ratio and Reynolds number, based on the tube diameter and bulk average axial velocity, were varied from 0.06 to 5.55 and from 35 to 330, respectively, to study their effects on the secondary flow patterns. Good agreement was found between the results obtained in the present work and those predicted previously. That the dominance of the torsion effect on the secondary flow is not limited to low Reynolds numbers as reported previously is pointed out. Moreover, the transformation of two recirculating cells into one cell is documented in detail.     

Heat transfer to non-Newtonian fluids in laminar flow through rectangular ducts

Heat transfer to non-newtonian fluids flowing laminarly through rectangular ducts is examined. The conservation equations of mass, momentum, and energy are solved numerically with the aid of a finite volume technique. The viscoelastic behavior of the fluid is represented by the Criminale-Ericksen-Filbey (CEF) constitutive equation. Secondary flows occur due to the elastic behavior of the fluid, and, consequently, heat transfer is strongly enhanced. It is observed that shear thinning yields negligible heat transfer enhancement effect, when compared with the secondary flow effect. Maximum heat transfer is shown to occur for some combinations of parameters. Thus, there are optimal combinations of aspect ratio and Reynolds numbers, which depend on the fluid's mechanical behavior. This result can be usefully explored in thermal designs of certain industrial processes.     

A finite element analysis of laminar unsteady flow of viscoelastic fluids through channels with non-uniform cross-sections

The effects of non-Newtonian behaviour of a fluid and unsteadiness on flow in a channel with non-uniform cross-section have been investigated. The rheological behaviour of the fluid is assumed to be described by the constitutive equation of a viscoelastic fluid obeying the Oldroyd-B model. The finite element method is used to analyse the flow. The novel features of the present method are the adoption of the velocity correction technique for the momentum equations and of the two-step explicit scheme for the extra stress equations. This approach makes the computational scheme simple in algorithmic structure, which therefore implies that the present technique is capable of handling large-scale problems. The scheme is completed by the introduction of balancing tensor diffusivity (wherever necessary) in the momentum equations. It is important to mention that the proper boundary condition for pressure (at the outlet) has been developed to solve the pressure Poisson equation, and then the results for velocity, pressure and extra stress fields have been computed for different values of the Weissenberg number, viscosity due to elasticity, etc. Finally, it is pertinent to point out that the present numerical scheme, along with the proper boundary condition for pressure developed here, demonstrates its versatility and suitability for analysing the unsteady flow of viscoelastic fluid through a channel with non-uniform cross-section.     

Analytical series solution for the fully developed forced convection duct flow with frictional heating and variable viscosity

Laminar forced convection flow of a liquid in the fully developed region of a circular duct with isothermal wall is analyzed. The effects of viscous dissipation as well as of temperature dependent viscosity are taken into account. The coupled momentum and energy equations are solved analytically by means of a power series method. Then, reference is made to the Poiseuille model for the temperature change of viscosity. For a fixed value of the axial pressure gradient along the duct, dual solutions are found for the velocity and temperature fields. Although dual solutions correspond to the same value of the axial pressure gradient, they lead in general to different values of the average fluid velocity, of the average fluid temperature and of the wall heat flux. It is shown that, for a given fluid and for a fixed duct radius, the absolute value of the axial pressure gradient has an upper bound above which no steady laminar solution can exist.     

A multigrid method for predicting periodically fully developed flow

The use of multigrid methods in complex fluid flow problems is still under development. In this paper a full multigrid procedure has been incorporated in a finite volume solution for predicting fully developed fluid flow in a streamwise periodic geometry. Steady computations in two-dimensional body fitted co-ordinates have shown considerable savings in computation time by this multigrid method.     

Radiation interactions in transient energy transfer in fully developed laminar flows

Analysis and numerical procedures are presented to investigate the transient radiative interactions of nongray absorbing-emitting species in laminar fully-developed flows between two parallel black plates. The particular species considered are OH, CO, CO₂, and H₂O and different mixtures of these species. Transient and steady-state results are obtained for the temperature distribution and bulk temperature for different plate spacings, wall temperatures, and pressures. Results, in general, indicate that the rate of radiative heating can be quite high during earlier times. This information is useful in designing thermal protection systems for transient operations.Nomenclature A band absorptance = A(u, ), cm^–1 - A ₀ band width parameter, cm^–1 - C ₀ correlation parameter, atm^–1-cm^–1 - C _p specific heat at constant pressure, kJ/kg-K = erg/gm-K - e Planck's function, (W-cm^–2)/cm^–1 - e ₀ Planck's function evaluated at wave number ₀ - e ₁, e ₂ emissive power of surfaces with temperatures T ₁ and T ₂, W-cm^–2 - H _1i , H ₁ gas property for the large path length limit - k thermal conductivity, erg/cm-s-K - K ₁ gas property for the optically thin limit - M large path length parameter, nondimensional - N optically thin parameter, nondimensional - P pressure, atm - q _R total radiative heat flux, W/cm² - q _Rw spectral radiation heat flux, (W-cm^–2)/cm^–1 - q _w wall het flux, W/cm² - S integrated intensity of a wide band, atm^–1-cm² - t time, s - T temperature, K - T ₁ wall temperature, K; T ₁=T _w - T _b bulk temperature, K - u nondimensional coordinate = SPy/A ₀ - u ₀ nondimensional path length = SPL/A ₀ - y transverse coordinate, cm - line structures parameter - nondimensional temperature - _b dimensionless bulk temperature - spectral absorption coefficient, cm^–1 - nondimensional coordinate = y/L=u/u ₀ - density, kg/m³ - nondimensional time - wave number, cm^–1 - ₀ wave number at the band center, cm^–1     

Method of moments for laminar dispersion in an oscillatory flow through curved channels with absorbing walls

The unsteady dispersion of a solute, when the fluid is driven through a curved channel with absorbing walls by an imposed pulsatile pressure gradient, is studied using the method of moments. The study examines the effect of oscillatory Reynolds number, amplitude/frequency of the pressure pulsation and boundary absorption on the longitudinal dispersion. The methodology involves a set of unsteady integral moment equations obtained by applying the Aris-Barton method of moments on the convective-diffusion equation for a curved channel. Central moments are obtained from the moment equations which are solved by a finite-difference implicit scheme. The effect of curvature and boundary absorption on the effective dispersion coefficient from the initial to the stationary stage of the oscillatory flow is studied. Amplitude of the effective dispersion coefficient is found to increase with curvature and decrease with frequency of the pressure pulsation. For large Peclet number and Schmidt number, the amplitude of the dispersion coefficient can be 1.6 times that in a straight channel at large times. Also, for large times, the amplitude of the dispersion coefficient is twice the amplitude of the dispersion coefficient as α, the frequency parameter changes from 0.5 to 1.0. The axial distributions of mean concentration are determined from the first four central moments by using the Hermite polynomial representation. The effect of curvature is to delay the stationary state and also the approach to normality of the concentration distribution. The study has importance in understanding the spreading of pollutants in tidal basins and natural current fields.     

Lyapunov analysis for fully developed homogeneous isotropic turbulence

The present work studies the isotropic and homogeneous turbulence for incompressible fluids through a specific Lyapunov analysis. The analysis consists in the calculation of the velocity fluctuation through the Lyapunov theory applied to the local deformation using the Navier-Stokes equations, and in the study of the mechanism of energy cascade through the finite scale Lyapunov analysis of the relative motion between two particles. The analysis provides an explanation for the mechanism of energy cascade, leads to the closure of the von Kármán-Howarth equation, and describes the statistics of the velocity difference. Several tests and numerical results are presented.     

Entropy generation analysis of fully developed laminar forced convection in a helical tube with uniform wall temperature

In the present study, fully developed laminar flow and heat transfer in a helically coiled tube with uniform wall temperature have been investigated analytically. Expressions involving relevant variables for entropy generation rate contributed to heat transfer and friction loss, and total entropy generation rate have been derived. The effect of various flow and coil parameters like Reynolds number, curvature ratio, coil pitch, etc. on the entropy generation rate has been studied for two fluids- air and water. The results of the present study have been compared to the corresponding entropy generation values of straight pipe. Investigating the results, some optimum values for Reynolds number have been proposed and compared with the optimum Reynolds numbers of laminar flow inside a coiled tube subjected to constant heat flux boundary condition.     

Electroviscous effects in steady fully developed flow of a power-law liquid through a cylindrical microchannel

Electroviscous effects in steady, fully developed, pressure-driven flow of power-law liquids through a uniform cylindrical microchannel have been investigated numerically by solving the Poisson–Boltzmann and the momentum equations using a finite difference method. The pipe wall is considered to have uniform surface charge density and the liquid is assumed to be a symmetric 1:1 electrolyte solution. Electroviscous resistance reduces the velocity adjacent to the wall, relative to the velocity on the axis. The effect is shown to be greater when the liquid is shear-thinning, and less when it is shear-thickening, than it is for Newtonian flow. For overlapping electrical double layers and elevated surface charge density, the electroviscous reduction in the near-wall velocity can form an almost stationary (zero shear) layer there when the liquid is shear-thinning. In that case, the liquid behaves approximately as if it is flowing through a channel of reduced diameter. The induced axial electrical field shows only a weak dependence on the power-law index with the dependence being greatest for shear-thinning liquids. This field exhibits a local maximum as surface charge density increases from zero, even though the corresponding electrokinetic resistance increases monotonically. The magnitude of the electroviscous effect on the apparent viscosity, as measured by the ratio of the apparent and physical consistency indices, decreases monotonically as the power-law index increases. Thus, overall, the electroviscous effect is stronger in shear-thinning, and weaker in shear-thickening liquids, than it is when the liquid is Newtonian.     

Numerical investigation of a laminar pulsating flow in a rectangular duct

A pulsating laminar flow of a viscous, incompressible liquid in a rectangular duct has been studied. The motion is induced under an imposed pulsating pressure difference. The problem is solved numerically. Different flow regimes are characterized by a non‐dimensional parameter based on the frequency (ω) of the imposed pressure gradient oscillations and the width of the duct (h). This, in fact, is the Reynolds number of the problem at hand. The induced velocity has a phase lag (shift) with respect to the imposed pressure oscillations, which varies from zero at very slow oscillations, to 90° at fast oscillations. The influence of the aspect ratio of the rectangular duct and the pulsating pressure gradient frequency on the phase lag, the amplitude of the induced oscillating velocity, and the wall shear were analyzed. Copyright © 1999 John Wiley & Sons, Ltd.     

Analytical study of mixed electroosmotic-pressure-driven flow in rectangular micro-channels

Operational state of many miniaturized devices deals with flow field in microchannels. Pressure-driven flow (PDF) and electroosmotic flow (EOF) can be recognized as the two most important types of the flow field in such channels. EOF has many advantages in comparison with PDF, such as being vibration free and not requiring any external mechanical pumps or moving parts. However, the disadvantages of this type of flow such as Joule heating, electrophoresis demixing, and not being suitable for mobile devices must be taken into consideration carefully. By using mixed electroosmotic/pressure-driven flow, the role of EOF in producing desired velocity profile will be reduced. In this way, the advantages of EOF can be exploited, and its disadvantages can be prevented. Induced pressure gradient can be utilized in order to control the separation in the system. Furthermore, in many complicated geometries such as T-shape microchannels, turns may induce pressure gradient to the electroosmotic velocity. While analytical formulas are completely essential for analysis and control of any industrial and laboratory microdevices, lack of such formulas in the literature for solving Poisson–Boltzmann equation and predicting electroosmotic velocity field in rectangular domains is evident. In the present study, first a novel method is proposed to solve Poisson–Boltzmann equation (PBE). Subsequently, this solution is utilized to find the electroosmotic and the mixed electroosmotic/pressure-driven velocity profile in a rectangular domain of the microchannels. To demonstrate the accuracy of the presented analytical method in solving PBE and finding electroosmotic velocity, a general nondimensional example is analyzed, and the results are compared with the solution of boundary element method. Additionally, the effects of different nondimensional parameters and also aspect ratio of channels on the electroosmotic part of the flow field will be investigated.