Convective and radiative heat transfer coefficients for individual human body segments

 Human thermal physiological and comfort models will soon be able to simulate both transient and spatial inhomogeneities in the thermal environment. With this increasing detail comes the need for anatomically specific convective and radiative heat transfer coefficients for the human body. The present study used an articulated thermal manikin with 16 body segments (head, chest, back, upper arms, forearms, hands, pelvis, upper legs, lower legs, feet) to generate radiative heat transfer coefficients as well as natural- and forced-mode convective coefficients. The tests were conducted across a range of wind speeds from still air to 5.0 m/s, representing atmospheric conditions typical of both indoors and outdoors. Both standing and seated postures were investigated, as were eight different wind azimuth angles. The radiative heat transfer coefficient measured for the whole-body was 4.5 W/m² per K for both the seated and standing cases, closely matching the generally accepted whole-body value of 4.7 W/m² per K. Similarly, the whole-body natural convection coefficient for the manikin fell within the mid-range of previously published values at 3.4 and 3.3 W/m² per K when standing and seated respectively. In the forced convective regime, heat transfer coefficients were higher for hands, feet and peripheral limbs compared to the central torso region. Wind direction had little effect on convective heat transfers from individual body segments. A general-purpose forced convection equation suitable for application to both seated and standing postures indoors was h _c=10.3v ^0.6 for the whole-body. Similar equations were generated for individual body segments in both seated and standing postures. Received: 21 May 1996/Accepted: 27 November 1996     

Computational fluid dynamics analysis of cyclist aerodynamics: Performance of different turbulence-modelling and boundary-layer modelling approaches

This study aims at assessing the accuracy of computational fluid dynamics (CFD) for applications in sports aerodynamics, for example for drag predictions of swimmers, cyclists or skiers, by evaluating the applied numerical modelling techniques by means of detailed validation experiments. In this study, a wind-tunnel experiment on a scale model of a cyclist (scale 1:2) is presented. Apart from three-component forces and moments, also high-resolution surface pressure measurements on the scale model’s surface, i.e. at 115 locations, are performed to provide detailed information on the flow field. These data are used to compare the performance of different turbulence-modelling techniques, such as steady Reynolds-averaged Navier–Stokes (RANS), with several k–ε and k–ω turbulence models, and unsteady large-eddy simulation (LES), and also boundary-layer modelling techniques, namely wall functions and low-Reynolds number modelling (LRNM). The commercial CFD code Fluent 6.3 is used for the simulations. The RANS shear-stress transport (SST) k–ω model shows the best overall performance, followed by the more computationally expensive LES. Furthermore, LRNM is clearly preferred over wall functions to model the boundary layer. This study showed that there are more accurate alternatives for evaluating flow around bluff bodies with CFD than the standard k–ε model combined with wall functions, which is often used in CFD studies in sports.     

Aerodynamic study of different cyclist positions: CFD analysis and full-scale wind-tunnel tests

Three different cyclist positions were evaluated with Computational Fluid Dynamics (CFD) and wind-tunnel experiments were used to provide reliable data to evaluate the accuracy of the CFD simulations. Specific features of this study are: (1) both steady Reynolds-averaged Navier–Stokes (RANS) and unsteady flow modelling, with more advanced turbulence modelling techniques (Large-Eddy Simulation – LES), were evaluated; (2) the boundary layer on the cyclist’s surface was resolved entirely with low-Reynolds number modelling, instead of modelling it with wall functions; (3) apart from drag measurements, also surface pressure measurements on the cyclist’s body were performed in the wind-tunnel experiment, which provided the basis for a more detailed evaluation of the predicted flow field by CFD. The results show that the simulated and measured drag areas differed about 11% (RANS) and 7% (LES), which is considered to be a close agreement in CFD studies. A fair agreement with wind-tunnel data was obtained for the predicted surface pressures, especially with LES. Despite the higher accuracy of LES, its much higher computational cost could make RANS more attractive for practical use in some situations. CFD is found to be a valuable tool to evaluate the drag of different cyclist positions and to investigate the influence of small adjustments in the cyclist’s position. A strong advantage of CFD is that detailed flow field information is obtained, which cannot easily be obtained from wind-tunnel tests. This detailed information allows more insight in the causes of the drag force and provides better guidance for position improvements.     

Experimental study on convective heat transfer coefficient of the human body

1. 1. The convective heat transfer coefficient of the human body is essential to predict convective heat loss from the body.

2. 2. The object of this paper is to calculate the convective heat transfer coefficient of the human body using heat flow meters and to estimate the thermally equivalent sphere and cylinder to the human body.

3. 3. The experimental formulae of the convective heat transfer coefficient for the whole body were obtained by regression analysis for natural, forced and mixed convection.

4. 4. Diameters of the thermally equivalent sphere and cylinder of the human body were calculated as 12.9 and 12.2 cm, respectively.

Mechanism of heat transfer in different regions of human body

The processes of heat transfer in a human body were studied with the use of a mathematical model. It has been shown that only conductive or only convective heat transfer may occur in different body areas. The rate of blood-mediated heat transfer in the presence of blood circulation is many times higher than heat transfer due to temperature gradient; therefore, the convective process prevails over the conductive process. The body core contains a variety of blood vessels, and the bulk of blood concentrates there in the norm. Hence, heat transfer in it is mainly convective. In surface tissues, where the rate of blood circulation is lower and the vasculature has certain specific features, heat transfer is mainly conductive. Hence, the core and surface tissues are absolutely different body zones in terms of heat transfer.     

高粘发酵体系不同搅拌桨的CFD模拟及发酵验证

搅拌桨是高好氧高黏度微生物发酵实现高效反应必不可少的因素之一,不同搅拌桨组合对发酵过程的影响十分重要。威兰胶是由产碱杆菌在高耗氧高粘度发酵体系下合成的胞外微生物多糖,广泛应用于水泥、石油、油墨、食品等行业中。本研究借助于计算流体力学(Computational fluid dynamics,CFD)的方法,以威兰胶发酵液体系为研究体系,研究了6种不同搅拌桨组合在反应器内流体速率分布、剪切速率、和气含率等参数。将模拟效果较好的3种组合用于威兰胶发酵过程。研究表明MB-4-6搅拌桨组合对改善发酵罐内部的溶氧及流场分布效果最明显,威兰胶产量水平提高了13%。同时在该组合下威兰胶的产品粘度得到有效提高。     

Development and validation of computational fluid dynamics models for prediction of heat transfer and thermal microenvironments of corals

We present Computational Fluid Dynamics (CFD) models of the coupled dynamics of water flow, heat transfer and irradiance in and around corals to predict temperatures experienced by corals. These models were validated against controlled laboratory experiments, under constant and transient irradiance, for hemispherical and branching corals. Our CFD models agree very well with experimental studies. A linear relationship between irradiance and coral surface warming was evident in both the simulation and experimental result agreeing with heat transfer theory. However, CFD models for the steady state simulation produced a better fit to the linear relationship than the experimental data, likely due to experimental error in the empirical measurements. The consistency of our modelling results with experimental observations demonstrates the applicability of CFD simulations, such as the models developed here, to coral bleaching studies. A study of the influence of coral skeletal porosity and skeletal bulk density on surface warming was also undertaken, demonstrating boundary layer behaviour, and interstitial flow magnitude and temperature profiles in coral cross sections. Our models compliment recent studies showing systematic changes in these parameters in some coral colonies and have utility in the prediction of coral bleaching.     

The heat transfer analysis of nanoparticle heat source in alanine tissue by molecular dynamics

The purpose of this study is to simulate the heat transfer problem when the 3-D Alanine tissue is heated by the gold nanoparticle in the field of molecular dynamics. In this paper, the Alanine molecule is adopted and its parameters are available in the GROMACS protein data bank. A computing algorithm is developed to evaluate the heat transfer phenomena in the nano-scale biological system based on the molecular dynamics and the protein data bank. The value of the thermal conductivity of Alanine is calculated from the autocorrelation function of the Green-Kubo formalism and this result has a roughly approximation with the bulk thermal conductivity reported by experimental data . Two kinds of problems are investigated in the paper. One is the Alanine tissue heated by the constant heat source and the other is by the time-varying heat source. The numerical results show that a temperature jump exists around the source and the temperature profiles drop to the environmental temperature within a very short distance. It concludes that only a small region around the nano-scale heat source is affected by the heated process. Therefore, the results of the nanoparticle-heated method could be applied to the clinical therapy of tumor, and the normal cells are destroyed only within a smaller region than those of chemotherapy or surgery.     

Effect of nanostructure on heat transfer between fluid and copper plate: a molecular dynamics simulation study

A molecular simulation is developed to study the effect of surface nanostructures on nanoscale flows. Based on this method, particles equation of motion is solved through the Verlet algorithm. Meanwhile, a physically sound method is applied to control the momentum and temperature of the simulation box. By adding an external force on the top copper plate according to the velocity difference between on-the-fly and desired velocities, simulations on convection of argon flows between two solid walls are performed. The top wall, which holds a higher temperature, moves at a constant velocity relative to bottom one along with the streamwise direction. These simulation results show that the nanostructures particularly affect fluid density oscillations adjacent to solid wall and nanostructures. In addition, these nanostructures also have significant effects on temperature and velocity distributions in simulation system.     

Airflow behavior changes in upper airway caused by different head and neck positions: Comparison by computational fluid dynamics

The feasibility of computational fluid dynamics (CFD) to evaluate airflow characteristics in different head and neck positions has not been established. This study compared the changes in volume and airflow behavior of the upper airway by CFD simulation to predict the influence of anatomical and physiological airway changes due to different head–neck positions on mechanical ventilation. One awake volunteer with no risk of difficult airway underwent computed tomography in neutral position, extension position (both head and neck extended), and sniffing position (head extended and neck flexed). Three-dimensional airway models of the upper airway were reconstructed. The total volume (V) and narrowest area (A_min) of the airway models were measured. CFD simulation with an Spalart–Allmaras model was performed to characterize airflow behavior in neutral, extension, and sniffing positions of closed-mouth and open-mouth ventilation. The comparison result for V was neutral <extension≈sniffing, and for A_min was neutral<extension<sniffing. A_min in sniffing position was nearly 3.0 times that in neutral position and 1.7 times that in extension position. The pressure drop and velocity increasing were more obvious in neutral than sniffing or extension position at the same airflow rate. In sniffing position, pressure differences decreased and velocity remained almost constant. Recirculation airflow was generated near the subglottic region in neutral and extension positions. Sniffing position improves airway patency by increasing airway volume and decreasing airway resistance, suggesting that sniffing position may be the optimal choice for mask ventilation.     

Computational fluid dynamics modeling of mass transfer behavior in a bioreactor for hairy root culture. I. Model development and experimental validation

A two-dimensional axisymmetric computational fluid dynamics (CFD) model based on a porous media model and a discrete population balance model was established to investigate the hydrodynamics and mass transfer behavior in an airlift bioreactor for hairy root culture.During the hairy root culture of Echinacea purpurea, liquid and gas velocity, gas holdup, mass transfer rate, as well as oxygen concentration distribution in the airlift bioreactor were simulated by this CFD model. Simulative results indicated that liquid flow and turbulence played a dominant role in oxygen mass transfer in the growth domain of the hairy root culture. The dissolved oxygen concentration in the hairy root clump increased from the bottom to the top of the bioreactor cultured with the hairy roots, which was verified by the experimental detection of dissolved oxygen concentration in the hairy root clump. This methodology provided insight understanding on the complex system of hairy root culture and will help to eventually guide the bioreactor design and process intensification of large-scale hairy root culture.     

Computational fluid dynamics modeling and analysis of the effect of 3-D distortion of the human aortic arch

An idealized CFD model and a realistic one were used to investigate the effect of the 3-D distortion of the aortic arch on the blood flow and its pathophysiological significance with respect to the pathogenesis of the aortic aneurysm. From the results of the flow simulations, the distortion of the centerline of the pipe was shown to affect significantly the flow structure. A right-handed vortex at the descending arch, and a left-handed one at the end of the arch tended to develop in the realistic model. But the secondary flow did not become a single helix. The top of the arch was the region where complex spatial and temporal WSS distributed. It was also observed that the direction of WSS had a significant circumferential component at the top of the arch.     

Computational fluid dynamics analysis of surgical adjustment of left ventricular assist device implantation to minimise stroke risk

Background. Currently, mechanical support is the most promising alternative to cardiac transplantation. Ventricular assist devices (VADs) were originally used to provide mechanical circulatory support in patients awaiting planned heart transplantation (‘bridge-to-transplantation’ therapy). The success of short-term bridge devices led to clinical trials evaluating the clinical suitability of long-term support (‘destination’ therapy) with left ventricular assist devices (LVADs). The first larger scale, randomised trial that tested long-term support with an LVAD reported a 44% reduction in the risk of stroke or death in patients with an LVAD. In spite of the success of LVADs as bridge-to-transplantation and long-term support, patients managed by these devices are still at risk of several adverse events. The most devastating complication is caused by embolisation of thrombi formed within the LVAD or inside the heart into the brain. Prevention of thrombi formation is attempted through anticoagulation management and by improving LVADs design; however, there is still significant occurrence of thromboembolic events in patients. Investigators have reported that the incidence of thromboembolic cerebral events ranges from 14% to 47% over a period of 6–12 months.

Methods and approach. An alternative method to reduce the incidence of cerebral embolisation is proposed by the co-authors, and the hypothesis is that it is possible to minimise the number of thrombi flowing into the carotid and vertebral arteries by an optimal placement of the LVAD outflow conduit, with or without the addition of aortic bypass connecting the ascending aorta and the innominate artery (IA), or left carotid artery. This paper presents the computational fluid dynamics (CFD) analysis of the aortic arch haemodynamics using a representative geometry of the human aortic arch with or without an alternative aortic bypass. In order to study the trajectory of the thrombi within the aortic arch bed, the CFD code, Fluent 6.3, is utilised to resolve the flow field and to solve the Lagrangian particle tracking of thrombi released randomly at the inlet of the LVAD cannula.

Results. Results are presented for simulations of thrombi in the range of 2–5 mm. The percentage of individual diameter as well as aggregate diameter thrombi flowing to the carotid and vertebral arteries as a function of LVAD conduit placement and aortic bypass implantation is reported. The influence of the LVAD conduit implantation and bypass reveals a nearly 50% variation in predicted cerebral embolism rates.

Conclusions. The adjustment of the location of the anastomosis of the LVAD outflow cannula as well as its angle of incidence plays a significant role in the level of thromboembolisms. By proper adjustment in this CFD study of a synthetic model of an aortic arch bed, we found that nearly a 50% reduction in cerebral embolism could be achieved for a configuration consisting of a shallow angle of implantation over a baseline normal incidence of the LVAD cannula. Within the limitations of our model, we have established that the LVAD implantation geometry is an important factor and should be taken into consideration when implanting an LVAD. It is possible that other parameters such as distance of the LVAD outflow cannula to the root of the IA could affect the thrombi embolisation probabilities. However, the results of this study suggest that the risk of stroke may be significantly reduced by as much as 50% by tailoring the VAD implantation by a simple surgical manoeuvre. The results of this line of research may ultimately lead to techniques that can be used to estimate the optimal LVAD configuration in a patient-specific manner by pre-operative imaging.     

Computational fluid dynamics as a tool for testing functional and ecological hypotheses in fossil taxa

Computational fluid dynamics is a method for simulating fluid flows that has been widely used in engineering for decades, and which also has applications for studying function and ecology in fossil taxa. However, despite the possible benefits of this approach, computational fluid dynamics has been used only rarely in palaeontology to date. The theoretical basis underlying the technique is outlined and the main steps involved in carrying out computer simulations of fluid flows are detailed. I also describe previous studies that have applied the method to fossils and discuss their potential for informing future research directions in palaeontology. Computational fluid dynamics can enable large‐scale comparative analyses, as well as exacting tests of hypotheses related to the function and ecology of ancient organisms. In this way, it could transform our understanding of many extinct fossil groups.     

A molecular dynamics boundary condition for heat exchange between walls and a fluid

In molecular dynamics simulations of heat transfer in micro channels, a lot of computation time is used when the wall molecules are explicitly simulated. To save computation time, implicit boundary conditions, such as the Maxwell conditions, can be used. With these boundary conditions, heat transfer is still a problem. In this work, we derive a new boundary condition based on a vibrating potential wall. The heat-transfer properties of the new boundary condition are shown to be comparable with those of the explicit wall. The computation time needed for the implicit boundary condition is very small compared with that needed for the explicit simulation.     

Computational fluid dynamics modeling of an inverted frustoconical shaking bioreactor for mammalian cell suspension culture

We previously developed an inverted frustoconical shaking bioreactor (IFSB) which had high mammalian cell culture performance when compared with a mechanically stirred tank reactor (STR) or a flat-bottom shaking bioreactor (FBSB). Here, we determined the mixing time (t) and volumetric oxygen transfer coefficient (k _La) of this IFSB at various speeds, and simulated the fluid hydrodynamics, including the shear stress and specific surface area, by computational fluid dynamics. The shortest mixing time was observed in a STR. The maximum kLa value of 12/h was achieved in the IFSB at an aeration rate of 4 L/h, demonstrating that our IFSB has enhanced oxygen transfer capabilities needed to meet the demands of mammalian cells. Simulation studies revealed a 3% greater specific surface area and a 21% lower shear strain in the IFSB compared to an FBSB under the same conditions. Additionally, the conical angle of the vessel, which significantly affected cell growth and recombinant protein production, was tested here. We conclude that, compared to the STR and FBSB, the IFSB has an increased liquid surface area for oxygen uptake and exhaust CO₂ stripping, an enhanced k _La for cell robust growth to a high cell density, and a lower shear stress to alleviate cell damage.