FLOW OF A TRAIN OF DEFORMABLE FLUID PARTICLES IN A TUBE

I.IntroductionWangandSkalakll]studiedforsolidspheresintubes.HymanandSkalak12--31studiedthemovementofacompactfluidspheres(non-deformable)intubes.Theyonlystudiedcompactdropswhoseradiusesarenotbiggerthan0.7.Pozrikidisl4]studiedthenlotiollofdeformabletluidparticlesintubes.Inthispaper,westudiedtwo-phasenowregime,consistingofaperiodictrainofequallyspacedliquidslugsmovingaxisymmetricallyinacontinuousliquidphasewhichisimmisciblewiththeslugphase(Fig.l(a)).Wereportnumericalsolutionsforthehydrodynamic…     

Extension of some numerical schemes to the analysis of gas and particle mixtures

In this paper, several numerical schemes are extended to obtain approximate solutions to the system of equations encountered in the analysis of multiphase mixtures of gas and particles. Both dense and dilute mixtures are studied, the gas is modelled as a perfect gas and the solid is considered incompressible. Although the tests employed throughout this work for studying the behaviour of the schemes are essentially one dimensional, the finite volume method developed permits its application to multidimensional problems in unstructured grids. Copyright © 2007 John Wiley & Sons, Ltd.     

两相流中圆柱状粒子阻力的研究

将圆柱状粒子和圆球粒子建立等效关系，由等效面积法得到了作用在圆柱状粒子上的阻力，然后通过有限体积法对圆柱状粒子的阻力进行了数值求解，经比较发现，在层流状态时由等效面积法得到的圆柱状粒子的阻力要乘上一个1．5—2的常系数．     

明渠气二相流的双流体模型

从解决平均二相流基本方程的封闭问题出发,分析了水利工程中气水二相流的特点,介绍了据此提出的可用于水利工程稀疏气泡流计算的双流体模型,并对明渠气水二相流进行了算例验证。     

Experimental and numerical investigations of pressure drop in a rectangular duct with modified power law fluids

Numerical solutions are presented for fully developed laminar flow for a modified power law fluid (MPL) in a rectangular duct. The solutions are applicable to pseudoplastic fluids over a wide shear rate range from Newtonian behavior at low shear rates, through a transition region, to power law behavior at higher shear rates. The analysis identified a dimensionless shear rate parameter which, for a given set of operating conditions, specifies where in the shear rate range a particular system is operating, i.e. in the Newtonian, transition, or power law regions. The numerical results of the friction factor times Reynolds number for the Newtonian and power law region are compared with previously published results showing agreement within 0.05% in the Newtonian region, and 0.9% and 5.1% in the power law region. Rheological flow curves were measured for three CMC-7H4 solutions and were found to be well represented by the MPL constitutive equation. The friction factor times Reynolds number values were measured in the transition region for which previous measurements were unavailable. Good agreement was found between experiment and calculation thus confirming the validity of the analysis.     

Extension of a high‐resolution scheme to 1D liquid–gas flow

This paper presents the extension of a high‐resolution conservative scheme to the one‐dimensional one‐pressure six‐equation two‐fluid flow model. Only mixtures of water and air have been considered in this study, both fluids have been characterized using simple equations of state, namely stiffened gas for the liquid phase and perfect gas for the gas phase. The resulting scheme is explicit and first‐order accurate in space and time. A second‐order version of the scheme has also been derived using the MUSCL strategy and slope limiters. Some numerical results show the good capabilities of this type of schemes in the solution of discontinuities in two‐fluid flow problems, all of them are based on water/air numerical benchmarks widely used in the two‐phase flow literature. Copyright © 2005 John Wiley & Sons, Ltd.     

The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

Direct numerical simulation of confined swirling jets

This study investigates the Lagrangian acceleration and velocity of fluid particles in swirling flows via direct numerical simulation. The intermittency characteristics of acceleration and velocity of fluid particles are investigated at different swirl numbers and Reynolds numbers. The flatness factor and trajectory curvature are used to analyse the effect of Lagrangian intermittency. The joint probability density function of Lagrangian acceleration and turbulence intensity is shown to explain the augmentation effect of Lagrangian intermittency by the strongly swirling levels under the relatively low intensity of turbulence. In addition, the correlation between the Lagrangian acceleration and the turbulence intensity is enhanced as the swirl level increases. It shows the important effect of swirl on the motion behaviour of fluid particles in the strongly swirling flows.     

Qualitative Questions in Fluid Mechanics

In this paper in honour of Professor Leen van Wijngaarden, some propositions about fluid mechanics are discussed. First, basic fluid mechanics research should be judged as much by its progress in clarifying the essential questions about the phenomena of fluid flow and in establishing general concepts, as by its contribution to the solutions of specific problems. In fact, the latter often contribute to the former. Both aspects attract good students to the subject.Second, researchers make more progress and are likely to impress a journal editor when they relate their problems to general physical and/or mathematical considerations, and when they analyse and present their results in a wide but fluid mechanically relevant context, for example through symmetry considerations, invariants (including dimensionless groups, scaling laws and topological constraints), differential properties (or jumps, wiggles and swirls), and through raising new questions and concepts of general significance from studies of specific flows.Lastly, decisions by organisations and individuals about future research directions also benefit from being considered in a wide conceptual framework.     

Viscosity of a Suspension with a Cubic Array of Spheres in a Shear Flow

Viscous fluid flow past an infinite periodic array of rigid spheres of the same radius is considered. A solution of the Stokes equations periodic in three variables is obtained for viscous incompressible flow with a linear velocity profile. The solution takes into account the hydrodynamic interaction of an infinite number of particles in the array. An expression for the effective viscosity of a suspension with a cubic array of particles is obtained.     

Advection upwinding splitting method to solve a compressible two‐fluid model

In this study, the advection upwinding splitting method (AUSM) is modified for the resolution of two‐phase mixtures with interfaces. The compressible two‐fluid model proposed by Saurel and Abgrall is chosen as the model equations. Dense and dilute phases are described in terms of the volume fraction and equations of state to represent multi‐phase mixtures. Test cases involving an air–water shock tube, water faucet, and dilute particulate turbulent flows through a 90° bend are used to verify the current work. It is shown that the AUSM based on flux differences (AUSMD) contains the mechanism to correctly capture the contact discontinuity and interfaces between phases. In addition, a successful application to dilute particulate turbulence flows by the AUSMD is demonstrated. Copyright © 2001 John Wiley & Sons, Ltd.     

Stability analysis of annular flow structure, using a discrete form of the “two-fluid model”

The differential form of the “two-fluid model” for annular flow, neglecting surface tension, is ill-posed, and it is not suited for examining the stability of the steady-state solutions with respect to the average film thickness. It is shown here that a discrete (difference) representation of the two-fluid model may lead to an appropriate criterion for the stability of the steady-state solutions. Exactly the same criterion is obtained from the requirement that the kinematic waves will propagate in the downstream direction. The suggested discrete form of the “two-fluid model” is used to perform transient simulation and for examining the system response to finite disturbances.     

Analytical and numerical investigations of the reflection of asymmetric nonstationary shock waves

The reflection of asymmetric nonstationary shock waves is analytically and numerically studied in this paper. An analytical approach, which is a combination of the shock dynamic and shock polar methods, is advanced to predict the reflection wave configurations. The numerical simulations are performed by the finite volume method based on the second-order MUSCL-Hancock scheme and the HLLC approximate Riemann solver, with the self-adaptive unstructured mesh. It is found that the transition between the overall regular reflection and overall Mach reflection in the asymmetric nonstationary reflection agrees with the detachment criterion, which is analogous to the reflection in pseudo-steady flows (i.e. shock reflection over a wedge). Some special reflection wave configurations, which have never been observed in steady or nonstationary shock reflections so far, are found to exist in this asymmetric reflection. Furthermore, the domains and boundaries of various overall reflection wave configurations are analytically predicted, and the effect of mis-synchronization is also discussed.     

颗粒存在对网栅后气相流动湍流特性的影响

应用 L DV测试技术对方管内网栅后的气固两相流动 ,在颗粒平均粒径 Dp =0 .2 1 mm,颗粒质量浓度为 0 .2 4%、0 .36 % ;Dp =0 .35 mm,颗粒质量浓度为 0 .1 2 %、0 .2 1 %、0 .335 % ;Dp =0 .6 mm,颗粒质量浓度为 0 .1 6 %、0 .2 45 %、0 .34 5 % ;Dp =0 .9mm,颗粒质量浓度为 0 .2 0 5 %、0 .30 %多种工况下进行了测量 ,得出各种工况下气流脉动速度、湍动能沿流动方向的衰减规律 ,通过与纯气流条件下的实验结果比较 ,分析了颗粒浓度及颗粒尺寸对网栅后气相流动湍流特性的影响。根据测量结果 ,提出了一个在有固相颗粒存在时 ,关于湍流模型方程中模型常数 C2的修定方法。     

群沙射流结构的实验研究

本文以90－150μm的标准细沙为研究对象，利用超声波速度剖面测量仪和数码摄录机进行了垂直入水群沙射流的实验研究，给出了的群沙入水后的运动规律。我们发现，群沙湍射入水的一霎那，底部瞬时产生反射波，由它带动的流速与群沙入水带动的流速同量级。     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

The numerical and experimental simulation of hypervelocity flow around the HYFLEX vehicle forebody

Numerical and experimental techniques are used to model the flow and pressure distribution around the forebody of the HYFLEX hypersonic flight vehicle. We compare numerical simulation results with modified Newtonian theory and flight data to determine the accuracy of the computational fluid dynamics (CFD) technique used. The numerical simulations closely match the trends in flight data, and show that real gas effects have a small but significant influence on the nose pressure distribution. We also present pressure results from a scale-model tested in a shock tunnel, and compare them with simulation results. For the shock tunnel experiment, the model was placed such that part of the upper surface was in a region of the test flow where nonuniformities were significant, and it was shown that the numerical simulation could adequately capture these experimental flow features. The binary scaling parameter (describing the similarity in species dissociation between flight and model) was used to design the scale-model tests in the shock tunnel, and its effectiveness is discussed. We find that matching the flight Mach number in the shock tunnel experiment is not critical for reproducing flight pressure data, so long as flight velocity is matched, and binary scaling is maintained. Received 11 June 1998 / Accepted 1 September 1998     

The impulsive response of sandwich beams: Analytical and numerical investigation of regimes of behaviour

An analytical model is developed to classify the impulsive response of sandwich beams based on the relative time-scales of core compression and the bending/stretching response of the sandwich beam. It is shown that an overlap in time scales leads to a coupled response and to the possibility of an enhanced shock resistance. Four regimes of behaviour are defined: decoupled responses with the sandwich core densifying partially or completely, and coupled responses with partial or full core densification. These regimes are marked on maps with axes chosen from the sandwich beam transverse core strength, the sandwich beam aspect ratio and the level of blast impulse. In addition to predicting the time-scales involved in the response of the sandwich beam, the analytical model is used to estimate the back face deflection, the degree of core compression and the magnitude of the support reactions. The predictions of the analytical model are compared with finite element (FE) simulations of impulsively loaded sandwich beams comprising an anisotropic foam core and elastic, ideally plastic face-sheets. The analytical and numerical predictions are in good agreement up to the end of core compression. However, the analytical model under-predicts the peak back face deflection and over-predicts the support reactions, especially for sandwich beams with high strength cores. The FE calculations are employed to construct design charts to select the optimum transverse core strength that either minimises the back face deflections or support reactions for a given sandwich beam aspect ratio or blast impulse. Typically, the value of the transverse core strength that minimises the back face deflection also minimises the support reactions. However, the optimal core strength depends on the level of blast impulse, with higher strength cores required for greater blasts.     

垂直下降管中两相流的流型以及液相粘度对流型的影响

对空气－油在垂直下降管中的流型进行了实验研究，采用的管径为２９ｍｍ，油和空气的折算流速分别达到４ｍ／ｓ和２０ｍ／ｓ，并借助于压降脉动分析和目测观察相结合的方法来进行流型的识别。研究表明，油气两相流的流型不同于低粘液体的两相流流动，通过实验研究并结合前人的研究成果，给出了液相粘度对流型转变的影响趋势。