BLOOD FLOW AND MACROMOLECULAR TRANSPORT IN CURVED BLOOD VESSELS

A numerical analysis of the steady/pulsatile flow and macromolecular (such as LDL and Albumin) transport in curved blood vessels was carried out. The computational results predict that the vortex of the secondary flow is time-dependent in the aortic arch. The concentration of macromolecule concentrates at the region of sharp curve, and the wall concentration at the outer part is higher than that at the inner part. Atherosclerosis and thrombosis are prone to develop in such regions with sharp flow.     

Computer Simulation of Non-Newtonian Flow and Mass Transport Through Coronary Arterial Stenosis

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｏｓｔｏｂｖｉｏｕｓｃｈａｒａｃｔｅｒｏｆａｔｈｅｒｏｓｃｌｅｒｏｓｉｓｉｓｔｈｅｌｏｃａｌｉｎｔｉｍａｌａｃｃｕｍｕｌａｔｉｏｎｏｆｌｏｗ＿ｄｅｎｓｉｔｙｌｉｐｏｐｒｏｔｅｉｎ(ＬＤＬ) ,ｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｆｉｂｅｒａｌｔｉｓｓｕｅａｎｄｔｈｅｓｔｅｎｏｓｉｓ.Ｏｎｏｎｅｈａｎｄ ,ｉｔｉｓｆｏｕｎｄｔｈａｔｔｈｅｄｉｓｅａｓｅｏｆｔｅｎｏｃｃｕｒｓａｔｔｈｅｃｏｍｐｌｅｘｇｅｏｍｅｔｒｙｒｅｇｉｏｎ ,ｓｕｃｈａｓｂｉｆｕｒｃａｔｉｏｎｚｏｎｅ,ｃｕｒｅ…     

A New Unsteady Three Dimensional Model for Macromolecular Transport and Water Filtration Across the Arterial Wall

ＩｎｔｒｏｄｕｃｔｉｏｎＡｔｈｅｒｏｇｅｎｅｓｉｓ,ｏｎｅｏｆｔｈｅｍａｊｏｒｃａｕｓｅｓｆｏｒｔｈｅｄｅａｔｈｏｆｐｅｏｐｌｅｉｎｗｅｓｔｅｒｎｃｏｕｎｔｒｉｅｓ,ｈａｓｔｗｏｐｒｏｍｉｎｅｎｔｐｒｏｐｅｒｔｉｅｓ.Ｆｉｒｓｔ,ｉｔｕｓｕａｌｌｙｔａｋｅｓｐｌａｃｅａｔｔｈｅｂｒａｎｃｈｅｓｏｒｃｕｒｖａｔｕｒｅｓｏｆｔｈｅｂｉｇａｒｔｅｒｉｅｓ,ｗｈｅｒｅｔｈｅｒｅｉｓａｎａｂｒｕｐｔｃｈａｎｇｅｏｆｂｌｏｏｄｆｌｏｗ .Ｓｅｃｏｎｄｌｙ ,ｔｈｅｒｅｉｓａｎａｃｃｕｍｕｌａｔｉｏｎｏｆｌｏｃａｌｆ…     

脉动流条件下动脉狭窄血管内脂质浓度极化现象的计算机数值模拟

在脉动流条件下,用计算机数值模拟的方法对低密度脂蛋白(LDL)在动脉狭窄血管段内的质量传输进行了研究。计算结果表明．无论是在定常流还是在脉动流条件下．LDL都将聚积于血管狭窄处峰口附近的流动分离点,LDL壁面浓度在此处最高。在脉动流条件下,LDL在血流受扰动区的聚积高于定常流的值;而且．流动分离点处LDL壁面浓度峰值覆盖的区域也宽于定常流。本文所揭示出的LDL在血管狭窄处的质量传输现象可能在动脉粥样硬化的局部性和动脉狭窄的形成中起着很重要的作用。     

Bifurcation of flow and mass transport in a curved blood vessel

A numerical analysis of flow and concentration fields of macromolecules in a, slightly curved blood vessel was carried out. Based on these results, the effect of the bifurcation of a flow on the mass transport in a curved blood vessel was discussed. The macromolecules turned out to be easier to deposit in the inner part of the curved blood vessel near the critical Dean number. Once the Dean number is higher than the critical number, the bifurcation of the flow appears. This bifurcation can prevent macromolecules from concentrating in the inner part of the curved blood vessel. This result is helpful for understanding the possible correlations between the blood dynamics and atherosclerosis. The project supported by National Natural Science Foundation of China (10002003), JSPS Postdoctoral Fellowship for Foreign Researcher and Foundation for University Teachers, the Ministry of Education     

动脉狭窄内低密度脂蛋白传输的数值研究：LDL的浓度极化现象

用计算机数值模拟的方法 ,对低密度脂蛋白 ( LDL)在动脉狭窄血管段内的质量传输进行了研究。计算结果表明 ,由于血管壁渗流的存在 ,LDL这样的脂质大分子会聚积在血管的内壁表面 ,发生一种工程上称为浓度极化的现象。LDL浓度在动脉狭窄口后的流动分离点出现峰值。该浓度峰值随雷诺数和动脉狭窄度的增加而呈逐渐下降的趋势。作者认为 ,该区域 LDL浓度的局部升高是引发动脉粥样硬化局部性和动脉狭窄产生的一个非常重要的原因。     

Effects of fluid recirculation on mass transfer from the arterial surface to flowing blood

The effect of disturbed flow on the mass transfer from arterial surface to flowing blood was studied numerically,and the results were compared with that of our previous work.The arterial wall was assumed to be viscoelastic and the blood was assumed to be incompressible and non-Newtonian fluid,which is more close to human arterial system.Numerical results indicated that the mass transfer from the arterial surface to flowing blood in regions of disturbed flow is positively related with the wall shear rates and it is significantly enhanced in regions of disturbed flow with a local minimum around the reattachment point which is higher than the average value of the downstream.Therefore,it may be implied that the accumulation of cholesterol or lipids within atheromatous plaques is not caused by the reduced efflux of cholesterol or lipids,but by the infiltration of the LDL(low-density lipoprotein) from the flowing blood to the arterial wall.     

On Recent Progress in Modelling and Simulations of Multi-scale Transfer of Mass,Momentum and Particles in Bio-medical Applications

We present a short overview of some of our most recent work that combines the mathematical modeling, advanced computer simulations and state-of-the-art experimental techniques of physical transport phenomena in various bio-medical applications. In the first example, we tackle predictions of complex blood flow patterns in the patient-specific vascular system (carotid artery bifurcation) and transfer of the so-called “bad” cholesterol (low-density lipoprotein, LDL) within the multi-layered artery wall. This two-way coupling between the blood flow and corresponding mass transfer of LDL within the artery wall is essential for predictions of regions where atherosclerosis can develop. It is demonstrated that a recently developed mathematical model, which takes into account the complex multi-layer arterial-wall structure, produced LDL profiles within the artery wall in good agreement with in-vivo experiments in rabbits, and it can be used for predictions of locations where the initial stage of development of atherosclerosis may take place. The second example includes a combination of pulsating blood flow and medical drug delivery and deposition controlled by external magnetic field gradients in the patient specific carotid artery bifurcation. The results of numerical simulations are compared with own PIV (Particle Image Velocimetry) and MRI (Magnetic Resonance Imaging) in the PDMS (silicon-based organic polymer) phantom. A very good agreement between simulations and experiments is obtained for different stages of the pulsating cycle. Application of the magnetic drug targeting resulted in an increase of up to ten fold in the efficiency of local deposition of the medical drug at desired locations. Finally, the LES (Large Eddy Simulation) of the aerosol distribution within the human respiratory system that includes up to eight bronchial generations is performed. A very good agreement between simulations and MRV (Magnetic Resonance Velocimetry) measurements is obtained. Magnetic steering of aerosols towards the left or right part of lungs proved to be possible, which can open new strategies for medical treatment of respiratory diseases.     

An investigation on near wall transport characteristics in an adiabatic upward gas–liquid two-phase slug flow

The near-wall transport characteristics, inclusive of mass transfer coefficient and wall shear stress, which have a great effect on gas–liquid two-phase flow induced internal corrosion of low alloy pipelines in vertical upward oil and gas mixing transport, have been both mechanistically and experimentally investigated in this paper. Based on the analyses on the hydrodynamic characteristics of an upward slug unit, the mass transfer in the near wall can be divided into four zones, Taylor bubble nose zone, falling liquid film zone, Taylor bubble wake zone and the remaining liquid slug zone; the wall shear stress can be divided into two zones, the positive wall shear stress zone associated with the falling liquid film and the negative wall shear stress zone associated with the liquid slug. Based on the conventional mass transfer and wall shear stress characteristics formulas of single phase liquid full-pipe turbulent flow, corrected normalized mass transfer coefficient formula and wall shear stress formula are proposed. The calculated results are in good agreement with the experimental data. The shear stress and the mass transfer coefficient in the near wall zone are increased with the increase of superficial gas velocity and decreased with the increase of superficial liquid velocity. The mass transfer coefficients in the falling liquid film zone and the wake zone of leading Taylor bubble are lager than those in the Taylor bubble nose zone and the remaining liquid slug zone, and the wall shear stress associated falling liquid film is larger than that associated the liquid slug. The mass transfer coefficient is within 10⁻³ m/s, and the wall shear stress below 10³ Pa. It can be concluded that the alternate wall shear stress due to upward gas–liquid slug flow is considered to be the major cause of the corrosion production film fatigue cracking.     

Two-phase annular flow and evaporative heat transfer in a microchannel

A physical and mathematical model has been developed to predict the two-phase flow and heat transfer in a microchannel with evaporative heat transfer. Sample solutions to the model were obtained for both constant wall temperature and constant wall heat flux conditions. Results are provided for evaporation rate, liquid film thickness, liquid and vapor phase pressure and temperature distributions. In addition to the sample calculations that were used to illustrate the transport characteristics, computations based on the current model were performed to generate results for comparisons with the experimental results of Qu and Mudawar (2004) where two different mass flow rates of the working fluid were used in the experiment. The comparisons of total pressure drops with the experimental data of Qu and Mudawar (2004) cover the wall heat flux range of 142.71-240 W/cm² with a total channel mass flux of 400.1 kg/m² s and also the wall heat flu range of 99.54-204.39 W/cm² with total channel mass flux of 401.9 kg/m² s. The calculated results from the current model match closely with those of Qu and Mudawar (2004).     

Heat transfer characteristics of developing gaseous slip-flow in rectangular microchannels with variable physical properties

The effects of variable physical properties on the flow and heat transfer characteristics of simultaneously developing slip-flow in rectangular microchannels with constant wall temperature are numerically investigated. A colocated finite-volume method is used in order to solve the mass, momentum and energy equations in their most general form. Various channel aspect ratios are studied at different Knudsen numbers. Simulations indicate that the constant physical property assumption can result in under/over-prediction of the local friction and heat transfer coefficients depending on the problem configuration. Density and thermo-physical property variations have significant effects on predicting flow and heat transfer characteristics in the developing and fully-developed regions. The degree of discrepancy varies for different cases depending on Knudsen number, aspect ratio and the temperature difference between the channel inlet and the wall. The results suggest that even low temperature differences can alter the friction and heat transfer coefficients considerably.     

Theoretical investigation of combustion characteristics in ram-jet dump combustor with side-inlet

The combustion flow of a sudden-expansioin dump combustor with injecting side-inlet is analysed using the SIMPLE-C algorithm and the Jones-Launder k-? two-equation turbulence model. The transport properties of velocity, turbulence kinetic energy, temperature and combustion efficiency as a function of the injected mass fraction and the number of side-inlet nozzles are solved in this paper. The axial velocities of the sudden-expansion dump combustor without injected side-inlet flow are solved first and found to be in good agreement with the experimental data of Moon and Rudinger. For a fixed value of the side-inlet number the wall temperature and combustion efficiency of the dump combustor are decreased when the injected mass fraction is increased. For a fixed value of the injected mass fraction the wall temperature and combustion efficiency are decreased when the number of side-inlet nozzles is increased.     

磁场和Hall电流对狭窄动脉中血液流动的影响

A micropolar model for blood simulating magnetohydrodynamic flow through a horizontally nonsymmetric but vertically symmetric artery with a mild stenosis is presented. To estimate the effect of the stenosis shape, a suitable geometry has been considered such that the horizontal shape of the stenosis can easily be changed just by varying a parameter referred to as the shape parameter. Flow parameters, such as velocity, the resistance to flow （the resistance impedance）, the wall shear stress distribution in the stenotic region, and its magnitude at the maximum height of the stenosis （stenosis throat）, have been computed for different shape parameters, the Hartmann number and the Hall parameter. This shows that the resistance to flow decreases with the increasing values of the parameter determining the stenosis shape and the Hail parameter, while it increases with the increasing Hartmann number. The wall shear stress and the shearing stress on the wall at the maximum height of the stenosis possess an inverse characteristic to the resistance to flow with respect to any given value of the Hartmann number and the Hall parameter. Finally, the effect of the Hartmann number and the Hall parameter on the horizontal velocity is examined.     

Entrainment for horizontal annular gas-liquid flow

Measurements of entrainment are presented for air and water flowing in horizontal 2.54 and 5.08 cm pipelines. After the initiation of atomization, entrainment increases with the third power of the gas velocity. At very high gas velocities a fully entrained oendition in reached for which further increases in the gas velocity do not cause a decrease in the flow rate of the wall film. Gas density bas a small effect provided comparisons are made at the same gas velocity rather than at the same mass flowrate. The results are interpreted by asauming that the rate of deposition of droplets on the wall film varies linearly with the concentration of droplets and that the rate of atomization of the wall film varies linearly with its flow rate.     

Unsteady flow with heat and mass transfer due to impingement of slot jet on an inclined plate

The unsteady laminar incompressible boundary layer flow due to a two-dimensional slot jet on a flat plate at an angle of attack has been studied. The unsteadiness in the flow field is due to the free stream velocity distribution or wall temperature (concentration) which varies with time. The governing partial differential equations in primitive variables have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The effect of the variation of the free stream velocity distribution with time is found to be more pronounced on the skin friction than on the heat or mass transfer. The Prandtl number and the variation of the wall temperature with time strongly affect the heat transfer. Similarly, the Schmidt number and the variation of the concentration at the wall with time strongly affect the mass transfer. Beyond a certain critical value of the viscous dissipation parameter, the plate gets heated instead of being cooled.     

Flow and heat/mass transfer in a wavy duct with various corrugation angles in two dimensional flow regimes

In this study, two dimensional heat/mass transfer characteristics and flow features were investigated in a rectangular wavy duct with various corrugation angles. The test duct had a width of 7.3 mm and a large aspect ratio of 7.3 to simulate two dimensional characteristics. The corrugation angles used were 100°, 115°, 130°, and 145°. Numerical analysis using the commercial code FLUENT, was used to analyze the flow features. In addition, the oil-lamp black method was used for flow visualization. Local heat/mass transfer coefficients on the corrugated walls were measured using a naphthalene sublimation technique. The Reynolds number, based on the duct hydraulic diameter, was varied from 700 to 5,000. The experimental results and numerical analysis showed interesting and detailed features in the wavy duct. Main flow impinged on upstream of a pressure wall, and the flow greatly enhanced heat/mass transfer. On a suction wall, however, flow separation and reattachment dominantly affected the heat/mass transfer characteristics on the wall. As the corrugation angle decreased (it means the duct has more sharp turn), the region of flow stagnation at the front part of the pressure wall became wider. Also, the position of flow reattachment on the suction wall moved upstream as the corrugation angle decreased. A high heat transfer rate appeared at the front part of the pressure wall due to main-flow impingement, and at the front part of the suction wall due to flow reattachment. The high heat/mass transfer region by the main-flow impingement and the circulation flow induced at a valley between the pressure and suction walls changed with the corrugation angle and the Reynolds number. As the corrugation angle decreased, the flow in the wavy duct changed to transition to turbulent flow earlier.     

Direct numerical simulation of a turbulent Couette-Poiseuille flow,part 2: Large- and very-large-scale motions

A direct numerical simulation dataset of a fully developed turbulent Couette-Poiseuille flow is analyzed to investigate the spatial organization of streamwise velocity-fluctuating u-structures on large and very large scales. Instantaneous and statistical flow fields show that negative-u structures with a small scale on a stationary bottom wall grow throughout the centerline due to the continuous positive mean shear, and they penetrate to the opposite moving wall. The development of an initial vortical structure related to negative-u structures on the bottom wall into a large-scale hairpin vortex packet with new hairpin vortices, which are created upstream and close to the wall, is consistent with the auto-generation process in a Poiseuille flow (Zhou et al., J. Fluid Mech., vol. 387, 1999, pp. 353–396). Although the initial vortical structure associated with positive-u structures on the top wall also grows toward the bottom wall, the spatial development of the structure is less coherent with weak strength due to the reduced mean shear near the top wall, resulting in less turbulent energy on the top wall. The continuous growth of the structures from a wall to the opposite wall explains the enhanced wall-normal transport of the streamwise turbulent kinetic energy near the centerline. Finally, an inspection of the time-evolving instantaneous fields and conditional averaged flow fields for the streamwise growth of a very long structure near the centerline exhibits that a streamwise concatenation of adjacent large-scale u-structures creates a very-large-scale structure near the channel centerline.     

An analytical theory of heated duct flows in supersonic combustors

One-dimensional analytical theory is developed for supersonic duct flow with variation of cross section, wall friction, heat addition, and relations between the inlet and outlet flow parameters are obtained. By introducing a selfsimilar parameter, effects of heat releasing, wall friction, and change in cross section area on the flow can be normalized and a self-similar solution of the flow equations can be found. Based on the result of self-similar solution, the sufficient and necessary condition for the occurrence of thermal choking is derived. A relation of the maximum heat addition leading to thermal choking of the duct flow is derived as functions of area ratio, wall friction, and mass addition, which is an extension of the classic Rayleigh flow theory, where the effects of wall friction and mass addition are not considered. The present work is expected to provide fundamentals for developing an integral analytical theory for ramjets and scramjets.     

Three-phase upward flow in a vertical pipe

In this paper we develop an approach to design a three-phase, gas–solid–liquid flow system that transports pneumatically scarified solid particles, including sticky ones, through a vertical pipe. The proposed system permits the introduction and maintenance of a liquid film that coats the pipe’s inner wall and acts as a lubricant that ensures sticky particles continue to move upward without permanently adhering to the pipe wall. The system’s operating conditions fall within the boundaries of the annular dispersed region on a typical flow pattern map of vertical flow of a gas–liquid mixture. High gas superficial velocities combined with low liquid superficial velocities characterize such a region. A combination of a modified one-dimensional, two-fluid annular dispersed flow model and a one-dimensional pneumatic conveying model is shown to describe this transport process satisfactorily. Solution of the combined models produces all the necessary design parameters including power requirements and superficial velocities of the two-fluid media needed to transport a given amount of solid particles. Results of model calculations are compared with rare three-phase flow data obtained prior to the development of the present model, by an independent experimental team that used the physical conditions of the present approach. Reasonable agreement justifies the use of the combined model for engineering design purposes.     

Wall Versus Centerline Polymer Injection in Turbulent Channel Flows

This experimental study compares the mean and turbulence characteristics of turbulent channel flows with polymer injection at the wall and at the centerline to assess the impact of the injection location on drag reduction. It also contrasts the drag reduction performance of a hydrolyzed polymer versus a non-ionic polymer under the same conditions. Wall injection of non-ionic and hydrolized polymers resulted in 23% and 9% larger drag reduction than corresponding centerline injection, respectively. In all cases, the polymer was structured and the presence of macromolecular polymer structures, even when concentrated mostly away from the wall, seemed to be able to affect the turbulence structure in the flow.