Lattice Boltzmann Simulation of Airflow and Heat Transfer in a Model Ward of a Hospital

Author(s):  
Md. Farhad Hasan ◽  
Taasnim Ahmed Himika ◽  
Md. Mamun Molla

In this research, a very popular alternative computational technique, the lattice Boltzmann method (LBM), has been used to simulate the indoor airflow and heat transfer in a model hospital ward. Different Reynolds numbers have been used to study the airflow pattern. Boundary conditions for velocity and temperature have also been discussed in detail. Several tests have been conducted for code validation. LBM is demonstrated through simulation in forced convection inside hospital ward with six beds for two different situations: ward without partition and ward with partition. Changes in average rate of heat transfer in terms of average Nusselt numbers have also been recorded for those situations. Average Nusselt numbers were found to differ for different cases. In terms of airflow, it has been found that, for various Reynolds numbers, airflow changes its pattern and leads to few recirculations for relatively higher Reynolds number but remains steady for low Reynolds number. It was observed that partition narrowed the channel for airflow and once the air overcame this barrier, it gets free space and recirculation appears more. For higher Reynolds number, the average rate of heat transfer increases and patients near the recirculation zone release maximum heat and will feel more comfortable.

2016 ◽  
Vol 2016 ◽  
pp. 1-15 ◽  
Author(s):  
Taasnim Ahmed Himika ◽  
Md. Farhad Hasan ◽  
Md. Mamun Molla

In the present investigation the airflow and heat transfer for mixed convection have been simulated for a model general ward of hospital with six beds and partitions using the Lattice Boltzmann Method (LBM). Three different Reynolds numbers 100, 250, and 350 have been considered. Bounce-back condition has been applied at the wall. Results have been represented in three different case studies and the changes have been discussed in terms of streamlines and isotherms. Code validation has also been included before going through the simulation process and it shows good agreement with previously published papers when the comparison is made on average Nusselt number. Results show that the pattern of indoor airflow is varied in each and every case study due to the effect of mixed convection flow and placement of partition. In addition, the changes in average rate of heat transfer indicate that patients closer to inlet get the most air and feel better and if any patient does not need much air, he or she should be kept near the outlet to avoid temperature related complications.


1999 ◽  
Vol 122 (2) ◽  
pp. 375-385 ◽  
Author(s):  
C. R. Hedlund ◽  
P. M. Ligrani

Local flow behavior and heat transfer results are presented from two swirl chambers, which model passages used to cool the leading edges of turbine blades in gas turbine engines. Flow results are obtained in an isothermal swirl chamber. Surface Nusselt number distributions are measured in a second swirl chamber (with a constant wall heat flux boundary condition) using infrared thermography in conjunction with thermocouples, energy balances, and in situ calibration procedures. In both cases, Reynolds numbers Re based on inlet duct characteristics range from 6000 to about 20,000. Bulk helical flow is produced in each chamber by two inlets, which are tangent to the swirl chamber circumference. Important changes to local and globally averaged surface Nusselt numbers, instantaneous flow structure from flow visualizations, and distributions of static pressure, total pressure, and circumferential velocity are observed throughout the swirl chambers as the Reynolds number increases. Of particular importance are increases of local surface Nusselt numbers (as well as ones globally averaged over the entire swirl chamber surface) with increasing Reynolds number. These are tied to increased advection, as well as important changes to vortex characteristics near the concave surfaces of the swirl chambers. Higher Re also give larger axial components of velocity, and increased turning of the flow from each inlet, which gives Go¨rtler vortex pair trajectories greater skewness as they are advected downstream of each inlet. [S0889-504X(00)00502-X]


Author(s):  
C. R. Hedlund ◽  
P. M. Ligrani

Local flow behavior and heat transfer results are presented from two swirl chambers, which model passages used to cool the leading edges of turbine blades in gas turbine engines. Flow results are obtained in an isothermal swirl chamber. Surface Nusselt number distributions are measured in a second swirl chamber (with a constant wall beat flux boundary condition) using infrared thermography, in conjunction with thermocouples, energy balances, and in situ calibration procedures. In both cases, Reynolds numbers Re based on inlet duct characteristics range from 6000 to about 20000. Bulk helical flow is produced in each chamber by two inlets which ore tangent to the swirl chamber circumference. Important changes to local and globally-averaged surface Nusselt numbers, instantaneous flow structure from flow visualizations, and distributions of static pressure, total pressure, and circumferential velocity are observed throughout the swirl chambers as the Reynolds number increases. Of particular importance are increases of local surface Nusselt numbers (as well as ones globally-averaged over the entire swirl chamber surface) with increasing Reynolds number. These are tiad to increased advection, as well as important changes to vortex characteristics near the concave surfaces of the swirl chambers. Higher Re also give larger axial components of velocity, and increased turning of the flow from each inlet, which gives Görtler vnrtex pair trajectories greater skewness as they are advected downstream of each inlet.


Author(s):  
Zhong Ren ◽  
Sneha Reddy Vanga ◽  
Nathan Rogers ◽  
Phil Ligrani ◽  
Keith Hollingsworth ◽  
...  

The present study provides new heat transfer data for both the surfaces of the full coverage effusion cooling plate within a double wall cooling test facility. To produce the cooling stream, a cold-side cross-flow supply for the effusion hole array is employed. Also utilized is a unique mainstream mesh heater, which provides transient thermal boundary conditions, after mainstream flow conditions are established. For the effusion cooled surface, presented are spatially-resolved distributions of surface adiabatic film cooling effectiveness, and surface heat transfer coefficients (measured using infrared thermography). For the coolant side, presented are spatially-resolved distributions of surface Nusselt numbers (measured using liquid crystal thermography). Of interest are the effects of streamwise development, blowing ratio, and Reynolds number. Streamwise hole spacing and spanwise hole spacing (normalized by effusion hole diameter) on the effusion plate are 15 and 4, respectively. Effusion hole diameter is 6.35 mm, effusion hole angle is 25 degrees, and effusion plate thickness is 3 hole diameters. Considered are overall effusion blowing ratios from 2.9 to 7.5, with subsonic, incompressible flow, and constant freestream velocity with streamwise development, for two different mainstream Reynolds numbers. For the hot side (mainstream) of the effusion film cooling test plate, results for two mainflow Reynolds numbers of about 145000 and 96000 show that the adiabatic cooling effectiveness is generally higher for the lower Reynolds number for a particular streamwise location and blowing ratio. The heat transfer coefficient is generally higher for the low Reynolds number flow. This is due to altered supply passage flow behavior, which causes a reduction in coolant lift-off of the film from the surface as coolant momentum, relative to mainstream momentum, decreases. For the coolant side of the effusion test plate, Nusselt numbers generally increase with blowing ratio, when compared at a particular streamwise location and mainflow Reynolds number.


2021 ◽  
Vol 2057 (1) ◽  
pp. 012026
Author(s):  
A V Barsukov ◽  
V V Terekhov ◽  
V I Terekhov

Abstract The results of numerical simulation of the separation flow in matrix channels by the RANS method are presented. The simulation is performed at the Reynolds number Re = 12600, determined by the mass-average velocity and the height of the channel. The distribution of the local Nusselt number is obtained for various Reynolds numbers in the range of 5÷15⋅103 and several rib angles. It is shown that the temperature distribution on the surface is highly nonuniform; in particular, the maximum heat transfer value is observed near the upper edge facets, in the vicinity of which the greatest velocity gradient is observed.


2009 ◽  
Vol 131 (11) ◽  
Author(s):  
Y.-C. Shih ◽  
J. M. Khodadadi ◽  
K.-H. Weng ◽  
A. Ahmed

The periodic state of laminar flow and heat transfer due to an insulated or isothermal rotating cylinder object in a square cavity is investigated computationally. A finite-volume-based computational methodology utilizing primitive variables is used. Various rotating objects (circle, square, and equilateral triangle) with different sizes are placed in the middle of a square cavity. A combination of a fixed computational grid and a sliding mesh was utilized for the square and triangle shapes. For the insulated and isothermal objects, the cavity is maintained as differentially heated and isothermal enclosures, respectively. Natural convection heat transfer is neglected. For a given shape of the object and a constant angular velocity, a range of rotating Reynolds numbers are covered for a Pr=5 fluid. The Reynolds numbers were selected so that the flow fields are not generally affected by the Taylor instabilities (Ta<1750). The periodic flow field, the interaction of the rotating objects with the recirculating vortices at the four corners, and the periodic channeling effect of the traversing vertices are clearly elucidated. The simulations of the dynamic flow fields were confirmed against experimental data obtained by particle image velocimetry. The corresponding thermal fields in relation to the evolving flow patterns and the skewness of the temperature contours in comparison to the conduction-only case were discussed. The skewness is observed to become more marked as the Reynolds number is lowered. Transient variations of the average Nusselt numbers of the respective systems show that for high Re numbers, a quasiperiodic behavior due to the onset of the Taylor instabilities is dominant, whereas for low Re numbers, periodicity of the system is clearly observed. Time-integrated average Nusselt numbers of the insulated and isothermal object systems were correlated with the rotational Reynolds number and shape of the object. For high Re numbers, the performance of the system is independent of the shape of the object. On the other hand, with lowering of the hydraulic diameter (i.e., bigger objects), the triangle and the circle exhibit the highest and lowest heat transfers, respectively. High intensity of the periodic channeling and not its frequency is identified as the cause of the observed enhancement.


1997 ◽  
Vol 119 (1) ◽  
pp. 61-67 ◽  
Author(s):  
J.-X. Chen ◽  
X. Gan ◽  
J. M. Owen

A superposed radial outflow of air is used to cool two disks that are rotating at equal and opposite speeds at rotational Reynolds numbers up to 1.2 × 106. One disk, which is heated up to 100°C, is instrumented with thermocouples and fluxmeters; the other disk, which is unheated, is made from transparent polycarbonate to allow the measurement of velocity using an LDA system. Measured Nusselt numbers and velocities are compared with computations made using an axisymmetric elliptic solver with a low-Reynolds-number k–ε turbulence model. Over the range of flow rates and rotational speeds tested, agreement between the computations and measurements is mainly good. As suggested by the Reynolds analogy, the Nusselt numbers for contrarotating disks increase strongly with rotational speed and weakly with flow rate; they are lower than the values obtained under equivalent conditions in a rotor–stator system.


1994 ◽  
Author(s):  
Jian-Xin Chen ◽  
Xiaopeng Gan ◽  
J. Michael Owen

This paper describes a combined experimental and computational study of the heat transfer from an electrically-heated disc rotating close to an unheated stator. A radial outflow of cooling air was used to remove heat from the disc, and local Nusselt numbers were measured, using fluxmeters at seven radial locations, for nondimensional flow rates up to C = 9680 and rotational Reynolds numbers up to Reφ = 1.2 × 106. Computations were carried out using an elliptic solver with a low-Reynolds-number k-ε turbulence model, and the agreement between the measured and computed velocities and Nusselt numbers was mainly good.


Author(s):  
Junsik Lee ◽  
Zhong Ren ◽  
Jacob Haegele ◽  
Geoff Potts ◽  
Jae Sik Jin ◽  
...  

Data which illustrate the effects of jet-to-target plate distance and Reynolds number on the heat transfer from an array of jets impinging on a flat plate are presented. Considered are Reynolds numbers Rej ranging from 8,200, to 52,000, with isentropic jet Mach numbers of approximately 0.1 to 0.2. Jet-to-target plate distances Z of 1.5D, 3.0D, 5.0D, and 8.0D are employed, where D is the impingement hole diameter. Steamwise and spanwise hole spacings are 8D. Local and spatially-averaged Nusselt numbers show strong dependence on the impingement jet Reynolds number for all situations examined. Experimental results also illustrate the dependence of local Nusselt numbers on normalized jet-to-target plate distance, especially for smaller values of this quantity. The observed variations are partially due to accumulating cross-flows produced as the jets advect downstream, as well as the interactions of the vortex structures which initially form around the jets, and then impact and interact as they advect away from stagnation points along the impingement target surface. The highest spatially-averaged Nusselt numbers are present for Z/D = 3.0 for Rej of 8,200, 20,900, and 30,000. When Rej = 52,000, spatially-averaged Nusselt numbers increase as Z/D decreases, with the highest value present at Z/D = 1.5.


2018 ◽  
Vol 29 (09) ◽  
pp. 1850079 ◽  
Author(s):  
Monireh Asadi Abchouyeh ◽  
Rasul Mohebbi ◽  
Omid Solaymani Fard

The aim of this work is to conduct numerical study of fluid flow and natural convection heat transfer by utilizing the nanofluid in a two-dimensional horizontal channel consisting of a sinusoidal obstacle by lattice Boltzmann method (LBM). The fluid in the channel is a water-based nanofluid containing Cuo nanoparticles. Thermal conductivity and nanofluid’s viscosity are calculated by Patel and Brinkman models, respectively. A wide range of parameters such as the Reynolds number ([Formula: see text]–400) and the solid volume fraction ranging ([Formula: see text]–0.05) at different non-dimensional amplitude of the wavy wall of the sinusoidal obstacle ([Formula: see text]–20) on the streamlines and temperature contours are investigated in the present study. In addition, the local and average Nusselt numbers are illustrated on lower wall of the channel. The sensitivity to the resolution and representation of the sinusoidal obstacle’s shape on flow field and heat transfer by LBM simulations are the main interest and innovation of this study. The results showed that increasing the solid volume fraction [Formula: see text] and Reynolds number Re leads to increase the average Nusselt numbers. The maximum average Nusselt number occurs when the Reynolds number and solid volume fraction are maximum and amplitude of the wavy wall is minimum. Also, by decreasing the [Formula: see text], the vortex shedding forms up at higher Reynolds number in the wake region downstream of the obstacle.


Sign in / Sign up

Export Citation Format

Share Document