Magnetohydrodynamic co-rotating flow in a vertical cylindrical container

ABSTRACT

The effect of an external axial magnetic field on the liquid metal flow produced by co-rotation of the top and bottom disks in a vertical cylindrical container with a vertical temperature gradient is numerically analyzed. The governing Navier–Stokes, energy, and potential equations along with appropriate boundary conditions are solved using the finite-volume method. Comparisons with the previous results were performed and found to be in excellent agreement. It was observed that the Reynolds number is increased, and the axisymmetric basic state loses stability for circular patterns of axisymmetric vortices and spiral waves. In the mixed convection case the axisymmetric mode disappears, giving an asymmetric mode m = 1. It was also found that the primary thresholds, Re_cr corresponding to modes m = 1 and 2, increase with an increase in Hartmann number (Ha). We can therefore conclude that when the magnitude of the magnetic field exceeds a certain value, the instability becomes a steady bifurcation. Finally, stability diagrams were established according to the numerical results of this investigation. These diagrams show the evolution of primary thresholds as a function of Hartmann number for various values of Richardson number.     

Stability of Swirling Flows with Heat Transfer in a Cylindrical Enclosure with CO/Counter-Rotating End Disks Under an Axial Magnetic Field

In this article, a numerical study of swirling flows with heat transfer generated by two rotating end disks (co- and counter-rotating) inside a cylindrical enclosure having an aspect ratio equal to 2, filled with a liquid metal, and submitted to a vertical temperature gradient and an axial magnetic field is studied. The governing Navier-Stokes, energy, and potential equations along with appropriate boundary conditions are solved by using the finite-volume method. The flow and temperature fields are presented by stream function and isotherms, respectively. This flow is very unstable and reveals a great richness of structures. In an oscillatory regime, results are presented for various values of the Hartmann number, Ha = 5, 10, 20 and 30, and Richardson numbers, Ri = 0, 0.5, 1, 2 and 4, in order to see their effects on the value of the critical Reynolds number, Re_cr. Stability diagrams are established according to the numerical results of this investigation. These diagrams show the dependence of Re_cr with the increase of Ha for various values of Ri. The flow between co-rotating end disks is very different from the flow between counter-rotating end disks. Finally, this study confirms the possibility of stabilization of a liquid metal flow by application of an axial magnetic field.     

Wall–bounded flow and heat transfer around a circular cylinder at low Reynolds and Hartmann numbers

A two–dimensional numerical simulation is performed following a finite volume approach to analyze the forced convection heat transfer for the hydromagnetic flow around a circular cylinder at low Reynolds numbers. The cylinder is placed within a rectangular channel subjected to externally applied magnetic fields and acted upon by the magnetohydrodynamic (MHD) flow of a viscous incompressible and electrically conductive fluid. The magnetic field is applied either along the streamwise or transverse directions. The simulation is carried out for the range of Reynolds number 10 ≤ Re ≤ 80 with Hartmann number 0 ≤ Ha ≤ 10 and for different Prandtl numbers, Pr = 0.02 (liquid metal), 0.71 (air), and 7 (water) for a blockage ratio β = 0.25. The flow is steady for the above range of conditions. Apart from the channel wall, the magnetic field provides additional stability to the flow as a result of which the recirculation region behind the obstacle reduces with increasing magnetic field strength for a particular Reynolds number. The rate of heat transfer is found almost invariant at low Re whereas it increases slightly for higher Re with the applied magnetic field. The heat transfer increases as usual with the Reynolds number for all Hartmann numbers. © 2012 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.21025     

Magnetohydrodynamic natural convection in a vertical cylindrical cavity with sinusoidal upper wall temperature

A series of numerical simulations were performed in order to study liquid metal MHD natural convection in a vertical cylindrical container with a sinusoidal temperature distribution at the upper wall and the other surfaces being adiabatic. Starting from the basic hydrodynamic case, the effect of vertical (axial) and horizontal magnetic fields is assessed. Depending on the magnitude of the Rayleigh and Hartmann numbers, both turbulent and laminar (azimuthally symmetric or not) flows are observed. The results show that the increase of Rayleigh number promotes heat transfer by convection while the increase of Hartmann number favors heat conduction. The vertical magnetic field reduces the Nusselt number more than the horizontal. The circulation patterns for the most convective cases are confined close to the top corner of the container with the simultaneous formation of a secondary flow pattern at the bottom corner, while for the more conductive cases only one circulation pattern exists covering the entire domain.     

Wall‐Confined Flow and Heat Transfer Around a Square Cylinder at Low Reynolds and Hartmann Numbers

This article discusses the results obtained through a two‐dimensional numerical simulation following a finite volume approach on the forced convection heat transfer for the hydromagnetic flow around a square cylinder at low Reynolds and Hartmann numbers. The magnetohydrodynamic (MHD) flow of a viscous incompressible and electrically conducting fluid is assumed to take place in a rectangular channel subjected to externally imposed magnetic fields and the cylinder is fixed within the channel. The magnetic fields may be applied either along the streamwise or transverse directions. Simulations are performed for the range of kinetic Reynolds number 10 ≤ Re ≤ 60 with Hartmann number 0 ≤ Ha ≤ 15 and for different thermal Prandtl numbers, Pr = 0.02 (liquid metal), 0.71 (air), and 7 (water) for a blockage ratio β = 0.25. A steady flow can be expected for the above range of conditions. Besides the channel wall, the magnetic field imparts additional stability to the flow as a consequence of which the recirculation region behind the obstacle reduces with increasing magnetic field strength for a particular Re. The critical Hartmann numbers for the complete suppression of flow separation in the case of a transversely applied magnetic field are computed. The rate of heat transfer is found almost invariant at low Re whereas it increases moderately for higher Re with the applied magnetic field. The heat transfer increases in general with the Reynolds number for all Hartmann numbers. Finally, the influence of obstacle shape on the thermohydrodynamic quantities is noted. © 2013 Wiley Periodicals, Inc. Heat Trans Asian Res, 43(5): 459–475, 2014; Published online 3 October 2013 in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.21091     

Simulation of a cavity insulated by a vertical single band cold air curtain

The cold air curtain used for cavity insulation in the present study is a fundamental idealization of a refrigerated air curtain in display cases. Flow and heat transfer characteristics of vertical single band cold air curtains are studied numerically to investigate the effects of mixed convection and height/width ratios. The results show that for a given Grashof number, the Richardson number must be less than a critical value to assure thermal insulation. From studies of the effect of height/width ratios, it is found that small height/width ratios yield an increased critical Richardson number (thus a reduced critical Reynolds number) and a decreased volumetric infiltration rate. For the cavity whose height/depth ratio equals 2 in the present study, the critical Richardson numbers are 0.62, 0.28 and 0.20 for the height/width ratios of 10, 15 and 20, respectively.     

Mixed convection flow and heat transfer across a square cylinder under the influence of aiding buoyancy at low Reynolds numbers

In this paper, mixed convection flow and heat transfer around a long cylinder of square cross-section under the influence of aiding buoyancy are investigated in the vertical unconfined configuration (Reynolds number, Re = 1–40 and Richardson number, Ri = 0–1). The semi-explicit finite volume method implemented on the collocated grid arrangement is used to solve the governing equations along with the appropriate boundary conditions. The onset of flow separation occurs between Re = 1–2, between Re = 2–3 and between Re = 3–4 for Ri = 0, 0.5 and 1, respectively. The flow is found to be steady for the range of conditions studied here. The friction, pressure and total drag coefficients are found to increase with Richardson number, i.e., as the influence of aiding buoyancy increases drag coefficients increase at the constant value of the Reynolds number. The temperature field around the obstacle is presented by isotherm contours at the Prandtl number of 0.7 (air). The local and average Nusselt numbers are calculated to give a detailed study of heat transfer over each surface of the square cylinder and an overall heat transfer rate and it is found that heat transfer increases with increase in Reynolds number and/or Richardson number. The simple expressions for the wake length and average cylinder Nusselt number are obtained for the range of conditions covered in this work.     

Acknowledgment

An experimental and numerical study has been carried out in order to investigate mixed and natural convection heat transfer in a two-dimensional enclosure. A discrete isothermal heat source is located at one of the vertical walls. Also, two ventilation ports are at the bottom and on top of the opposite wall. A forced flow condition was imposed by providing an inlet of air at the bottom port. A Mach–Zehnder interferometer was used to visualize the temperature field within the enclosure and to determine the local and average heat transfer characteristics of the heat source. Five heater positions on the vertical wall and different Rayleigh numbers (4.5 × 10⁵ to 1.15 × 10⁶) and Reynolds numbers (120 to 1600) were considered in the experiments. A finite volume code has been developed based on the SIMPLE algorithm and hybrid discretization scheme for the numerical study. It is observed that the interaction of natural convection with the forced flow leads to various flow fields depending on the Richardson number, Reynolds number and the heater position. Also, results show different trends for variation of the average Nusselt number with the heater position at low and high Reynolds numbers. An optimum position for the heat source, at which the maximum heat transfer is achieved, exists for high Reynolds numbers and has been found to be at the middle of the vertical wall.     

Three-Dimensional Natural Convection in the Horizontal Bridgman Configuration Under Various Wall Electrical Conductivity and Magnetic Field

A numerical investigation of the three-dimensional natural convection of a liquid metal contained in the horizontal Bridgman configuration, having an aspect ratio equal to 5 and submitted to an external magnetic field in either the longitudinal or vertical direction, is presented. The numerical approach is based on the finite-volume approximation. A computer program based on the SIMPLER algorithm is developed. The effect of a magnetic field provides a notable change on the flow and thermal structures. The strongest stabilization of the convection flow is found when the magnetic field is oriented vertically. Also, wall electrical conductivity has an effect on the average Nusselt number. A good agreement between our numerical simulations and experimental data found in the literature is obtained.     

Two-dimensional laminar fluid flow and heat transfer in a channel with a built-in heated square cylinder

A numerical investigation was conducted to analyze the unsteady flow field and heat transfer characteristics in a horizontal channel with a built-in heated square cylinder. Hydrodynamic behavior and heat transfer results are obtained by the solution of the complete Navier–Stokes and energy equations using a control volume finite element method (CVFEM) adapted to the staggered grid. The Computation was made for two channel blockage ratios (β=1/4 and 1/8), different Reynolds and Richardson numbers ranging from 62 to 200 and from 0 to 0.1 respectively at Pr=0.71. The flow is found to be unstable when the Richardson number crosses the critical value of 0.13. The results are presented to show the effects of the blockage ratio, the Reynolds and the Richardson numbers on the flow pattern and the heat transfer from the square cylinder. Heat transfer correlation are obtained through forced and mixed convection.     

Lattice Boltzmann simulation of MHD mixed convection in a two‐sided lid‐driven square cavity

In this paper the effects of a magnetic field on mixed convection flow in a two‐sided lid‐driven cavity have been analyzed by the lattice Boltzmann method (LBM). The Hartmann number varied from Ha = 0 to 100. The study has been conducted for different Richardson numbers (Ri) from 0.01 to 100 while the direction of the magnetic field was investigated in the x‐direction. Consequences demonstrate that the heat transfer augments with an increment of the Richardson number for different Hartmann numbers for two cases. The heat transfer declines with the growth of the magnetic field for various Richardson numbers for two cases. The difference between the values of heat transfer for the two cases at variant parameters is negligible but the trend of fluid flow for the two cases is multifarious. © 2011 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.20402     

Convection and instability of thermocapillary flow in a liquid bridge subject to a non-uniform rotating magnetic field

A three-dimensional liquid bridge is considered in this study to numerically investigate the effects of an external non-uniform rotating magnetic field (RMF) on the thermocapillary flow in semiconductor melt under microgravity. Simulations are carried out to examine the convection and instability features of the thermocapillary flow over a range of Marangoni numbers (Ma = 15–50) under a non-uniform RMF. The present results show that applying an external non-uniform RMF enhances the maximum tangential velocity and depresses the maximum axial velocity. As a consequence, an approximately axisymmetric flow is maintained in the melt under the effect of the non-uniform RMF, which is beneficial for growing high quality crystal. Further investigation of the thermocapillary flow subject to different non-uniform RMFs (corresponding to Taylor numbers Ta = 3.8 × 10²–1.86 × 10⁴ and Rotating Reynolds number Re_ω = 2.2 × 10⁴) reveals that the thermocapillary convection may undergo a transition from the approximately axisymmetric steady flow to a periodically oscillatory flow for Ma above a critical value. The critical Ma generally increases with the intensity of the non-uniform RMF.     

Analysis of the flow and heat transfer characteristics for assisting incompressible laminar flow past an open cavity

A numerical study has been carried out to analyze the effects of mixed convective assisting flow past three-dimensional open cavity over a wide range of Reynolds (100–1000) and Richardson (0.001–10) numbers. The vertical walls in the inflow and outflow sides are isothermal while all other walls are adiabatic. The cavity is assumed to be cubic in geometry and the flow is laminar. A direct numerical simulation is undertaken to investigate the flow structure, the heat transfer characteristics and the complex interaction between the induced stream flow at ambient temperature and the buoyancy-induced flow from the heated wall. It is found that the flow becomes stable at moderate Grashof number and exhibits a three-dimensional structure, while for high Richardson number the mixed convection effects come into play and push the recirculating zone further upstream and the flow may becomes unstable.     

Experimental investigation of forced and mixed convection heat transfer in a foam-filled horizontal rectangular channel

An experimental study was performed to investigate the heat transfer characteristics of the mixed convection flow through a horizontal rectangular channel where open-cell metal foams of different pore densities (10, 20 and 30 PPI) were situated. A uniform heat flux was applied at all of the bounding walls of the channel. For each of three values of the uniform heat flux, temperatures were measured on the entire surfaces of the walls. Results for the average and local Nusselt numbers are presented as functions of the Reynolds and Richardson numbers. The Reynolds number based on the channel height of the rectangular channel was varied from 600 to 33000, while the Richardson number ranged from 0.02 to 103, extending over forced, mixed and natural convection. Second important parameter that influences the heat transfer is the aspect ratio of the foams. Three different aspect ratios (AR) as 0.25, 0.5 and 1 are tested. Based on the experimental data, new empirical correlations have been constructed to link the Nusselt number. The results of all cases were compared to that of the empty channel and the literature. We found that our results were in agreement with those that are mentioned in the literature.     

Three dimensional lattice Boltzmann simulation for mixed convection of nanofluids in the presence of magnetic field

In the present study, a three dimensional thermal lattice Boltzmann model was developed to investigate the flow dynamics and mixed convection heat transfer of Al₂O₃/water nanofluid in a cubic cavity in the presence of magnetic field. The model was first validated with previous numerical and experimental results. Satisfactory agreement was obtained. Then the effects of Rayleigh number, nanoparticle volume fraction, Hartmann number and Richardson number on nanofluid flow dynamics and heat transfer were examined. Numerical results indicate that adding nanoparticles to pure water leads to heat transfer enhancement for low Rayleigh numbers. However, this enhancement might be weakened and even reversed for high Rayleigh numbers. In addition, the results show the external applied magnetic field has an effect of suppressing the convective heat transfer in the cavity. Moreover, the results demonstrate that the Richardson number in mixed convection has significant influences on both streamlines and temperature field.     

Mixed convection heat transfer in a differentially heated square enclosure with a conductive rotating circular cylinder at different vertical locations

In this investigation, a numerical simulation using a finite volume scheme is carried out for a laminar steady mixed convection problem in a two-dimensional square enclosure of width and height (L), with a rotating circular cylinder of radius (R = 0.2 L) enclosed inside it. The solution is performed to analyze mixed convection in this enclosure where the left side wall is subjected to an isothermal temperature higher than the opposite right side wall. The upper and lower enclosure walls are considered adiabatic. The enclosure under study is filled with air with Prandtl number is taken as 0.71. Fluid flow and thermal fields and the average Nusselt number are presented for the Richardson numbers ranging as 0, 1, 5 and 10, while Reynolds number ranging as 50, 100, 200 and 300. The effects of various locations and solid-fluid thermal conductivity ratios on the heat transport process are studied in the present work. The results of the present investigation explain that increase in the Richardson and Reynolds numbers has a significant role on the flow and temperature fields and the rotating cylinder locations have an important effect in enhancing convection heat transfer in the square enclosure. The results explain also, that the average Nusselt number value increases as the Reynolds and Richardson numbers increase and the convection phenomenon is strongly affected by these parameters. The results showed a good agreement with further published works.     

Numerical study of fluid flow and heat transfer in microchannel cooling passages

Numerical investigation was conducted for fluid flow and heat transfer in microchannel cooling passages. Effects of viscosity and thermal conductivity variations on characteristics of fluid flow and heat transfer were taken into account in theoretical modeling. Two-dimensional simulation was performed for low Reynolds number flow of liquid water in a 100 μm single channel subjected to localized heat flux boundary conditions. The velocity field was highly coupled with temperature distribution and distorted through the variations of viscosity and thermal conductivity. The induced cross-flow velocity had a marked contribution to the convection. The heat transfer enhancement due to viscosity-variation was pronounced, though the axial conduction introduced by thermal-conductivity-variation was insignificant unless for the cases with very low Reynolds numbers.     

Numerical simulation of double-diffusive mixed convective flow in rectangular enclosure with insulated moving lid

The present study is concerned with the mixed convection in a rectangular lid-driven cavity under the combined buoyancy effects of thermal and mass diffusion. Double-diffusive convective flow in a rectangular enclosure with moving upper surface is studied numerically. Both upper and lower surfaces are being insulated and impermeable. Constant different temperatures and concentration are imposed along the vertical walls of the enclosure, steady state laminar regime is considered. The transport equations for continuity, momentum, energy and spices transfer are solved. The numerical results are reported for the effect of Richardson number, Lewis number, and buoyancy ratio on the iso-contours of stream line, temperature, and concentration. In addition, the predicted results for both local and average Nusselt and Sherwood numbers are presented and discussed for various parametric conditions. This study was done for 0.1 ≤ Le ≤ 50 and Prandtl number Pr = 0.7. Through out the study the Grashof number and aspect ratio are kept constant at 10⁴ and 2 respectively and ?10 ≤ N ≤ 10, while Richardson number has been varied from 0.01 to 10 to simulate forced convection dominated flow, mixed convection and natural convection dominated flow.     

Numerical prediction of flow and heat transfer of power-law fluids in a plane channel with a built-in heated square cylinder

Structure of unsteady laminar flow and heat transfer of power-law fluids in two-dimensional horizontal plane channel with a built-in heated square cylinder is studied numerically. The governing equations are solved using a control volume finite element method (CVFEM) adapted to the staggered grid. Computations are performed over a range of Reynolds and Richardson numbers from Re = 20 to 200 and from Ri = 0 to 8, respectively at fixed Prandtl number Pr = 50 and blockage ratio value β′ = 1/8. Three different values of the power-law index (n = 0.5, 1 and 1.4) are considered in this study to show its effect on the value of the critical Reynolds number defining the transition between two different flow regimes (symmetrical and periodic flows), the variations of Strouhal number, drag and lift coefficients and the heat transfer from the square cylinder as function of Reynolds number. Heat transfer correlations are obtained through forced convection. A discussion about the buoyancy effect on the flow pattern and the heat transfer for different power-law index is also presented.     

Direct and inverse mixed convections in an enclosure with ventilation ports

The numerical study presented in this work describes the direct and inverse mixed convection problems in a slot-ventilated enclosure subjected to an unknown heat flux on one side. Particularly, the interaction of internal natural convection with the cold ventilated flow leads to various flow fields depending on the Richardson number, Reynolds number, and the functional form of the imposed boundary heat flux. Fluid and heat transport structures across the enclosure are visualized by the streamlines and heatlines, respectively. Subsequently, an iterative conjugate gradient method is applied such that the gradient of the cost function is introduced when the appropriate sensitivity and adjoint problems are defined for a domain of arbitrary geometries. In this approach, no a priori information is needed about the unknown boundary heat fluxes to be determined. The accuracy of the heat flux profile solutions is shown to depend strongly on the values of Reynolds number and flux functional forms. Effects of measurement errors on the accuracy of estimation are also investigated. The present work is significant for the flow control simultaneously involving the natural convection and forced convection.