SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θmodel of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θmodel in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θmodel has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θmodel was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.     

A K-ε TWO-EQUATION TURBULENCE MODEL FOR THE SOLID-LIQUID TWO-PHASE FLOWS

A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and ε equations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.     

A CLOSED SYSTEM OF EQUATIONS FOR DENSE TWO-PHASE FLOW AND EXPRESSIONS OF SHEARING STRESS OF DISPERSED PHA’E AT A WALL

In this paper.each of the two phases in dense two-phase flow is considered ascontinuous medium and the fimdamental equations for two-phase flow are described inEulerian form.The generalized constitutive relation of the Bingham fluid is applied to thedispersed phase with the analysis of physical mechanism of of dense two-phase flow.Theshearing stress of dispersed phase at a wall is used to give a boundary condition.Then amathematical model for dense two-phase flow is obtained.In addition.the expressions ofshearing stress of dispersed phase at a wall is derived according to the fundamental,model ofthe frictional collision between dispersed-phase particles and the wall.     

Similarity solutions of vertical plane wall plume based on finite analytic method

The turbulent flow of vertical plane wall plume with concentration variation was studied with the finite analytical method. The k-epsilon model with the effect of buoyancy on turbulent kinetic energy and its dissipation rate was adopted. There were similarity solutions in the uniform environment for the system of equations including the equation of continuity, the equation of momentum along the flow direction and concentration, and equations of k, epsilon. The finite analytic method was applied to obtain the similarity solution. The calculated data of velocity, relative density difference, the kinetic energy of turbulence and its dissipation rate distribution for vertical plane plumes are in good agreement with the experimental data at the turbulent Schmidt number equal to 1.0. The variations of their maximum value along the direction of main flow were also given. It shows that the present model is good, i.e., the effect of buoyancy on turbulent kinetic energy and its dissipation rate should be taken into account, and the finite analytic method is effective.     

Improved subgrid scale model for dense turbulent solid-liquid two-phase flows

The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules. Assuming that the solid-phase velocity distributions obey the Maxwell equations, the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow, the solid-particle‘s governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations. Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale (SGS) model, in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor, is proposed to model the two-phase governing equations by applying dimension analyses. Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls, the velocity and pressure fields, and the volumetric concentration are calculated. The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.     

A note on a new theory of turbulence

In this paper,the creation and annihilaton of turbulent eddies are described aselementary particles in the quantum field theory.An elementary particle may be consideredas a solid entity as it exists in quantum theory,but a turbulent eddy is often changed in sizeand shape with time due to its energy dissipation in a turbulent field.Therefore,in order toapply the method of the quantum field theory to the turbulent field by analogy,the entity ofthe same eddy should be defined firstly.According to the linearized theory,the turbulenteddies with the similarity character in lime duration may be considered as the entity of thesame eddy,and the creation and annihilation of turbulent eddies without the similarcharacters are related to the interaction term φ_2 in equation(2.6).Then,the creationoperator and annihilation operator similar to those in the quantum field theory are used todescribe the state of turbulent eddy field.Finally,a“Schrǒdinger”equation of turbulenteddies is formulated based upon the nonline     

Mixture flow of particles and power-law fluid in round peristaltic tube

The erythrocyte and blood flowing in the blood vessel can be treated as the two-phase flow of the mixture of particles and a power-law fluid in a peristaltic tube.In the present work, the peristaltic transport of a power-law fluid and the suspension of particles in a tube is investigated by a perturbation method using the long wavelength approximation. The influence of different parameters on the velocity profile and streamlines is explored. Results show that there is a deflection of the flow field when the power-law index n = 0.5 or 1.5 compared with the Newtonian fluid where the trapping zone is symmetric to a certain cross section. The flux rate and reflux of the material are identified,and the conditions under which the reflux appears are determined. Moreover, a reflux phenomenon occurs near the wall. The trapping zone is related to not only the tube geometry and the flow flux but also the fluid properties. Both the length and width of the trapping zone increase with an increase in θ or φ. The trapping zone is more difficult to produce in the shear-thinning fluid than the shear-thickening fluid.     

Transient growth in turbulent particle-laden channel flow

Linear transient growth of optimal perturbations in particle-laden turbulent channel flow is investigated in this work.The problem is formulated in the framework of a Eulerian-Eulerian approach,employing two-way coupling between fine particles and fluid flow.The model is first validated in laminar cases,after which the transient growth of coherent perturbations in turbulent channel flow is investigated,where the mean particle concentration distribution is obtained by direct numerical simulation.It is shown that the optimal small-scale structures for particles are streamwise streaks just below the optimal streamwise velocity streaks,as was previously found in numerical simulations of particle-laden channel flow.This indicates that the optimal growth of perturbations is a dominant mechanism for the distribution of particles in the near-wall region.The current study also considers the transient growth of small-and large-scale perturbations at relatively high Reynolds numbers,which reveals that the optimal large-scale structures for particles are in the near-wall region while the optimal large-scale structures for fluid enter the outer region.     

Unsteady heat and mass transfer in MHD flow over an oscillatory stretching surface with Soret and Dufour effects

In this paper, we study the unsteady coupled heat and mass transfer of two-dimensional MHD fluid over a moving oscillatory stretching surface with Soret and Dufour effects. Viscous dissipation effects are adopted in the energy equation. A uniform magnetic field is applied vertically to the flow direction. The governing equations are reduced to non-linear coupled partial differential equations and solved by means of homotopy analysis method （HAM）. The effects of some physical parameters such as magnetic parameter, Dufour number, Soret number, the Prandtl num- ber and the ratio of the oscillation frequency of the sheet to its stretching rate on the flow and heat transfer characteristics are illustrated and analyzed.     

Modelling of pressure–strain correlation in compressible turbulent flow

Previous studies carried out in the early 1990s conjectured that the main compressible effects could be associated with the dilatational effects of velocity fluctuation. Later, it was shown that the main compressibility effect came from the reduced pressure-strain term due to reduced pressure fluctuations. Although better understanding of the compressible turbulence is generally achieved with the increased DNS and experimental research effort, there are still some discrepancies among these recent findings. Analysis of the DNS and experimental data suggests that some of the discrepancies are apparent if the compressible effect is related to the turbulent Mach number, Mt. From the comparison of two classes of compressible flow, homogenous shear flow and inhomogeneous shear flow （mixing layer）, we found that the effect of compressibility on both classes of shear flow can be characterized in three categories corresponding to three regions of turbulent Mach numbers： the low-Mr, the moderate-Mr and high-Mr regions. In these three regions the effect of compressibility on the growth rate of the turbulent mixing layer thickness is rather different. A simple approach to the reduced pressure-strain effect may not necessarily reduce the mixing-layer growth rate, and may even cause an increase in the growth rate. The present work develops a new second-moment model for the compressible turbulence through the introduction of some blending functions of Mt to account for the compressibility effects on the flow. The model has been successfully applied to the compressible mixing layers.     

A new method for probing the turbulence scalar spectrum by ultrasonic sonar scanning

 A new method is proposed to obtain a turbulent scalar spectrum and the energy dissipation rate in turbulent flow from ultrasonic frequency scanning data. A scanning sonar with frequency varying from 0.5 MHz to 5 MHz has been used to directly probe the energy dissipation rate ɛ and the three-dimensional scalar spectrum E _θ(k). Experiments were conducted in a laboratory open-channel flow in clear water with Reynolds numbers varying from 1.2×10⁵ to 6.9×10⁵. Good agreement is found between measured spectra and those predicted by the Batchelor theory. The energy dissipation rates compare favourably with those obtained from acoustic Doppler velocimeter measurements. Received: 20 March 1997/Accepted: 27 September 1997     

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate ɛ_τ) equations were considered. The emphasis of this paper is focused on the effects of the ɛ_τ-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate of change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters, given by /ɛ, _τ/ɛ_τ, , Sk/ɛ and G/ɛ, becoming constant. Here, and _τ are the production of turbulent kinetic energy k and temperature variance , respectively, ɛ and ɛ_τ are their respective dissipation rates, R is the mixed time scale ratio, G is the buoyant production of k and S is the mean shear gradient. Calculations show that the ɛ_τ-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular ɛ_τ-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the ɛ-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence. Received 21 April 2000 and accepted 21 February 2001     

Turbulence kinetic energy budget in bubbly flows in a vertical duct

Understanding turbulence kinetic energy (TKE) budget in gas–liquid two-phase bubbly flows is indispensable to develop and improve turbulence models for the bubbly flows. In this study, a molecular tagging velocimetry based on photobleaching reaction was applied to turbulent bubbly flows with sub-millimeter bubbles in a vertical square duct to examine the applicability of the k–ε models to the bubbly flows. Effects of bubbles on TKE budget are discussed and a priori tests of the standard and low Reynolds number k–ε models are carried out to examine the applicability of these models to the bubbly flows. The conclusions obtained are as follows: (1) The photobleaching molecular tagging velocimetry is of use for validating turbulence models. (2) The bubbles increase the liquid velocity gradient in the near wall region, and therefore, enhance the production and dissipation rates of TKE. (3) The k–ε models can reasonably evaluate the production rate of TKE in the bubbly flows. (4) The modulations of diffusion due to the bubbles have different characteristics from the diffusion enhancement due to shear-induced turbulence. Hence, the k–ε models fail in evaluating the diffusion rate in the near wall region in the bubbly flows. (5) The k–ε models represent the trends of the production, dissipation, and diffusion rates of ε in the bubbly flow, although more accurate experimental data are required for quantitative validation of the ε equation.     

Improvement of cold-wire response for measurement of temperature dissipation

This work aimed at improving fine-scale measurements using cold-wire anemometry. The dissipation ɛ _θ of the temperature variance was measured on the axis of a heated turbulent round jet. The measurements were performed with a constant current anemometer (CCA) operating fine Pt–10%Rh wires at very low overheat. The CCA developed for this purpose allowed the use of the current injection method in order to estimate the time constant of the wire. In the first part of the paper, it is shown that the time constants obtained for two wire diameters −d=1.2 and d=0.58 μm – compare well with those measured at the same time using two other methods (laser excitation and pulsed wire). Moreover, for these two wires, the estimated time constants were in good agreement with those obtained from a semi-empirical relation. In the second part of the paper, a compensation procedure – post-processing filtering – was developed in order to improved the frequency response of the cold-wire probes. The measurements carried out on the axis of the jet (Re _D =16 500, Re _λ  ≃ 167) showed that the frequency response of the 1.2 μm wire was significantly improved. In fact, the spectral characteristics of the compensated signal obtained with the 1.2 μm wire compared fairly well with those from the 0.58 μm wire. Moreover, the results indicated that the compensation procedure must be applied when the cut-off frequency of the cold-wire f _c is lower than two times the Kolmogorov frequency f _K. In the case where f _c ≃ 0.6f _K, the compensation procedure can reduce the error in the estimate of ɛ _θ by more than 20%. When f _c ≃ 2f _K, the effect of the compensation is reduced to about 5%. Received: 3 November 2000/Accepted: 23 March 2001     

On the Shear Number Effect in Stratified Shear Flow

The influence of the shear number on the turbulence evolution in a stably stratified fluid is investigated using direct numerical simulations on grids with up to 512 × 256 × 256 points. The shear number SK/ε is the ratio of a turbulence time scale K/ε to the shear time scale 1/S. Simulations are performed at two initial values of the Reynolds number Re _Λ= 44.72 and Re _Λ= 89.44. When the shear number is increased from small to moderate values, the nondimensional growth rate γ= (1/SK)dK/dt of the turbulent kinetic energy K increases since the shear forcing and its associated turbulence production is larger. However, a further increase of the shear number from moderate to large values results in a reduction of the growth rate γ and the turbulent kinetic energy K shows long-time decay for sufficiently large values of the shear number. The inhibition of turbulence growth at large shear numbers occurs for both initial values of the Reynolds number and can be explained with the predominance of linear effects over nonlinear effects when the shear number is sufficiently high. It is found that, at the higher initial value of the Reynolds number, the reduction of the growth rate occurs at a higher value of the shear number. The shear number is found to affect spectral space dynamics. Turbulent transport coefficients decrease with increasing shear number. Received 23 June 1998 and accepted 25 February 1999     

On the correlation between enstrophy and energy dissipation rate in a turbulent wake

All three components of the vorticity fluctuation have been measured simultaneously in a turbulent wake using a new eight-sensor vorticity probe. The vorticity fluctuation spectra agree reasonably well with those from a direct numerical simulation of a turbulent channel flow at high wavenumbers. The similarity between the instantaneous energy dissipation rate ε and the instantaneous enstrophy ω² is examined using spectra and probability density functions. The correlation between ω² and ε is evaluated in some detail. The homogeneous value of ε is strongly correlated with ω². The full value of ε and, more especially its isotropic value, are less well correlated with the enstrophy. Conditional averaging indicates that high enstrophy regions are associated with high energy dissipation rate regions.     

On a Generalized Nonlinear K–ε Model and the Use of Extended Thermodynamics in Turbulence

A resent extension of the nonlinear K–ε model is critically discussed from a basic theoretical standpoint. While it was said in the paper that this model was formulated to incorporate relaxation effects, it will be shown that the model is incapable of describing one of the most basic such turbulent flows as is obvious but is described for clarity. It will be shown in detail that this generalized nonlinear K–ε model yields erroneous results for the Reynolds stress tensor when the mean strains are set to zero in a turbulent flow – the return-to-isotropy problem which is one of the most elementary relaxational turbulent flows. It is clear that K–ε type models cannot describe relaxation effects. While their general formalism can describe relaxation effects, the nonlinear K–ε model – which the paper is centered on – cannot. The deviatoric part of the Reynolds stress tensor is predicted to be zero when it actually only gradually relaxes to zero. Since this model was formulated by using the extended thermodynamics, it too will be critically assessed. It will be argued that there is an unsubstantial physical basis for the use of extended thermodynamics in turbulence. The role of Material Frame-Indifference and the implications for future research in turbulence modeling are also discussed. Received 19 February 1998 and accepted 23 October 1998     

On near-wall hot-wire measurements

A specially constructed hot-wire probe was used to obtain very near-wall velocity measurements in both a fully developed turbulent channel flow and flat plate boundary layer flow. The near-wall hot-wire probe, having been calibrated in a specially constructed laminar flow calibration rig, was used to measure the mean streamwise velocity profile, distributions of streamwise and spanwise intensities of turbulence and turbulence kinetic energy k in the viscous sublayer and beyond; these distributions compare very favorably with available DNS results obtained for channel flow. While low Reynolds number effects were clearly evident for the channel flow, these effects are much less distinct for the boundary layer flow. By assuming the dissipating range of eddy sizes to be statistically isotropic and the validity of Taylor's hypothesis, the dissipation rate ɛ _iso in the very near-wall viscous sublayer region and beyond was determined for both the channel and boundary layer flows. It was found that if the convective velocity U _c in Taylor's hypothesis was assumed to be equal to the mean velocity Uˉ at the point of measurement, the value of (ɛ⁺ _iso)₁ thus obtained agrees well with that of (ɛ ⁺)_DNS for y ⁺ ≥ 80 for channel flow; this suggests the validity of assuming U _c=Uˉ and local isotropy for large values of y ⁺. However, if U _c was assumed to be 10.6u _τ , the value of (ɛ⁺ _iso)₂ thus obtained was found to compare reasonably well with the distribution of (ɛ⁺ _iso)_DNS for y ⁺≤ 15. Received: 31 May 1999/Accepted: 20 December 1999