高温激波管化学非平衡流动数值模拟研究

为了预测氢氧定容燃烧驱动的高温激波管性能,需要准确分析激波管非定常化学非平衡流动过程.本文在破膜前的驱动段定容燃烧以及破膜后的化学非平衡流动数值模拟中,引入双时间步长方法,发展高温激波管化学非平衡流动数值模拟方法,该方法在时间上具有二阶精度.计算结果与目前存在的激波管流动解析解以及零维化学反应系统的数值解进行了比较,吻合较好.对于典型高温激波管状态,采用有限体积方法离散准一维流动Euler控制方程,并通过将流动过程和化学反应动力学过程耦合求解,获得了激波管内部的化学非平衡流动特征.     

Stabilization of flow boiling in a micro tube with air injection

Air injection as a stabilization method is evaluated for flow boiling in a micro tube. Pyrex glass tube coated by ITO film is employed as a test tube for flow visualization with water as a working fluid. Air bubble and liquid slug lengths are controlled by changing air and liquid mass velocities. Wall temperatures and inlet/outlet pressures show very large fluctuations during flow boiling without air injection. Severe reverse flow is also observed from flow visualization. On the other hand, wall temperature and inlet/outlet pressures as well as visualized flow patterns become very stable with air injection. In addition, much higher heat transfer coefficients are obtained for air injected cases. It is observed from the flow visualization that the flow becomes much stable and shows regular patterns.     

Effect of surfactant concentration on saturated flow boiling in vertical narrow annular channels

Saturated flow boiling of environmentally acceptable nonionic surfactant solutions of Alkyl (8–16) was compared to that of pure water. The concentration of surfactant solutions was in the range of 100–1000 ppm. The liquid flowed in an annular gap of 2.5 and 4.4 mm between two vertical tubes. The heat was transferred from the inner heated tube to two-phase flow in the range of mass flux from 5 to 18 kg/m² s and heat flux from 40 to 200 kW/m². Boiling curves of water were found to be heat flux and channel gap size dependent but essentially mass flux independent. An addition of surfactant to the water produced a large number of bubbles of small diameter, which, at high heat fluxes, tend to cover the entire heater surface with a vapor blanket. It was found that the heat transfer increased at low values of relative surfactant concentration C/C₀, reaches a maximum close to the value of C/C₀ = 1 (where C₀ = 300 ppm is the critical micelle concentration) and decreased with further increase in the amount of additive. The dependence of the maximal values of the relative heat transfer enhancement, obtained at the value of relative concentration of C/C₀ = 1, on the boiling number Bo may be presented as single curve for both gap sizes and the whole range of considered concentrations.     

A method of concurrent thermographic–photographic visualization of flow boiling in a minichannel

A method is developed to capture the distribution of surface temperature while simultaneously imaging the bubble motions in diabatic flow boiling in a horizontal minichannel. Liquid crystal thermography is used to obtain highly resolved surface temperature measurements on the uniformly heated upper surface of the channel. High-speed images of the flow field are acquired simultaneously and are overlaid with the thermal images. The local surface temperature and heat transfer coefficient can be analyzed with the knowledge of the nucleation site density and location, and bubble motion and size evolution. The horizontal channel is 1.2 mm high × 23 mm wide × 357 mm long, and the working fluids are Novec 649 and R-11. Optical access is through a machined glass plate which forms the bottom of the channel. The top surface is an electrically heated 76 μm-thick Hastelloy foil held in place by a water-cooled aluminum and glass frame. The heat loss resulting from this construction is computed using a conduction model in Fluent. The model is driven by temperature measurements on the foil, glass plate and aluminum frame. This model produces a corrected value for the local surface heat flux and enables the computation of the bulk fluid temperature and heat transfer coefficient along the channel. The streamwise evolution of the heat transfer coefficient for single-phase laminar flow is compared to theoretical values for a uniform-flux boundary condition. Examples of the use of the facility for visualizing subcooled two-phase flows are presented. These examples include measurements of the surface temperature distribution around active nucleation sites and the construction of boiling curves for locations along the test surface. Points on the curve can be associated with specific image sequences so that the role of mechanisms such as nucleation and the sliding of confined bubbles may be discerned.     

Effect of elastic wall motion on oscillatory flow of micropolar fluid in an annular tube

Summary Oscillatory flow of a micropolar fluid in an annular tube is investigated. The outer wall of the tube is taken to be elastic and the variation in the diameter of the elastic wall due to pulsatile nature of pressure gradient is assumed to be small. The wall motion is governed by a tube law. The nonlinear equations governing the fluid flow and the tube law are solved using perturbation analysis. The steady-streaming phenomenon due to the interaction of convected inertia with viscous effects is studied. The analysis, is carried out for zero mean flow rate. It presents the effects of the elastic nature of the wall combined with micropolar fluid parameters on the mean pressure gradient and wall shear stress for different catheter sizes and frequency parameters. It is found that the effect of micropolarity is of considerable importance for small steady-streaming Reynolds number. Also, it is observed that the relationship between mean pressure gradient and the flow rate depends on the amplitude of the diameter variation, flow rate waveforms and the phase difference between them.     

二维水平通道内流动沸腾换热的格子Boltzmann模拟

利用格子Boltzmann方法模拟二维水平通道内水的流动沸腾过程,获得不同壁面过热度下流型特点和不同因素对换热过程的影响规律。结果表明,随着壁面过热度升高,流道内流型依次经历从泡状流、弹状流到反环流的转变,平均热流密度和平均换热系数先增大后减小。入口流速降低会使流道内出现受限气泡流,核态沸腾受到抑制。提高入口流速能够有效促进气泡脱离,壁面平均换热系数随入口流速增大而增大,但增长速率有所减小。减小通道宽度有利于汽化现象发生,核态沸腾得到强化,壁面平均换热系数有所提高。     

Interstitial fluid flow: simulation of mechanical environment of cells in the interosseous membrane

In vitro experiments have shown that subtle fluid flow environment plays a significant role in living biological tissues,while there is no in vivo practical dynamical measurement of the interstitial fluid flow velocity.On the basis of a new finding that capillaries and collagen fibrils in the interosseous membrane form a parallel array,we set up a porous media model simulating the flow field with FLUENT software,studied the shear stress on interstitial cells’ surface due to the interstitial fluid flow,and analyzed the effect of flow on protein space distribution around the cells.The numerical simulation results show that the parallel nature of capillaries could lead to directional interstitial fluid flow in the direction of capillaries.Interstitial fluid flow would induce shear stress on the membrane of interstitial cells,up to 30 Pa or so,which reaches or exceeds the threshold values of cells’ biological response observed in vitro.Interstitial fluid flow would induce nonuniform spacial distribution of secretion protein of mast cells.Shear tress on cells could be affected by capillary parameters such as the distance between the adjacent capillaries,blood pressure and the permeability coefficient of capillary’s wall.The interstitial pressure and the interstitial porosity could also affect the shear stress on cells.In conclusion,numerical simulation provides an effective way for in vivo dynamic interstitial velocity research,helps to set up the vivid subtle interstitial flow environment of cells,and is beneficial to understanding the physiological functions of interstitial fluid flow.     

Numerical investigations of asymmetric flow of a micropolar fluid between two porous disks

Numerical solution is presented for the two- dimensional flow of a micropolar fluid between two porous coaxial disks of different permeability for a range of Reynolds number Re （-300≤ Re 〈 0） and permeability parameter A （1.0≤A ≤2.0）. The main flow is superimposed by the injection at the surfaces of the two disks. Von Karman＇s similarity transformations are used to reduce the governing equations of motion to a set of non-linear coupled ordinary differential equations （ODEs） in dimensionless form. An algorithm based on the finite difference method is employed to solve these ODEs and Richardson＇s extrapolation is used to obtain higher order accuracy. The results indicate that the parameters Re and A have a strong influence on the velocity and microrotation profiles, shear stresses at the disks and the position of the viscous/shear layer. The micropolar material constants cl, c2, c3 have profound effect on microrotation as compared to their effect on streamwise and axial velocity profiles. The results of micropolar fluids are compared with the results for Newtonian fluids.     

Solutions of xanthan gum/guar gum mixtures: shear rheology,porous media flow,and solids transport in annular flow

Mixtures of xanthan and guar gum in aqueous solution were studied in two flow situations: simple shear and porous media. In addition, solids transport in vertical annular flow of sand suspensions was explored. The zero shear rate viscosity of the solutions displayed a pronounced synergy: the viscosity of the mixture is higher than that of the polymer solutions in a wide range of relative concentrations of the two polymers, in agreement with previous literature. However, at relatively high shear rates, the viscosity approaches the value of the more viscous xanthan gum solutions at mass fractions of xanthan gum between 0.1 and 0.15, and the degree of synergy substantially decreases. Stress relaxation experiments in simple shear indicate that the polymer mixtures exhibit a well-defined yield stress after relaxation that is absent in solutions of pure polymers. In porous media flow experiments, a synergistic behavior mimicking the shear flow results was obtained for the polymer mixtures at low shear rates. However, at a critical shear rate, the apparent viscosity in porous media flows exceeds the shear viscosity due to the elongational nature of flow in the pores. The solids transport capacity in annular flows is well-represented by trends in shear viscosity and stress relaxation behavior. However, the lack of viscosity synergy at high shear rates limits the applicability of the mixtures as a way to improve solids suspension capacity in annular flows.     

Rotating flow of a third grade fluid in a porous space with Hall current

This study investigates the rotating magnetohydrodynamic (MHD) flow of a third-grade fluid in a porous space. Modified Darcy's law has been utilized for the flow modeling. The Hall effects are taken into consideration. The basic equations governing the flow are reduced to a highly nonlinear ordinary differential equation. This equation has been solved analytically by employing the homotopy analysis method (HAM). The effects of the various interesting parameters on the velocity distribution have been discussed.     

Peristaltic flow and heat transfer in a vertical porous annulus, with long wave approximation

In this paper, we study the interaction of peristalsis with heat transfer for the flow of a viscous fluid in a vertical porous annular region between two concentric tubes. Long wavelength approximation (that is, the wavelength of the peristaltic wave is large in comparison with the radius of the tube) is used to linearise the governing equations. Using the perturbation method, the solutions are obtained for the velocity and the temperature fields. Also, the closed form expressions are derived for the pressure-flow relationship and the heat transfer at the wall. The effect of pressure drop on flux is observed to be almost negligible for peristaltic waves of large amplitude; however, the mean flux is found to increase by 10-12% as the free convection parameter increases from 1 to 2. Also, the heat transfer at the wall is affected significantly by the amplitude of the peristaltic wave. This warrants further study on the effects of peristalsis on the flow and heat transfer characteristics.     

Two-dimensional wall temperature measurements and heat transfer enhancement for top-heated horizontal channels with flow boiling

Two-dimensional (circumferential and axial) wall temperature distributions were measured for top-heated coolant channels with internal geometries that include smooth walls, spiral fins and both twisted tape and spiral fins. Freon-71 was the working fluid. The flow regimes studied were single-phase, subcooled flow boiling, and stratified flow boiling. The inside diameter of all test sections was near 10.0 mm. Circumferentially averaged heat transfer coefficients at several axial locations were obtained for selected coolant channels for a volumetric flow rate of 4.738 x 10⁻⁵m³/s, 0.19 MPa (absolute) exit pressure, and 22.2°C inlet subcooling. Overall (averaged over the entire channel) heat transfer coefficients were compared for the various channel geometries. This comparison showed that the channel with large-pitch spiral fins had higher heat transfer coefficients at all power levels. However, the results appear to indicate that if the twist ratio (ratio of the twisted tape period to the inside diameter) is decreased, the configuration employing both fins and a twisted tape will have had greater enhancements.     

An analytical solution for steady flow of a Quemada fluid in a circular tube

An analytical solution is obtained for steady flow of Quemada-type fluids in a circular tube driven by a constant pressure gradient. Expressions are derived for velocity distribution and for volumetric flow rate as a function of pressure gradient or wall shear stress.     

Non-isothermal Poiseuille flow and its stability in a vertical annulus filled with porous medium

In this article, the non-isothermal Poiseuille flow and its stability in a vertical annulus filled with porous medium are investigated. The flow is induced by external pressure gradient and buoyancy force due to linearly varying inner wall temperature. The non-Darcy model along with Boussinesq approximation has been used. The Chebyshev spectral-collocation method has been adopted to solve the governing equations related to basic flow as well as its stability. Special attention is given to understand the effect of curvature parameter of the annular geometry on the flow, heat transfer rate and stability of the stably stratified flow. A comprehensive numerical experiment indicates that reducing gap between two concentric cylinders decreases the heat transfer rate as well as the maximum magnitude of the flow velocity. It stabilizes the flow which has been shown through stability analysis. Furthermore, appropriateness of the Forchheimer term in the momentum equation has been examined by investigating the flow regime as well as its stability in the presence and absence of Forchheimer term. Finally, it has been found from the energy analysis at critical point that the thermal-buoyant instability is the only mode of instability for the considered range of different parameters.     

一种新型螺旋内槽管的气液固三相流数值模拟研究

针对一种新型螺旋内槽管,采用先进的计算流体力学（CFD）数值模拟方法,对管内的气（天然气）-液（水）-固（水合物）三相流流动特性进行了模拟研究。模型采用欧拉-欧拉-欧拉三流体模型结合颗粒动力学的理论,考察了不同的表观速度（0.3 m/s,0.5 m/s,0.7 m/s）,水合物粒径（500 μm,750 μm,1000 μm）,气泡大小（10 μm,100 μm,1000 μm）,螺距（400 mm,800 mm）,螺纹头数（12,20）及螺纹旋向对于管内三相流动特性的影响。通过数值计算,由于气液固三相间的密度差,在螺旋内槽的作用下,水合物和天然气在管中心位置聚集,同时管壁处的含量减小。流体表观流速和气泡越大,壁面处的水合物和天然气的体积分数越小;由于天然气的密度小于水合物和水的密度,天然气更多集中在管中心,越靠近管壁含量越少;颗粒的粒径越大,壁面处的水合物含量越少,而对于天然气的分布则影响不大;螺距越小,螺纹头数越多,螺旋流强度越大,气液固三相分离效果越好,壁面处的水合物和天然气的含量越小;同时,螺纹旋向的改变对于三相的分离效果影响较小。     

Numerical simulation of streamwise fluidelastic instability of tube bundles subjected to two-phase cross flow

This work aims to develop and validate a numerical model to simulate the flow-structure interaction in tube bundles subjected to two-phase flow. The model utilizes a mixture multiphase module in which a drift flux formulation is used to account for the slip between the phases. Two methods of numerical flow-structure interaction are used to predict the onset of fluidelastic instability (FEI) in the streamwise direction for a two-phase air–water flow mixture in parallel triangular tube bundles. These models are the hybrid analytical-flow field model and the direct numerical flow/structure coupling model. This work investigates the effects of void fractions in the range of 20% to 80% and several pitch-to-diameter ratios (P/D) in the range of 1.3 to 1.7. The results of the fluidelastic forces and the stability threshold are validated against the experimental data available in the literature and show an excellent agreement. The streamwise FEI threshold shows a significant dependency on the pitch-to-diameter ratio while the void fraction exhibits a lesser effect. Generally, the stability threshold increases as the pitch-to-diameter ratio increases. The model that was developed paves the way for devising of more reliable prediction tools for FEI in steam generators.     

Flow pattern and flow characteristics for counter-current two-phase flow in a vertical round tube with wire-coil inserts

Flow pattern, void fraction and slug rise velocity on counter-current two-phase flow in a vertical round tube with wire-coil inserts are experimentally studied. Flow pattern and slug rise velocity are measured visually with a video camera. The void fraction is measured by the quick-closing valve method. Four kinds of coils with different coil pitches and coil diameters are used as inserts. The presence of wire-coil inserts induces disturbance into gas and liquid flows so that the shape and motion of gas slug or bubbles in a wire-coil inserted tube are quite different from those observed in a smooth tube without insert. The bubbly flow occurs in the low gas superficial velocity region in the wire-coil inserted tube, while the slug or churn/annular flow only appears in the smooth tube without insert over the all test range. The measured slug rise velocity in the wire-coil inserted tube is higher than that in the smooth tube. With modified mean flow velocity calculated with core area, the slug rise velocity in wire-coil tube inserted is in good agreement with Nicklin's correlation. The void fraction in a wire-coil inserted tube is lower than that in a smooth tube in the range of high gas superficial velocities. By introducing a simple assumption on considering the effective flowing area, the measured void fractions in a wire-coil inserted tube are in relatively good agreement with the predicted result based on the drift flux model proposed by others with the correlation for slug rise velocity given by others when the coil pitch is dense.     

Numerical simulation of blood flow and interstitial fluid pressure in solid tumor microcirculation based on tumor-induced angiogenesis

A coupled intravascular–transvascular–interstitial fluid flow model is developed to study the distributions of blood flow and interstitial fluid pressure in solid tumor microcirculation based on a tumor-induced microvascular network. This is generated from a 2D nine-point discrete mathematical model of tumor angiogenesis and contains two parent vessels. Blood flow through the microvascular network and interstitial fluid flow in tumor tissues are performed by the extended Poiseuille’s law and Darcy’s law, respectively, transvascular flow is described by Starling’s law; effects of the vascular permeability and the interstitial hydraulic conductivity are also considered. The simulation results predict the heterogeneous blood supply, interstitial hypertension and low convection on the inside of the tumor, which are consistent with physiological observed facts. These results may provide beneficial information for anti-angiogenesis treatment of tumor and further clinical research. The project supported by the National Natural Science Foundation of China (10372026).