The bubble-slug flow pattern transition and instabilities of void fraction waves

The propagation of spontaneous void fraction disturbances in a nitrogen-water flow has been studied through the statistical analysis of conductivity probe signals, for void fractions ranging from 0.1 to 0.5 and including the bubble-slug transition. The power spectral density function and the standard deviation of the void fraction have been computed for each probe, as well as the system phase factor (related to the wave velocity), the coherence function and the system gain factor between each pair of consecutive probes as functions of frequency. For bubble flow, the results are compatible with the results obtained by other authors. The transition from bubble to slug flow is associated with void fraction wave instabilities. Two kinds of instabilities seem to occur simultaneously: amplitude increase (system gain factor > 1) and wave-breaking.     

Measurement accuracy of a mono-fiber optical probe in a bubbly flow

The measurement accuracy of a mono-fiber optical probe is studied experimentally using isolated bubbles rising freely in a still liquid. The dwell time of the probe tip within the gas phase, which is obtained from both the optical probe signal and high-speed visualization, is compared with the value expected for a non-perturbed bubble. The difference originates mainly from the intrusive nature of the optical probe, which modifies the bubble behavior when it comes into contact with the probe tip. This interaction increases the dwell time if the bubble is pierced by the probe near its pole, and shortens it for piercing near the equator. The mean dwell time, obtained by averaging for various piercing locations, is shortened and the local void fraction indicated by the probe is thus underestimated. It is shown that the void fraction error can be correlated with a modified Weber number, and this correlation is helpful for sensor selection and for uncertainty estimate. In addition, the distribution of gas dwell time usually differs from the response expected for an ideal probe. This deviation results from the dependence of the dwell time error on the piercing location. The dwell time distribution can be used to infer the dependence of the dwell time on the piercing location. Finally, the deformation of long fibers during the bubble-probe interaction significantly increases the measurement error. Observed results are consistent with data of Andreotti (2009), which were measured in an airlift flow, suggesting that present results are applicable also to the case of moving liquid. Conclusions of this study could be applied also to conductivity probes or more generally to the interaction of a bubble with any kind of thin, intrusive sensor or fiber.     

Modeling and Simulation of Local Thermal Non-equilibrium Effects in Porous Media with Small Thermal Conductivity

The two-equation model in porous media can describe the local thermal non-equilibrium (LTNE) effects between fluid and solid at REV scale, with the temperature differences in a solid particle neglected. A multi-scale model has been proposed in this study. In the model, the temperature differences in a solid particle are considered by the coupling of the fluid energy equation at REV scale with the heat conduction equation of a solid particle at pore scale. The experiments were conducted to verify the model and numerical strategy. The multi-scale model is more suitable than the two-equation model to predict the LTNE effects in porous media with small thermal conductivity. The effects of particle diameter, mass flow rate, and solid material on the LTNE effects have been investigated numerically when cryogenic nitrogen flows through the porous bed with small thermal conductivity. The results indicate that the temperature difference between solid center and fluid has the same trend at different particle diameters and mass flow rates, while the time to reach the local thermal equilibrium is affected by solid diameter dramatically. The results also show that the temperature difference between solid center and surface is much greater than that between solid surface and fluid. The values of \( \rho {\text{c}} \) for different materials have important influence on the time to reach the local thermal equilibrium between solid and fluid.     

Effect of temperature gradient within a solid particle on the rotation and oscillation modes in solid-dispersed two-phase flows

Liquid–solid two-phase flow with heat transfer is simulated, and the effect of temperature gradient within a solid particle on the particle behaviour and heat transfer is studied. The interaction between fluid and particles is considered with our original immersed solid approach on a rectangular grid system. The local heat flux at the fluid–solid interface is described with an anisotropic heat conductivity matrix, and the governing equation of temperature is time-updated with an implicit treatment for the diffusion term. The method is applied to a 2-D natural convection flow of a relatively low Rayleigh number including multiple particles. Heat transfer and particle behaviours are studied for different solid heat conductivities (ratio to the fluid conductivity ranging between 10⁻³ and 10³) and solid volume fractions. Under a condition of relatively low heat conductivity ratio, the particles show a simple circulating flow. By increasing the heat conductivity ratio, a transition of the particulate flow is observed to oscillation mode around the domain centre due to the buoyancy force as a restitution force. The oscillation period is found to vary with the heat conductivity ratio, and it is related to the time scales for the heat transfer via fluid and solid.     

A Directed Percolation Model for Clogging in a Porous Medium with Small Inhomogeneities

We propose a microscopic model based on directed percolation for the process of mechanical clogging of a porous medium by particles suspended in a fluid flow. Under appropriate conditions the deposited particles may form fractal clusters. A criterion for the occurrence of fractal clogging is presented. It links together the particle size and the pore size distribution. The effect of microscopic inhomogeneities is studied inside and outside the critical region using Monte Carlo calculations in two dimensions. The critical exponents remain unchanged because the perturbation induced by these inhomogeneities is irrelevant. The percolation threshold is found to shift to higher values almost linearly with increasing size of obstacles. For size distributed obstacles the arithmetic mean of the distribution is the only significant parameter which determines the shift. Type and broadness of the distribution have no influence. Also the percolation probability depends only on the mean even outside the critical region for all values of the occupation probability. Occupying the same fraction of the porous matrix, large obstacles cause more particles to deposit than small ones.     

Drift-flux correlation for gas-liquid two-phase flow in a horizontal pipe

A drift-flux correlation has been often used to predict void fraction of gas-liquid two-phase flow in a horizontal channel due to its simplicity and practicality. The drift-flux correlation includes two important drift-flux parameters, namely, the distribution parameter and void-fraction-weighted-mean drift velocity. In this study, an extensive literature survey for horizontal two-phase flow is conducted to establish void fraction database and to acquire existing drift-flux correlations. A total of 566 data is collected from 12 data sources and 4 flow-regime-dependent and 1 flow-regime-independent drift-flux correlations are identified. The predictive capability of the existing drift-flux correlations is assessed using the collected data. It is pointed out that the drift velocity determined by a regression analysis may include a significant error due to a compensation error between distribution parameter and drift velocity. In this study, a simple flow-regime-independent drift-flux correlation is developed. In the modeling approach, the void-fraction-weighted mean drift velocity is approximated to be 0 m/s, whereas the distribution parameter is given as a simple function of the ratio of non-dimensional superficial gas velocity to non-dimensional mixture volumetric flux. The newly developed correlation shows an excellent predictive capability of void fraction for horizontal two-phase flow. Mean absolute error (or bias), standard deviation (random error), mean relative deviation and mean absolute relative deviation of the correlation are 0.0487, 0.0985, 0.0758 and 0.206, respectively. The prediction accuracy of the correlation is similar to the correlation of Chexal et al. (1991), which was formulated based on the drift-flux parameters by means of many cascading constitutive relationships with numerous empirical parameters.     

The Effect of Strong Heterogeneity on the Onset of Convection in a Porous Medium: 2D/3D Localization and Spatially Correlated Random Permeability Fields

The effect of strong heterogeneity on the onset of convection induced by a vertical density gradient in a saturated porous medium governed by Darcy’s law is investigated. The general case, where there is heterogeneity in both the vertical and horizontal directions, and where there is heterogeneity in permeability, thermal conductivity, and applied temperature gradient, is considered. A computer package has been developed to implement an algorithm giving a criterion for instability, and this is now employed to investigate the case where there is two-dimensional variation in a horizontal plane and the case where the variation is generated by a log-normal distribution. In the latter case, spatially correlated fields with known stochastic properties are generated, and the results are analyzed in a statistical framework. We now test cases that are representative of natural, field-scale, geologic conditions—both in terms of the correlated structure and the much larger standard deviation of the permeability distribution.     

Numerical simulations of high pressure carbon dioxide fluid fluidized beds

Hydrodynamics of carbon dioxide fluid-particle mixtures are predicted using a low density ratio-based kinetic theory of granular flow in high pressure carbon dioxide fluid fluidized beds. A coexistence of particle waves and particle aggregates exists along bed height. The threshold to identify the occurrence of particle aggregates is suggested based on standard deviation of solid volume fractions in aggregative fluidization. The existence time fractions and frequencies of particle aggregates are predicted along axial direction. The effect of carbon dioxide fluid temperature and pressure on volume fraction and velocity distributions are analyzed at different inlet carbon dioxide velocities and particle densities in high pressure carbon dioxide fluidized beds. Simulated results indicate that the carbon dioxide-particles fluidization transits from particulate to aggregative states with the increase of inlet carbon dioxide velocities. The computed fluid volume fractions and heterogeneity indexes are close to the measurements in a high pressure carbon dioxide fluidized bed.     

Study of two-phase flows in reduced gravity using ground based experiments

Experimental studies have been carried out to support the development of a framework of the two-fluid model along with an interfacial area transport equation applicable to reduced gravity two-phase flows. The experimental study simulates the reduced gravity condition in ground based facilities by using two immiscible liquids of similar density namely, water as the continuous phase and Therminol 59^® as the dispersed phase. We have acquired a total of eleven data sets in the bubbly flow and bubbly to slug flow transition regimes. These flow conditions have area-averaged void (volume) fractions ranging from 3 to 30% and channel Reynolds number for the continuous phase between 2,900 and 8,800. Flow visualization has been performed and a flow regime map developed which is compared with relevant bubbly to slug flow regime transition criteria. The comparison shows that the transition boundary is well predicted by the criterion based on critical void fraction. The value of the critical void fraction at transition was experimentally determined to be approximately 25%. In addition, important two-phase flow local parameters, including the void fraction, interfacial area concentration, droplet number frequency and droplet velocity, have been acquired at two axial locations using state-of-the-art multi-sensor conductivity probe. The radial profiles and axial development of the two-phase flow parameters show that the coalescence mechanism is enhanced by either increasing the continuous or dispersed phase Reynolds number. Evidence of turbulence induced particle interaction mechanism is highlighted. The data presented in this paper clearly show the marked differences in terms of bubble (droplet) size, phase distribution and phase interaction in two-phase flow between normal and reduced gravity conditions.     

Experimental study of enhanced heat transfer by addition of CuO nanoparticle

An energy storage system has been designed to study the thermal characteristics of paraffin wax with an embedded nano size copper oxide (CuO) particle. This paper presents studies conducted on phase transition times, heat fraction as well as heat transfer characteristics of paraffin wax as phase change material (PCM) embedded with CuO nanoparticles. 40?nm mean size CuO particles of 2, 5 and 10% by weight were dispersed in PCM for this study. Experiments were performed on a heat exchanger with 1.5–10?l/min of heat transfer fluid (HTF) flow. Time-based variations of the temperature distributions are revealed from the results of observations of melting and solidification curves. The results strongly suggested that the thermal conductivity enhances 6, 6.7 and 7.8% in liquid state and in dynamic viscosity it enhances by 5, 14 and 30% with increasing mass fraction of the CNEPs. The thermal conductivity ratio of the composites can be augmented by a factor up to 1.3. The heat transfer coefficient during solidification increased about 78% for the maximum flow rate. The analysis of experimental results reveals that the addition of copper oxide nanoparticles to the paraffin wax enhances both the conduction and natural convection very effectively in composites and in paraffin wax. The paraffin wax-based composites have great potential for energy storage applications like industrial waste heat recovery, solar thermal applications and solar based dynamic space power generation with optimal fraction of copper oxide nanoparticles.     

Free Convection Heat Transfer From A Sphere In A Porous Medium Using A Thermal Non-equilibrium Model

This work studies the free convection heat transfer from a sphere with constant wall temperature embedded in a fluid-saturated porous medium using a thermal non-equilibrium model. The governing equations are transformed into boundary-layer partial differential equations by the coordinate transform, and the obtained governing equations are then solved by the cubic spline collocation method. The temperature distributions for fluid and solid phases are shown for different values of the porosity scaled thermal conductivity ratio, the interphase heat transfer parameter, and the streamwise coordinate. The effects of the porosity scaled thermal conductivity ratio and the interphase heat transfer parameter between solid and fluid phases on the local Nusselt numbers for fluid and solid phases are examined. Results show the local Nusset number for the porous medium can be increased by increasing the porosity scaled thermal conductivity ratio. Moreover, the thermal non-equilibrium effect is more significant for low values of the porosity scaled thermal conductivity ratio or the interphase heat transfer parameter.     

基于人工神经网络的非均匀固相应力模型

最小多尺度理论EMMS已经被引入多相质点网格法MP-PIC中,建立了非均匀EMMS固相应力模型.但现有的非均匀固相应力模型计算中,中间步骤繁琐且花费时间长.采用人工拟合的方式能获得非均匀固相应力表达式,但需要人为确定拟合变量和拟合函数,且针对于非均匀固相应力这种高度非线性函数所得到的拟合精度不高、用于MP-PIC模拟的结果相比原EMMS固相应力模型结果存在偏差.针对上述问题,本文提出通过机器学习的方法,规避对固相体积分数的局部分布情况的表征,并提出和建立能考虑颗粒浓度详细分布的人工神经网络ANN固相应力模型.首先,基于局部颗粒浓度和颗粒非均匀分布指数建立了双变量的ANN固相应力模型;进一步将当前网格及其周边网格颗粒浓度组成的序列来详细表征颗粒浓度分布,并建立相应的ANN固相应力模型.然后,将两种模型与EMMS固相应力模型进行了对比并测试了网格分辨率和粗化率对模型的影响.研究表明:基于ANN固相应力模型的模拟结果对EMMS固相应力模型结果有较高的还原度,同时具有一定的网格分辨率无关性和粗化率无关性.     

颗粒物质的分形特征与其物理力学性质的关系探讨

从颗粒流的一般本构方程出发,利用室内大型相对密度仪和直接剪切仪,针对粒径分形分布对孔隙率和内摩擦角的影响进行了实验研究。结果表明,在相同实体密度、弹性恢复系数和几何平均中值粒径条件下,随粒径分维数的增大,孔隙率增大,而内摩擦角、颗粒浓度、剪切应力系数与正应力系数均减小。进而分析了粒径分形特性对颗粒流动行为的影响,得出粒径分维数小的颗粒系统剪切应力大,系统的流动性能差,鲁棒性强;反之,分维数大的颗粒系统剪切应力小,系统的流动性能好,敏感性强。对系统组构的分形特征影响系统颗粒物质崩塌流动行为的机制进行了初步探索,为颗粒流系统的研究提供了新思路。     

Linear and non-linear double diffusive convection in a rotating porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and non-linear stability analyses. The Darcy model that includes the time derivative and Coriolis terms is employed as momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal and solute diffusions that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number and Taylor number on the stability of the system is investigated. The non-linear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out.     

Efficient stochastic FEM for flow in heterogeneous porous media. Part 1: random Gaussian conductivity coefficients

This paper is concerned with the development of efficient iterative methods for solving the linear system of equations arising from stochastic FEMs for single‐phase fluid flow in porous media. It is assumed that the conductivity coefficient varies randomly in space according to some given correlation function and is approximated using a truncated Karhunen–Loève expansion. Distinct discretizations of the deterministic and stochastic spaces are required for implementations of the stochastic FEM. In this paper, the deterministic space is discretized using classical finite elements and the stochastic space using a polynomial chaos expansion. The highly structured linear systems which result from this discretization mean that Krylov subspace iterative solvers are extremely effective. The performance of a range of preconditioned iterative methods is investigated and evaluated in terms of robustness with respect to mesh size and variability of the conductivity coefficient. An efficient symmetric block Gauss–Seidel preconditioner is proposed for problems in which the conductivity coefficient has a large standard deviation.The companion paper, herein, referred to as Part 2, considers the situation in which Gaussian random fields are transformed into lognormal ones by projecting the truncated Karhunen–Loève expansion onto a polynomial chaos basis. This results in a stochastic nonlinear problem because the random fields are represented using polynomial chaos containing terms that are generally nonlinear in the random variables. Copyright © 2013 John Wiley & Sons, Ltd.     

Local Thermal Non-equilibrium and Heterogeneity Effects on the Onset of Convection in an Internally Heated Porous Medium

The effect of local thermal non-equilibrium on the onset of convection in a porous medium consisting of two horizontal layers, each internally heated, is studied analytically. Linear stability theory is applied. Variations of permeability, fluid thermal conductivity, solid thermal conductivity, source strength in the solid and fluid phases, interphase heat-transfer coefficient and porosity are considered. It is found that heterogeneity of permeability, fluid thermal conductivity and source strength in the fluid phase have a major effect; heterogeneity of interphase heat-transfer coefficient and porosity have a lesser effect, while heterogeneity of solid thermal conductivity and source strength in the solid phase are relatively unimportant.     

The flow of closely fitting particles in tapered tubes

The flow of rigid spheres, truncated cones and elastic incompressible spheres in tapered tubes is investigated assuming that the Reynolds equation is valid in the fluid and the linear theory of elasticity is applicable in the solid. It is shown that leading terms in the asymptotic expansion of pressure drop in terms of minimum fluid film thickness for neutrally buoyant rigid spheres and truncated cones are of higher order of magnitude compared to the corresponding terms for the flow of these particles in circular cylindrical tubes. The effect of taper angle on pressure drop is reduced in the case of soft elastic particles because of particle deformations and significant velocities at the particle surface.     

Experimental measurement of thermophysical properties of oxide-water nano-fluids down to ice-point

This paper presents the measurement of the thermal conductivity and the dynamic viscosity of Al₂O₃-water (1-4% particle volume fraction) and TiO₂-water (1-6% particle volume fraction) nano-fluids carried out at atmospheric pressure in the temperature range from 1 to 40 °C, which is particularly interesting for the application of nano-fluids as thermal medium in refrigeration and air-conditioning.The thermal conductivity measurement was performed by using a Transient Hot Disk TPS 2500S apparatus instrumented with a 7577 probe (2.001 mm in radius) having a maximum uncertainty (k = 2) lower than ±5.0% of the reading. The dynamic viscosity measurement and the rheological analysis were carried out by a rotating disc type rheometer Haake Mars II instrumented with a single cone probe (60 mm in diameter and 1° angle) having a maximum uncertainty (k = 2) lower than ±5.0% of the reading.The thermal conductivity measurements of the tested nano-fluids show a great sensitivity to particle volume fraction and temperature and a weak sensitivity to cluster average size: TiO₂-water and Al₂O₃-water nano-fluids show a thermal conductivity enhancement (with reference to pure water) from −2 to 16% and from −2 to 23% respectively.TiO₂-water and Al₂O₃-water nano-fluids exhibit a Newtonian behaviour in all the investigated ranges of temperature and nano-particle volume fraction. The relative viscosity shows a great sensitivity to particle volume fraction and cluster average size and no sensitivity to temperature: TiO₂-water and Al₂O₃-water nano-fluids show a dynamic viscosity increase with respect to pure water from 17 to 210% and from 15 to 150% respectively.Al₂O₃-water nano-fluid seems to be more promising as thermal medium than TiO₂-water nano-fluid, particularly at low thermal level (between ambient temperature and ice point) where TiO₂-water is not suitable showing worse performance than pure water.Present experimental measurements were compared both with available measurements carried out by different researchers and computational models for thermophysical properties of suspensions.     

The distribution of solid particles suspended in a turbulent flow: A stochastic approach

Basic fluid mechanics and stochastic theories are applied to show that the concentration distribution of suspended solid particles in a direction normal to the mean streamlines of a two-dimensional turbulent flow is greatly influenced by the lift force exerted on them in the vicinity of the wall. Analytic solution shows that, when the direction of the mean flow is horizontal, the probability density functionp (y, t) for random displacements of the particles will have a maximum value at a point from the wall where the perpendicular component of the lift force precisely balances particle gravity. Interpretation of experimental observations is presented using this theory.