Numerical Computation of Fluid Flow and Heat Transfer in Microchannels

2005 ◽  
Author(s):  
Ningli Liu ◽  
Rene Chevray ◽  
Gerald A. Domoto ◽  
Elias Panides

A finite difference numerical approach for solving slightly compressible, time-dependent, viscous laminar flow is presented in this study. Simplified system of Navier-Stokes equations and energy equation are employed in the study in order to perform more efficient numerical calculations. Fluid flow and heat transfer phenomena in two dimensional microchannels are illustrated numerically in this paper. This numerical approach provides a complete numerical simulation of the development of the fluid flow and the temperature profiles through multi-dimensional microchannels.

2006 ◽  
Vol 129 (1) ◽  
pp. 106-115 ◽  
Author(s):  
A. B. Rahimi ◽  
R. Saleh

The unsteady viscous flow and heat transfer in the vicinity of an axisymmetric stagnation point of an infinite rotating circular cylinder with transpiration U0 are investigated when the angular velocity and wall temperature or wall heat flux all vary arbitrarily with time. The free stream is steady and with a strain rate of Γ. An exact solution of the Navier-Stokes equations and energy equation is derived in this problem. A reduction of these equations is obtained by the use of appropriate transformations for the most general case when the transpiration rate is also time-dependent but results are presented only for uniform values of this quantity. The general self-similar solution is obtained when the angular velocity of the cylinder and its wall temperature or its wall heat flux vary as specified time-dependent functions. In particular, the cylinder may rotate with constant speed, with exponentially increasing/decreasing angular velocity, with harmonically varying rotation speed, or with accelerating/decelerating oscillatory angular speed. For self-similar flow, the surface temperature or its surface heat flux must have the same types of behavior as the cylinder motion. For completeness, sample semi-similar solutions of the unsteady Navier-Stokes equations have been obtained numerically using a finite-difference scheme. Some of these solutions are presented for special cases when the time-dependent rotation velocity of the cylinder is, for example, a step-function. All the solutions above are presented for Reynolds numbers, Re=Γa2∕2υ, ranging from 0.1 to 1000 for different values of Prandtl number and for selected values of dimensionless transpiration rate, S=U0∕Γa, where a is cylinder radius and υ is kinematic viscosity of the fluid. Dimensionless shear stresses corresponding to all the cases increase with the increase of Reynolds number and suction rate. The maximum value of the shear stress increases with increasing oscillation frequency and amplitude. An interesting result is obtained in which a cylinder rotating with certain exponential angular velocity function and at particular value of Reynolds number is azimuthally stress-free. Heat transfer is independent of cylinder rotation and its coefficient increases with the increasing suction rate, Reynolds number, and Prandtl number. Interesting means of cooling and heating processes of cylinder surface are obtained using different rates of transpiration.


2003 ◽  
Author(s):  
Tien-Chien Jen ◽  
Tuan-Zhou Yan ◽  
S. H. Chan

A three-dimensional computational model is developed to analyze fluid flow in a semi-porous channel. In order to understand the developing fluid flow and heat transfer process inside the semi-porous channels, the conventional Navier-Stokes equations for gas channel, and volume-averaged Navier-Stokes equations for porous media layer are adopted individually in this study. Conservation of mass, momentum and energy equations are solved numerically in a coupled gas and porous media domain in a channel using the vorticity-velocity method with power law scheme. Detailed development of axial velocity, secondary flow and temperature fields at various axial positions in the entrance region are presented. The friction factor and Nusselt number are presented as a function of axial position, and the effects of the size of porous media inside semi-porous channel are also analyzed in the present study.


Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 756
Author(s):  
Federico Lluesma-Rodríguez ◽  
Francisco Álcantara-Ávila ◽  
María Jezabel Pérez-Quiles ◽  
Sergio Hoyas

One numerical method was designed to solve the time-dependent, three-dimensional, incompressible Navier–Stokes equations in turbulent thermal channel flows. Its originality lies in the use of several well-known methods to discretize the problem and its parallel nature. Vorticy-Laplacian of velocity formulation has been used, so pressure has been removed from the system. Heat is modeled as a passive scalar. Any other quantity modeled as passive scalar can be very easily studied, including several of them at the same time. These methods have been successfully used for extensive direct numerical simulations of passive thermal flow for several boundary conditions.


1978 ◽  
Vol 100 (4) ◽  
pp. 565-571 ◽  
Author(s):  
B. E. Launder ◽  
T. H. Massey

A scheme for handling the numerical analysis of viscous flow and heat transfer in tube banks is presented. It involves the use of a cylindrical network of nodes in the vicinity of the tubes with a Cartesian mesh covering the remainder of the flow domain. The approach has been incorporated into the numerical solving algorithm for the Navier Stokes equations of Gasman, et al. [8]. A number of demonstration calculations is presented including a numerical simulation of the staggered square bank for which Bergelin and co-workers [4, 9] have reported experimental results for pressure drop and heat transfer rate. Agreement between predicted and measured characteristics is satisfactory when account is taken of end and entry effects that are present in the experiments but necessarily omitted from the calculations. Indeed the close agreement of the laminar predictions with measurements extends to Reynolds numbers in excess of 1000, a level at which it has hitherto been supposed that turbulent motion in the fluid made a substantial contribution to friction and heat transfer.


Author(s):  
A. D. D. Craik

Classes of exact solutions of the Navier–Stokes equations for incompressible fluid flow are explored. These have spatially-uniform velocity gradients at each instant, but often display complex temporal behaviour. Particular illustrative cases are described and related to previously-known solutions.


2021 ◽  
Vol 406 ◽  
pp. 66-77
Author(s):  
Abdel Halim Zitouni ◽  
Pierre Spiteri ◽  
Mouloud Aissani ◽  
Younes Benkheda

In this work, the heat transfer by conduction and convection mode and effect of fluid flow on the morphology of the weld pool and the welding properties is investigated during Tungsten Inert Gas (TIG) process. In the first part, a computation code under Fortran was elaborated to solve the equations resulting from the finite difference discretization of the heat equation, taking into account the liquid-solid phase change with the associated boundary conditions. In order to calculate the velocity field during welding, the Navier-Stokes equations in the melt zone were simplified and solved considering their stream-vorticity formulation. A mathematical model was developed to study the effect of the melted liquid movement on the weld pool. The evolution of the fraction volume of the liquid and the thermal fields promoted the determination of the molten zone (MZ) and the Heat Affected Zone (HAT) dimensions, which seems to be in good agreement with literature.


Author(s):  
V. Dakshina Murty

A numerical method based on the finite elements is applied to the cooling of pulse detonation tube using heat pipe technology. Towards this end, the fluid flow and heat transfer in the wick are modeled as flow in a porous medium. The flow is described using the so called Darcy Brinkman model which has close resemblance to the Navier-Stokes equations. It is found that for Darcy numbers less than 0.0001 the results are indistinguishable from regular Darcy flows. The shape of the heat pipe is that of a fin with the proportion of the length of the evaporator section being varied. In this study two values of this ratio have been used, namely 1 and 0.5.


1990 ◽  
Author(s):  
B. L. Lapworth ◽  
J. W. Chew

Numerical solutions of the Reynolds-averaged Navier-Stokes equations have been used to model the influence of cobs and a bolt cover on the flow and heat transfer in a rotating cavity with an imposed radial outflow of air. Axisymmetric turbulent flow is assumed using a mixing length turbulence model. Calculations for the non-plane discs are compared with plane disc calculations and also with the available experimental data. The calculated flow structures show good agreement with the experimentally observed trends. For the cobbed and plane discs, Nusselt numbers are calculated for a combination of flow rates and rotational speeds; these show some discrepancies with the experiments, although the calculations exhibit the more consistent trend. Further calculations indicate that differences in thermal boundary conditions have a greater influence on Nusselt number than differences in disc geometry. The influence of the bolt cover on the heat transfer has also been modelled, although comparative measurements are not available.


Sign in / Sign up

Export Citation Format

Share Document