Numerical study on mixed convection in a horizontal flow past a square porous cavity using UNIFAES scheme

Mixed convection on the flow past a heated length and past a porous cavity located in a horizontal wall bounding a saturated porous medium is numerically simulated. The cavity is heated from below. The steady-state regime is studied for several intensities of the buoyancy effects due to temperature variations. The influences of Péclet and Rayleigh numbers on the flow pattern and the temperature distributions are examined. Local and global Nusselt numbers are reported for the heated surface. The convective-diffusive fluxes at the volume boundaries are represented using the UNIFAES, Unified Finite Approach Exponential-type Scheme, with the Power-Law approximation to reduce the computing time. The conditions established by Rivas for the quadratic order of accuracy of the central differencing to be maintained in irregular grids are shown to be extensible to other quadratic schemes, including UNIFAES, so that accuracy estimates could be obtained.

Laminar forced convection to a pseudoplastic fluid food in circular and annular ducts

In this manuscript we investigated experimentally the steady-state heat transfer to an important pseudoplastic fluid food, the soursop juice, flowing in laminar regime through circular and concentric annular ducts. The mean convection heat transfer coefficients, determined by measuring the bulk temperatures before and after the heating sections with constant temperatures of the tube walls, were used to correlate simple new empiric expressions to estimate the average Nusselt number in the thermal entrance of the considered geometries. In addition, the thermophysical properties of the tested fluid food, as well as the rheological behavior, being essential for the heat transfer analyses, were experimentally determined. (c) 2006 Elsevier Ltd. All rights reserved.

Nonlinear parameter estimation in laminar forced convection within a circular sector tube

In this article we examine an inverse heat convection problem of estimating unknown parameters of a parameterized variable boundary heat flux. The physical problem is a hydrodynamically developed, thermally developing, three-dimensional steady state laminar flow of a Newtonian fluid inside a circular sector duct, insulated in the flat walls and subject to unknown wall heat flux at the curved wall. Results are presented for polynomial and sinusoidal trial functions, and the unknown parameters as well as surface heat fluxes are determined. Depending on the nature of the flow, on the position of experimental points the inverse problem sometimes could not be solved. Therefore, an identification condition is defined to specify a condition under which the inverse problem can be solved. Once the parameters have been computed it is possible to obtain the statistical significance of the inverse problem solution. Therefore, approximate confidence bounds based on standard statistical linear procedure, for the estimated parameters, are analyzed and presented.

Magnetohydrodynamic Effects On Mixed Convection Flows In Channels And Ducts

This work focuses on magnetohydrodynamic (MHD) mixed convection flow of electrically conducting fluids enclosed in simple 1D and 2D geometries in steady periodic regime. In particular, in Chapter one a short overview is given about the history of MHD, with reference to papers available in literature, and a listing of some of its most common technological applications, whereas Chapter two deals with the analytical formulation of the MHD problem, starting from the fluid dynamic and energy equations and adding the effects of an external imposed magnetic field using the Ohm's law and the definition of the Lorentz force. Moreover a description of the various kinds of boundary conditions is given, with particular emphasis given to their practical realization. Chapter three, four and five describe the solution procedure of mixed convective flows with MHD effects. In all cases a uniform parallel magnetic field is supposed to be present in the whole fluid domain transverse with respect to the velocity field. The steady-periodic regime will be analyzed, where the periodicity is induced by wall temperature boundary conditions, which vary in time with a sinusoidal law. Local balance equations of momentum, energy and charge will be solved analytically and numerically using as parameters either geometrical ratios or material properties. In particular, in Chapter three the solution method for the mixed convective flow in a 1D vertical parallel channel with MHD effects is illustrated. The influence of a transverse magnetic field will be studied in the steady periodic regime induced by an oscillating wall temperature. Analytical and numerical solutions will be provided in terms of velocity and temperature profiles, wall friction factors and average heat fluxes for several values of the governing parameters. In Chapter four the 2D problem of the mixed convective flow in a vertical round pipe with MHD effects is analyzed. Again, a transverse magnetic field influences the steady periodic regime induced by the oscillating wall temperature of the wall. A numerical solution is presented, obtained using a finite element approach, and as a result velocity and temperature profiles, wall friction factors and average heat fluxes are derived for several values of the Hartmann and Prandtl numbers. In Chapter five the 2D problem of the mixed convective flow in a vertical rectangular duct with MHD effects is discussed. As seen in the previous chapters, a transverse magnetic field influences the steady periodic regime induced by the oscillating wall temperature of the four walls. The numerical solution obtained using a finite element approach is presented, and a collection of results, including velocity and temperature profiles, wall friction factors and average heat fluxes, is provided for several values of, among other parameters, the duct aspect ratio. A comparison with analytical solutions is also provided, as a proof of the validity of the numerical method. Chapter six is the concluding chapter, where some reflections on the MHD effects on mixed convection flow will be made, in agreement with the experience and the results gathered in the analyses presented in the previous chapters. In the appendices special auxiliary functions and FORTRAN program listings are reported, to support the formulations used in the solution chapters.

Convecção mista em cavidades com fontes de calor aletadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Heat transfer in laminar flow of non-Newtonian fluids in ducts of elliptical section

Laminar forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumption used in this work is a laminar flow of a power flow inside elliptical tube, under a boundary condition of first kind with constant physical properties and negligible axial heat diffusion (high Peclet number). To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number and the average Nusselt number for various cross-section aspect ratios. (C) 2006 Elsevier. SAS. All rights reserved.

LAMINAR THERMALLY DEVELOPING FLOW INSIDE RIGHT-ANGULARLY TRIANGULAR DUCTS

An analytical approach based on the generalized integral transform technique is presented, for the solution of laminar forced convection within the thermal entry region of ducts with arbitrarily shaped cross-sections. The analysis is illustrated through consideration of a right triangular duct subjected to constant wall temperature boundary condition. Critical comparisons are made with results available in the literature, from direct numerical approaches. Numerical results for dimensionless average temperature and Nusselt numbers are presented for different apex angles.

Thermally developing forced convection of non-Newtonian fluids inside elliptical ducts

Laminar-forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumptions used in this work are a non-Newtonian fluid, laminar flow, constant physical properties, and negligible axial heat diffusion (high Peclet number). Most of the previous research in elliptical ducts deal mainly with aspects of fully developed laminar flow forced convection, such as velocity profile, maximum velocity, pressure drop, and heat transfer quantities. In this work, we examine heat transfer in a hydrodynamically developed, thermally developing laminar forced convection flow of fluid inside an elliptical tube under a second kind of a boundary condition. To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform, where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number, and the average Nusselt number for various cross-section aspect ratios.

Integral transform solution for the forced convection of herschel-bulkley fluids in circular tubes and parallel-plates ducts

ABSTRACT: The thermal entry region in laminar forced convection of Herschel-Bulkley fluids is solved analytically through the integral transform technique, for both circular and parallel-plates ducts, which are maintained at a prescribed wall temperature or at a prescribed wall heat flux. The local Nusselt numbers are obtained with high accuracy in both developing and fully-developed thermal regions, and critical comparisons with previously reported numerical results are performed.

Simulação numérica da convecção mista em cavidade preenchida com meio poroso heterogêneo e homogêneo

In this work is presented mixed convection heat transfer inside a lid-driven cavity heated from below and filled with heterogeneous and homogeneous porous medium. In the heterogeneous approach, the solid domain is represented by heat conductive equally spaced blocks; the fluid phase surrounds the blocks being limited by the cavity walls. The homogeneous or pore-continuum approach is characterized by the cavity porosity and permeability. Generalized mass, momentum and energy conservation equations are obtained in dimensionless form to represent both the continuum and the pore-continuum models. The numerical solution is obtained via the finite volume method. QUICK interpolation scheme is set for numerical treatment of the advection terms and SIMPLE algorithm is applied for pressure-velocity coupling. Aiming the laminar regime, the flow parameters are kept in the range of 102≤Re≤103 and 103≤Ra≤106 for both the heterogeneous and homogeneous approaches. In the tested configurations for the continuous model, 9, 16, 36, and 64 blocks are considered for each combination of Re and Ra being the microscopic porosity set as constant φ=0,64 . For the pore-continuum model the Darcy number (Da) is set according to the number of blocks in the heterogeneous cavity and the φ. Numerical results of the comparative study between the microscopic and macroscopic approaches are presented. As a result, average Nusselt number equations for the continuum and the pore continuum models as a function of Ra and Re are obtained.

Aplicação da técnica de elementos finitos na estratificação de armazenadores térmicos com barreiras de meios porosos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Convecção natural em cavidades triangulares

Pós-graduação em Engenharia Mecânica - FEG

Heat transfer performance investigation of nanofluids flow in pipe

Different types of base fluids, such as water, engine oil, kerosene, ethanol, methanol, ethylene glycol etc. are usually used to increase the heat transfer performance in many engineering applications. But these conventional heat transfer fluids have often several limitations. One of those major limitations is that the thermal conductivity of each of these base fluids is very low and this results a lower heat transfer rate in thermal engineering systems. Such limitation also affects the performance of different equipments used in different heat transfer process industries. To overcome such an important drawback, researchers over the years have considered a new generation heat transfer fluid, simply known as nanofluid with higher thermal conductivity. This new generation heat transfer fluid is a mixture of nanometre-size particles and different base fluids. Different researchers suggest that adding spherical or cylindrical shape of uniform/non-uniform nanoparticles into a base fluid can remarkably increase the thermal conductivity of nanofluid. Such augmentation of thermal conductivity could play a more significant role in enhancing the heat transfer rate than that of the base fluid. Nanoparticles diameters used in nanofluid are usually considered to be less than or equal to 100 nm and the nanoparticles concentration usually varies from 5% to 10%. Different researchers mentioned that the smaller nanoparticles concentration with size diameter of 100 nm could enhance the heat transfer rate more significantly compared to that of base fluids. But it is not obvious what effect it will have on the heat transfer performance when nanofluids contain small size nanoparticles of less than 100 nm with different concentrations. Besides, the effect of static and moving nanoparticles on the heat transfer of nanofluid is not known too. The idea of moving nanoparticles brings the effect of Brownian motion of nanoparticles on the heat transfer. The aim of this work is, therefore, to investigate the heat transfer performance of nanofluid using a combination of smaller size of nanoparticles with different concentrations considering the Brownian motion of nanoparticles. A horizontal pipe has been considered as a physical system within which the above mentioned nanofluid performances are investigated under transition to turbulent flow conditions. Three different types of numerical models, such as single phase model, Eulerian-Eulerian multi-phase mixture model and Eulerian-Lagrangian discrete phase model have been used while investigating the performance of nanofluids. The most commonly used model is single phase model which is based on the assumption that nanofluids behave like a conventional fluid. The other two models are used when the interaction between solid and fluid particles is considered. However, two different phases, such as fluid and solid phases is also considered in the Eulerian-Eulerian multi-phase mixture model. Thus, these phases create a fluid-solid mixture. But, two phases in the Eulerian-Lagrangian discrete phase model are independent. One of them is a solid phase and the other one is a fluid phase. In addition, RANS (Reynolds Average Navier Stokes) based Standard κ-ω and SST κ-ω transitional models have been used for the simulation of transitional flow. While the RANS based Standard κ-ϵ, Realizable κ-ϵ and RNG κ-ϵ turbulent models are used for the simulation of turbulent flow. Hydrodynamic as well as temperature behaviour of transition to turbulent flows of nanofluids through the horizontal pipe is studied under a uniform heat flux boundary condition applied to the wall with temperature dependent thermo-physical properties for both water and nanofluids. Numerical results characterising the performances of velocity and temperature fields are presented in terms of velocity and temperature contours, turbulent kinetic energy contours, surface temperature, local and average Nusselt numbers, Darcy friction factor, thermal performance factor and total entropy generation. New correlations are also proposed for the calculation of average Nusselt number for both the single and multi-phase models. Result reveals that the combination of small size of nanoparticles and higher nanoparticles concentrations with the Brownian motion of nanoparticles shows higher heat transfer enhancement and thermal performance factor than those of water. Literature suggests that the use of nanofluids flow in an inclined pipe at transition to turbulent regimes has been ignored despite its significance in real-life applications. Therefore, a particular investigation has been carried out in this thesis with a view to understand the heat transfer behaviour and performance of an inclined pipe under transition flow condition. It is found that the heat transfer rate decreases with the increase of a pipe inclination angle. Also, a higher heat transfer rate is found for a horizontal pipe under forced convection than that of an inclined pipe under mixed convection.

Development of a numerical platform for the modeling and optimal control of liquid metal flows

The main purpose of this work is to develop a numerical platform for the turbulence modeling and optimal control of liquid metal flows. Thanks to their interesting thermal properties, liquid metals are widely studied as coolants for heat transfer applications in the nuclear context. However, due to their low Prandtl numbers, the standard turbulence models commonly used for coolants as air or water are inadequate. Advanced turbulence models able to capture the anisotropy in the flow and heat transfer are then necessary. In this thesis, a new anisotropic four-parameter turbulence model is presented and validated. The proposed model is based on explicit algebraic models and solves four additional transport equations for dynamical and thermal turbulent variables. For the validation of the model, several flow configurations are considered for different Reynolds and Prandtl numbers, namely fully developed flows in a plane channel and cylindrical pipe, and forced and mixed convection in a backward-facing step geometry. Since buoyancy effects cannot be neglected in liquid metals-cooled fast reactors, the second aim of this work is to provide mathematical and numerical tools for the simulation and optimization of liquid metals in mixed and natural convection. Optimal control problems for turbulent buoyant flows are studied and analyzed with the Lagrange multipliers method. Numerical algorithms for optimal control problems are integrated into the numerical platform and several simulations are performed to show the robustness, consistency, and feasibility of the method.

Forced convection to laminar flow of liquid egg yolk in circular and annular ducts

The steady-state heat transfer in laminar flow of liquid egg yolk - an important pseudoplastic fluid food - in circular and concentric annular ducts was experimentally investigated. The average convection heat transfer coefficients, determined by measuring temperatures before and after heating sections with constant temperatures at the tube wall, were used to obtain simple new empirical expressions to estimate the Nusselt numbers for fully established flows at the thermal entrance of the considered geometries. The comparisons with existing correlations for Newtonian and non-Newtonian fluids resulted in excellent agreement. The main contribution of this work is to supply practical and easily applicable correlations, which are, especially for the case of annulus, rather scarce and extensively required in the design of heat transfer operations dealing with similar shear-thinning products. In addition, the experimental results may support existing theoretical analyses.