Integrating thermal-concentration smoothed profile with lattice Boltzmann methods for simulating sedimentation of nonisothermal circular particles

A thermal-concentration smoothed profile-lattice Boltzmann method is proposed to study the effect of the concentration field on the dynamic behavior of nonisothermal cylindrical particles during the sedimentation process. The velocity, temperature, and concentration equations are solved using the lattice Boltzmann method. Moreover, the smoothed profile method is employed to enforce the nonslip boundary condition as well as constant temperature and constant concentration boundary conditions at the particles surfaces. Moreover, the Boussinesq approximation is used to couple the velocities, temperatures, and concentrations fields. The proposed combined method is validated by comparing the present numerical results with those found in the literature, showing good consistency. Then, the effect of the concentration buoyancy on the behavior of nonisothermal particles is discussed. In addition, the effect of Prandtl, Schmidt, and thermal Grashof numbers on the settling process is investigated. The results show that, by adding the effect of concentration, the maximum settling velocity of hot particles is reduced more relative to the cold ones; accordingly, the cold particles are settled faster than the hot ones. Finally, the sedimentation of two particles in a container at high thermal Grashof is investigated. It is shown that, at high thermal Grashof, there is an intense competition between the buoyancy force and gravity for the hot particles. The buoyancy flow generated leads to the reversal of the drafting-kissing-tumbling motion of the hot particles, making the particles move upward.     

Non-Darcy effects on convective boundary layer flow past a semi-infinite vertical plate in saturated porous media

The composite effects of viscosity, porosity, buoyancy parameter, thermal conductivity ratio and non-Darcy effects of Brinkman friction and Forscheimmer quadratic drag on the mixed convection boundary layer flow past a semi-infinite plate in a fully-saturated porous regime are theoretically and numerically investigated using Keller’s implicit finite-difference technique and a double-shooting Runge-Kutta method. The Brinkman Forcheimer-extended Darcy model is implemented in the hydrodynamic boundary layer equation. The effects of the various non-dimensional thermofluid parameters, viz Grashof number, Darcy number, and Forchheimer number, and also porosity, thermal conductivity and viscosity parameters on the velocity and temperature fields are discussed. Computations for both numerical schemes are made where possible and found to be in excellent agreement.     

Mixed convection adjacent to inclined flat surfaces embedded in a porous medium

An analysis is performed to study the flow and heat transfer characteristics of laminar mixed convection boundary layer flows from inclined (including horizontal and vertical) surfaces embedded in a saturated porous medium with constant aiding external flows and uniform surface temperature. Both the streamwise and normal components of the buoyancy forces are retained in the momentum equations. Nondimensionalization of the boundary layer equations results in the following three governing parameter: (1)Gr/Re, the ratio of the Grashof number to the Reynolds number; (2)Pe _x =Re _x Pr, the Peclet number; (3) φ, the angle of inclination from the horizontal. The resulting nonsimilar equations are solved by an efficient implicit finite-difference scheme. Numerical results are presented for flows with different values ofGr/Re in the range of 0 to 50, over a wide range of the Peclet numbersPe _x, and various values of φ ranging from 0 to 90 degrees. It is found that the local surface heat transfer rate increases with increasing the local Peclet number. In addition, as the plate is tilted from the horizontal to the vertical orientation, the local Nusselt number increases for a given Peclet number and the effect of the buoyancy force on the surface heat transfer rate increases.     

The effect of buoyancy forces on laminar convection with viscous dissipation in a vertical channel

The interaction of the viscous dissipation effect with the presence of buoyancy forces is investigated for laminar-flow heat transfer in a parallel-plate vertical channel. One of the channel walls is considered as isothermal with a prescribed temperature, while the other wall is considered as insulated. The velocity field is assumed to be parallel. The velocity field, the temperature field and the Nusselt number are obtained by a perturbation series method which employs the ratio between the Grashof number and the Reynolds number as the perturbation parameter. The radius of convergence of the perturbation series is estimated. Received on 10 December 1997     

Higher approximations to the free convection flow from a heated vertical flat plate

This paper is concerned with the problem of obtaining higher approximations for the free convection from a heated vertical flat plate to that represented by the well known solution of Schmidt and Beckmann. For large Grashof number, the perturbation problem is a singular one and the method of matched asymptotic expansions is used to construct inner and outer expansions for the velocity and temperature distributions. The small perturbation parameterε is the inverse of the fourth root of the Grashof number and the expansions are shown to involve only integral powers ofε. The first three terms in the expansion are calculated and numerical results are presented for the velocity, temperature, skin friction and heat transfer. The agreement with experiment is found to be excellent, and the theory fully explains the discrepancies which exist between boundary layer theory and experiment.     

Equilibrium characteristics of the secondary convection in a vertical fluid layer between two flat plates 1

The nonlinear stability of the natural convection in a vertical fluid layer between two flat plates with different temperatures is investigated by a direct method to find the equilibrium states of the secondary convection. We confine ourselves to two-dimensional flows and assume that the aspect ratio of the fluid layer is very large. Since the Prantl number is assumed to be very small, the buoyancy effect caused by temperature disturbances is negligible. As a result we obtained a neutral surface of the energy of the fundamental mode of the secondary convection. It is concluded that there is no finite amplitude instability below the critical Grashof number derived from linear stability theory, and that both the unstable equilibrium solution (threshold amplitude solution) and the stable equilibrium solution (finite amplitude solution) are found outside the neutral curve of the linear stability. Our results are almost consistent with those of Nagata and Busse (1983), but are more accurate and more thorough.     

Stability properties of the vertical boundary layers in differentially heated cavities

In the present study, the two-dimensional (2-D) stability properties of the vertical boundary layers in a cavity that is differentially heated over two opposing vertical walls is considered. The study is performed by introducing artificial, controlled perturbations at the base of the vertical boundary layer along the hot cavity wall and by following the evolution of these disturbances. For small initial perturbations, the evolution is governed by linear effects. This method accurately predicts the frequency of the bifurcation, which occurs for (much) larger Rayleigh numbers. Convective instability sets in for Rayleigh numbers much smaller than those at which the absolute instability (i.e., the bifurcation) occurs, and these Rayleigh numbers are in reasonable agreement with those for the boundary layer along a plate. The absolute instability does not result from the first wave which becomes unstable. For small Prandtl numbers (≤ 2), the unstable waves which lead to the absolute instability are shear-driven, and a single frequency is introduced in the flow after the bifurcation. For larger Prandtl numbers, the unstable waves are buoyancy driven and no single-frequency unsteady flow is observed after the bifurcation.     

Mixed convection flow in vertical channel with boundary conditions of third kind in presence of heat source/sink

The effects of viscous dissipation and heat source/sink on fully developed mixed convection for the laminar flow in a parallel-plate vertical channel are investigated.The plate exchanges heat with an external fluid.Both conditions of equal and different reference temperatures of the external fluid are considered.First,the simple cases of the negligible Brinkman number or the negligible Grashof number are solved analytically.Then,the combined effects of buoyancy forces and viscous dissipation in the presence of heat source/sink are analyzed by a perturbation series method valid for small values of the perturbation parameter.To relax the conditions on the perturbation parameter,the velocity and temperature fields are solved by using the Runge-Kutta fourth-order method with the shooting technique.The velocity,temperature,skin friction,and Nusselt numbers at the plates are discussed numerically and presented through graphs.     

Heat and mass transfer on a MHD third grade fluid with partial slip flow past an infinite vertical insulated porous plate in a porous medium

The influence of third grade, partial slip and other thermophysical parameters on the steady flow, heat and mass transfer of viscoelastic third grade fluid past an infinite vertical insulated plate subject to suction across the boundary layer has been investigated. The space occupying the fluid is porous. The momentum equation is characterized by a highly nonlinear boundary value problem in which the order of the differential equation exceeds the number of available boundary conditions. An efficient numerical scheme of midpoint technique with Richardson’s extrapolation is employed to solve the governing system of coupled nonlinear equations of momentum, energy and concentration. Numerical calculations were carried out for different values of various interesting non-dimensional quantities in the slip flow regime with heat and mass transfer and were shown with the aid of figures. The values of the wall shear stress, the local rate of heat and mass transfers were obtained and tabulated. The analysis shows that as the fluid becomes more shear thickening, the momentum boundary layer decreases but the thermal boundary layer increases; the magnetic field strength is found to decrease with an increasing temperature distribution when the porous plate is insulated. The consequences of increasing the permeability parameter and Schmidt number decrease both the momentum and concentration boundary layer thicknesses respectively whereas an increase in the thermal Grashof number gives rise to the thermal boundary layer thickness.     

Effects of thermal buoyancy and variable thermal conductivity on the MHD flow and heat transfer in a power-law fluid past a vertical stretching sheet in the presence of a non-uniform heat source

The paper considers the flow of a power-law fluid past a vertical stretching sheet. Effects of variable thermal conductivity and non-uniform heat source/sink on the heat transfer are addressed. The thermal conductivity is assumed to vary linearly with temperature. Similarity transformation is used to convert the governing partial differential equations into a set of coupled, non-linear ordinary differential equations. Two different types of boundary heating are considered, namely Prescribed power-law Surface Temperature (PST) and Prescribed power-law Heat Flux (PHF). Shooting method is used to obtain the numerical solution for the resulting boundary value problems. The effects of Chandrasekhar number, Grashof number, Prandtl number, non-uniform heat source/sink parameters, wall temperature parameter and variable thermal conductivity parameter on the dynamics are shown graphically in several plots. The skin friction and heat transfer coefficients are tabulated for a range of values of the parameters. Present study reveals that in a gravity affected flow buoyancy effect has a significant say in the control of flow and heat transfer.     

Hydromagnetic free convection flow with induced magnetic field effects

An exact solution is presented for the hydromagnetic natural convection boundary layer flow past an infinite vertical flat plate under the influence of a transverse magnetic field with magnetic induction effects included. The transformed ordinary differential equations are solved exactly, under physically appropriate boundary conditions. Closed-form expressions are obtained for the non-dimensional velocity (u), non-dimensional induced magnetic field component (B _x ) and wall frictional shearing stress i.e. skin friction function (τ _x ) as functions of dimensionless transverse coordinate (η), Grashof free convection number (G _r ) and the Hartmann number (M). The bulk temperature in the boundary layer (Θ) is also evaluated and shown to be purely a function of M. The Rayleigh flow distribution (R) is derived and found to be a function of both Hartmann number (M) and the buoyant diffusivity parameter (ϑ ^*). The influence of Grashof number on velocity, induced magnetic field and wall shear stress profiles is computed. The response of Rayleigh flow distribution to Grashof numbers ranging from 2 to 200 is also discussed as is the influence of Hartmann number on the bulk temperature. Rayleigh flow is demonstrated to become stable with respect to the width of the boundary layer region and intensifies with greater magnetic field i.e. larger Hartman number M, for constant buoyant diffusivity parameter ϑ ^*. The induced magnetic field (B _x ), is elevated in the vicinity of the plate surface with a rise in free convection (buoyancy) parameter G _r , but is reduced over the central zone of the boundary layer regime. Applications of the study include laminar magneto-aerodynamics, materials processing and MHD propulsion thermo-fluid dynamics.     

Combined forced and free convective flow of water at 4°C past a semi-infinite vertical plate

Effects of buoyancy forces on forced and free convective flow of water at 4°C past a semi-infinite vertical plate at constant temperature are studied. Flow is assumed to be vertically upwards. Similarity solutions are derived and the resulting equations are solved numerically on a computer. Velocity and temperature profiles are shown graphically and numerical values of the skin friction and the rate of heat transfer are entered in tables. It is observed that the skin friction and the Nusselt number increase with increasing Gr/Re², where Gr is the Grashof number and Re is the Reynolds number     

Vortex shedding around a heated square cylinder under the influence of buoyancy

The influence of buoyancy on vortex shedding and heat transfer from a cylinder of square cross-section exposed to a horizontal stream has been studied.Unsteady Navier-Stokes and energy equations are solved numerically using a control volume approach. Flow field has been analysed for a wide range of Reynolds number (which is based on the cross-sectional height of the cylinder) and Grashof number with Richardson number between 0 to 1. Our results show that the centerline symmetry of the wake is lost and the cylinder experiences a downwards lift when the buoyancy effect is considered. Vortex shedding suppression doesnt occur in the present case in which the cylinder is exposed to a horizontal cross-flow. Heat transfer from the cylinder increases due to increase in Reynolds number and Grashof number.     

Non-darcy natural convection on a vertical cylinder in a saturated porous medium

The steady free convection boundary layer flow of non-Darcy fluid along an isothermal vertical cylinder embedded in a saturated porous medium using the Ergun model has been studied. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme developed by Keller. It is found that the heat transfer is strongly affected by the modified Grashof number which characterizes the non-Darcy fluid, and the curvature parameter. Also the heat transfer is found to be more than that of the flat plate.     

Buoyancy effects on mixed convection heat and mass transfer in a duct with sudden expansions

The present paper deals with the study of heat and mass transfer by mixed convection in an inclined duct preceded with a double step expansion. The control volume based finite element method was used to solve the set of non-dimensional equations for the vorticity, stream function, energy and species conservation. Numerical simulations are carried out for different combinations of the Lewis number, thermal and mass diffusion Grashof numbers for different inclinations. Streamlines, temperature and concentration distributions are presented and discussed. The results show the effect of the secondary flow induced by buoyancy forces and the presence of the double step expansion on the heat and mass transfer mechanism. It is found that the recirculation vortices induced by the expansion can be present along the channel and the flow structure can be wavy. For the vertical orientation, asymmetric fields are observed for the different simulated cases.     

Mixed Convection in a Darcy–Brinkman Porous Medium with a Constant Convective Thermal Boundary Condition

The characteristics of the boundary layer flow past a plane surface adjacent to a saturated Darcy–Brinkman porous medium are investigated in this paper. The flow is driven by an external free stream moving with constant velocity. The surface is heated with a convective boundary condition with constant heat transfer coefficient. The problem is non-similar and is investigated numerically by a finite difference method. The problem is governed by four non-dimensional parameters, that is, the convective Darcy number, the convective Grashof number, the Prandtl number, and the axial distance along the plate. The influence of these parameters on the results is investigated, and the results are presented in tables and figures. The Darcy term and the Grashof term in the momentum equation contradict each other and this contradiction makes the problem complicated. However, the wall shear stress and the wall temperature increase continuously along the plate and the wall temperature always tends to 1.     

Combined forced and free flow in a vertical rectangular duct with prescribed wall heat flux

In this paper, combined forced and free convection is studied in a vertical rectangular duct with a prescribed uniform wall heat flux (H2 boundary condition). A different heat flux value for each plane wall is considered; the condition of a uniform wall heat flux throughout the duct results as a special case. The local momentum and energy balance equations are written in a dimensionless form and solved numerically, by means of a Galerkin finite element method. The numerical solution gives the dimensionless velocity and temperature distributions, together with the values of the Fanning friction factor, of the Nusselt number, of the momentum flux correction factor and of the kinetic energy correction factor. These dimensionless parameters are reported as functions of the aspect ratio and of the ratio between the Grashof number, Gr, and the Reynolds number, Re. The threshold values of Gr/Re for the onset of flow reversal are evaluated.     

Buoyancy effects on the mass transfer in absorption with a nonabsorbable gas

This paper describes the vapor side buoyancy effects on the mass transfer in absorption in the presence of a nonabsorbable gas. Experimental results on a diffusion-absorption refrigerator (DAR) indicate that the vapor side buoyancy effects on the mass transfer are significant when the density of the nonabsorbable gas is significantly less than that of the absorbable gas. A rectangular enclosure absorption problem is first solved to demonstrate the buoyancy effects without the presence of a forced flow. Then, mixed convection heat transfer in a circular pipe is simulated in such a way as to be analogous to the mixed convection mass-transfer problem in the DAR absorber. Finally, the vapor side mixed convection absorption between parallel plates is simulated including the effects of the absorbed mass on the mass balance. The Sherwood number dependence on the mass transfer Grashof number and Reynolds number as well as the effects of the suction boundary conditions are discussed. Each of these simulations had individual limitations, but, taken together, they illuminate the major aspects of the absorption physics.     

Three dimensional unsteady convection and mass transfer flow through porous medium

The problem of three dimensional unsteady convection flow through a porous medium, with effect of mass transfer bounded by an infinite vertical porous plate is discussed, when the suction at the plate is transverse sinusoidal and the plate temperature oscillates in time about a constant mean. Assuming the free stream velocity to be uniform, approximate solutions are obtained for the flow field, the temperature field, the skin-friction and the rate of heat transfer. The dependence of solution on Pr (Prandtl number), Gr (Grashof number based on temperature), Gc (modified Grashof number based on concentration difference), Sc (Schimdt number), the frequency and the permeability parameter is also investigated.