NUMERICAL SIMULATION OF 3-D DENSE SOLID-LIQUID TWO-PHASE TURBULENT FLOW IN A NON-CLOGGING MUD PUMP

A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to TECPLOT as the post-processor, the simulation results isvisualized. The results show the main flow characteristics: There exists backflow and aberrantvelocities at inlet area and a relative velocity slip between two phases; A jet-wake flow pattern isdiscerned around the shroud-suction side area; The relative velocity vector of solid phase iscloser to the pressure surface than that of liquid phase and the trend is more obvious with theincrease of diameter; The kinetic energy of turbulence k and the dissipation rate e reach theirpeaks at the corner of pressure and suction surface. The simulation results show a good agreementwith the experimental flow features in the impeller channel, which prove the turbulent model used isvalid and provide a theoretical design basis to non-clogging pumps.

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.

Deduction and Validation of an Eulerian-Eulerian Model for Turbulent Dilute Two-Phase Flows by Means of the Phase Indicator Function-Disperse Elements Probability Density Function

A statistical formalism overcoming some conceptual and practical difficulties arising in existing two-phase flow (2PHF) mathematical modelling has been applied to propose a model for dilute 2PHF turbulent Hows. Phase interaction terms with a clear physical meaning enter the equations and the formalism provides some guidelines for the avoidance of closure assumptions or the rational approximation of these terms. Continuous phase averaged continuity, momentum, turbulent kinetic energy and turbulence dissipation rate equations have been rigorously and systematically obtained in a single step. These equations display a structure similar to that for single-phase flows. It is also assumed that dispersed phase dynamics is well described by a probability density function (pdf) equation and Eulerian continuity, momentum and fluctuating kinetic energy equations for the dispersed phase are deduced. An extension of the standard k-e turbulence model for the continuous phase is used. A gradient transport model is adopted for the dispersed phase fluctuating fluxes of momentum and kinetic energy at the non-colliding, large inertia limit. This model is then used to predict the behaviour of three axisymmetric turbulent jets of air laden with solid particles varying in size and concentration. Qualitative and quantitative numerical predictions compare reasonably well with the three different sets of experimental results, studying the influence of particle size, loading ratio and flow confinement velocity.

A USM-Θ two-phase turbulence model for simulating dense gas-particle flows

A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.

Effects of the Reynolds number on the mean skin friction decomposition in turbulent channel flows

As the Reynolds number increases, the skin friction has been identified as the dominant drag in many practical applications. In the present paper, the effects of the Reynolds number on the mean skin friction decomposition in turbulent channel flows up to Reτ= 5 200 are investigated based on two different methods, i.e., the FukagataIwamoto-Kasagi(FIK) identity(FUKAGATA, K., IWAMOTO, K., and KASAGI, N.Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows.Physics of Fluids, 14(11), L73–L76(2002)) and the Renard-Deck(RD) identity(DECK,S., RENARD, N., LARAUFIE, R., and WEISS, P.′E. Large-scale contribution to mean wall shear stress in high-Reynolds-number flat-plate boundary layers up to Reθ= 13 650.Journal of Fluid Mechanics, 743, 202–248(2014)). The direct numerical simulation(DNS) data provided by Lee and Moser(LEE, M. and MOSER, R. D. Direct numerical simulation of turbulent channel flow up to Reτ≈ 5 200. Journal of Fluid Mechanics,774, 395–415(2015)) are used. For these two skin friction decomposition methods, their decomposed constituents are discussed and compared for different Reynolds numbers.The integrands of the decomposed constituents are locally analyzed across the boundary layer to assess the actions associated with the inhomogeneity and multi-scale nature of turbulent motion. The scaling of the decomposed constituents and their integrands are presented. In addition, the boundary layer is divided into three sub-regions to evaluate the contributive proportion of each sub-region with an increase in the Reynolds number.

Drag reduction of turbulent channel flows over an anisotropic porous wall with reduced spanwise permeability

The direct numerical simulation (DNS) is carried out for the incompressible viscous turbulent flows over an anisotropic porous wall. Effects of the anisotropic porous wall on turbulence modifications as well as on the turbulent drag reduction are investigated. The simulation is carried out at a friction Reynolds number of 180, which is based on the averaged friction velocity at the interface between the porous medium and the clear fluid domain. The depth of the porous layer ranges from 0.9 to 54 viscous units. The permeability in the spanwise direction is set to be lower than the other directions in the present simulation. The maximum drag reduction obtained is about 15.3% which occurs for a depth of 9 viscous units. The increasing of drag is addressed when the depth of the porous layer is more than 25 wall units. The thinner porous layer restricts the spanwise extension of the streamwise vortices which suppresses the bursting events near the wall. However, for the thicker porous layer, the wall-normal fluctuations are enhanced due to the weakening of the wall-blocking effect which can trigger strong turbulent structures near the wall.

Numerical study of the turbulent channel flow under space-dependent electromagnetic force control at different Reynolds numbers

In this paper, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with the direct numerical simulation(DNS) methods for different Reynolds numbers. A formulation is derived to express the relation between the drag and the Reynolds shear stress. With the application of optimal electromagnetic force, the in-depth relations among characteristic structures in the flow field, mean Reynolds shear stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The regular quasi-streamwise vortex structures, which appear in the flow field, have the same period with that of the electromagnetic force.These structures suppress the random velocity fluctuations, which leads to the absolute value of mean Reynolds shear stress decreasing and the distribution of that moving away from the wall. Moreover, the wave number of optimal electromagnetic force increases,and the scale of the regular quasi-streamwise vortex structures decreases as the Reynolds number increases. Therefore, the rate of drag reduction decreases with the increase in the Reynolds number since the scale of the regular quasi-streamwise vortex structures decreases.

EXPERIMENTAL STUDY ON TURBULENT FEATURES IN THE NEGATIVE TRANSPORT REGION OF ASYMMETRIC PLANE CHANNEL FLOW

Turbulent features of streamwise and vertical components of velocity in the negative transport region of asymmetric plane channel flow have been studied experimentally in details. Experiments show that turbulent fluctuations in negative transport region are suppressed, and their probability distributions are far from Gaussian. Besides, the skewness factors attain their negative maxima at the position of the maximum mean velocity, whereas the flatness factors attain their positive maxima at the same position.

The nature of near-wall convection velocity in turbulent channel flow

A novel notion of turbulent structure the local cascade structure-is introduced to study the convection phenomenon in a turbulent channel flow. A space-time cross-correlation method is used to calculate the convection velocity. It is found that there are two characteristic convection speeds near the wall, one associated with small-scale streaks of a lower speed and another with streamwise vortices and hairpin vortices of a higher speed. The new concept of turbulent structure is powerful to illustrate the dominant role of coherent structures in the near-wall convection, and to reveal also the nature of the convection-the propagation of patterns of velocity fluctuations-which is scale-dependent.

Preferential orientation of tracer spheroids in turbulent channel flow

Axis-symmetric spheroids, such as rod-like and disk-like particles, have been found to orient preferentially in near-wall turbulence by both experiment and numerical simulation. In current work we examined the orientation of inertialess spheroids in a turbulent channel flow at medium friction Reynolds number Reτ=100 given based on the half of channel height. Both elongated prolate spheroid and flat oblate spheroid are considered and further compared with the reference case of spherical particle. The statistical results show that in near wall region the prolate spheroids tend to align in the streamwise direction while the oblate spheroids prefer to orient in the wallnormal direction, which are consistent with earlier observation in low Reynolds number (Reτ=180)wall turbulence. Around the channel center we found that the orientation of spheroids is not fully isotropic, even though the fluid vorticity are almost isotropic. The mechanism that gives rise to such particle orientations in wall-turbulence has been found to be related to fluid Lagrangian stretching and compression (Zhao and Andersson 2016). Therefore, we computed the left Cauchy-Green strain tensor along Lagrangian trajectories of tracer spheroids in current flow field and analyzed the fluid Lagrangian stretching and compression. The results indicated that, similar to the earlier observations, the directions of the Lagrangian stretching and compression in near-wall region are in the streamwise and wall-normal directions, respectively. Furthermore, cross over the channel the prolate spheroids aligned with the direction of Lagrangian stretching but oblate spheroids oriented with the direction of Lagrangian compression. The weak anisotropy of orientations of fluid Lagrangian stretching and compression observed at the channel center could be the reason for the aforementioned modest anisotropic orientation of spheroids in channel central region.

Numerical Research on the Fiber Suspensions in a Turbulent T-shaped Branching Channel Flow

The concentration and orientation of fiber in a turbulent T-shaped branching channel flow are investi-gated numerically. The Reynolds averaged Navier-Stokes equations together with the Reynolds stress turbulent model are solved for the mean flow field and the turbulent kinetic energy. The fluctuating velocities of the fluid are assumed as a random variable with Gaussian distribution whose variance is related to the turbulent kinetic energy. The slender-body theory is used to simulate the fiber motion based on the known mean and fluctuating velocities of the fluid. The results show that at low Reynolds number, fiber concentration is high in the flow separation regions, and fiber orientation throughout the channel is widely distributed with a slight preference of aligning along the horizontal axis. With increasing of Re, the high concentration region disappears, and fiber orientation becomes ho-mogeneous without any preferred direction. At high Reynolds number, fiber concentration increases gradually along the flow direction. The differences in the distribution of concentration and orientation between different fiber aspect ratio are evident only at low Re. Both Re and fiber aspect ratio have small effect on the variance of orientation angle.

Mechanism of controlling turbulent channel flow with the effect of spanwise Lorentz force distribution

A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwise Lorentz force is selected, of which the control effects, such as flow characters, mean Reynolds stress, and drag reductions, at different parameters of amplitude A and wave number k_x are discussed. The results indicate that the control effects vary with the parameter A and k_x. With the increase of A, the drag reduction rate D_r first increases and then decreases rapidly at low k_x,and slowly at high k_x. The low drag reduction（or even drag increase） is due to a weak suppression or even the enhancements of the random velocity fluctuation and mean Reynolds stress. The efficient drag reduction is due to the quasi-streamwise vortex structure induced by Lorentz force, which contributes to suppressing the random velocity fluctuation and mean Reynolds stress, and the negative vorticity improves the distribution of streamwise velocity. Therefore, the optimal control effect with a drag reduction of up to 58% can be obtained.

Comparison of a Full Second-Order Moment Model and an Algebraic Stress Two-Phase Turbulence Model for Simulating Bubble-Liquid Flows in a Bubble Column

A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.

Numerical simulation of the effect of void fraction and inlet velocity on two-phase turbulence in bubble-liquid flows

There are contradicted opinions on whether bubbles enhance or reduce the liquid turbulence. In this paper, the effect of void fraction and inlet velocity on the bubble-liquid two-phase turbulence of the multiple bubble-liquid jets in a two-dimensional channel is studied by using the two-phase second-order moment turbulence model. The results confirm the phenomena observed in experiments and reported in references that at a low void fraction and low inlet velocities the bubbles enhance the liquid turbulence, whereas at a high void fraction and high inlet velocities the bubbles reduce the liquid turbulence.

Analysis of dynamic of two-phase flow in small channel based on phase space reconstruction combined with data reduction sub-frequency band wavelet

A new method of nonlinear analysis is established by combining phase space reconstruction and data reduction sub-frequency band wavelet. This method is applied to two types of chaotic dynamic systems(Lorenz and Rssler) to examine the anti-noise ability for complex systems. Results show that the nonlinear dynamic system analysis method resists noise and reveals the internal dynamics of a weak signal from noise pollution. On this basis, the vertical upward gas–liquid two-phase flow in a 2 mm × 0.81 mm small rectangular channel is investigated. The frequency and energy distributions of the main oscillation mode are revealed by analyzing the time–frequency spectra of the pressure signals of different flow patterns. The positive power spectral density of singular-value frequency entropy and the damping ratio are extracted to characterize the evolution of flow patterns and achieve accurate recognition of different vertical upward gas–liquid flow patterns(bubbly flow:100%, slug flow: 92%, churn flow: 96%, annular flow: 100%). The proposed analysis method will enrich the dynamics theory of multi-phase flow in small channel.

Multi-scale analysis of subgrid stress and energy dissipation in turbulent channel flow

In present study, the subgrid scale （SGS） stress and dissipation for multiscale formulation of large eddy simulation are analyzed using the data of turbulent channel flow at Ret = 180 obtained by direct numerical simulation. It is found that the small scale SGS stress is much smaller than the large scale SGS stress for all the stress components. The dominant contributor to large scale SGS stress is the cross stress between small scale and subgrid scale motions, while the cross stress between large scale and subgrid scale motions make major contributions to small scale SGS stress. The energy transfer from resolved large scales to subgrid scales is mainly caused by SGS Reynolds stress, while that between resolved small scales and subgrid scales are mainly due to the cross stress. The multiscale formulation of SGS models are evaluated a priori, and it is found that the small- small model is superior to other variants in terms of SGS dissipation.

Observed characteristics of flow,water mass,and turbulent mixing in the Preparis Channel

Preparis Channel is the very important exchange path of energy and materials between the northern Bay of Bengal and Andaman Sea(AS).A set of hydrographic measurements,a microstructure profiler,and a deep mooring were used to determine the characteristics of water masses,turbulent mixing,and flows in the Preparis Channel.The unprecedented short-term mooring data reveal that a deep current in the deep narrow passage(below 400 m)of the Preparis Channel flows toward the Bay of Bengal(BoB)with a mean along-stream velocity of 25.26 cm/s at depth of 540 m;above the deep current,there are a relatively weak current flows toward the AS with a mean along-stream velocity of 15.46 cm/s between 500 m and 520 m,and another weak current flows toward the BoB between 430 m and 500 m.Thus,a sandwiched vertical structure of deep currents(below 400 m)is present in the Preparis Channel.The volume transport below 400 m is 0.06 Sv(1 Sv=106 m^(3)/s)from the AS to the BoB.In the upper layer(shallower than 300 m),the sea water of the AS is relatively warmer and fresher than that in the BoB,indicating a strong exchange through the channel.Microstructure profiler observations reveal that the turbulent diffusivity in the upper layer of the Preparis Channel reaches O(10−4 m^(2)/s),one order larger than that in the interior of the BoB and over the continental slope of the northern AS.We speculate that energetic high-mode internal tides in the Preparis Channel contribute to elevated turbulent mixing.In addition,a local“hotspot”of turbidity is identified at the deep mooring site,at depth of about 100 m,which corresponds to the location of elevated turbulent mixing in the Preparis Channel.

On the capability of the curvilinear immersed boundary method in predicting near-wall turbulence of turbulent channel flows

The immersed boundary method has been widely used for simulating flows over complex geometries.However,its accuracy in predicting the statistics of near-wall turbulence has not been fully tested.In this work,we evaluate the capability of the curvilinear immersed boundary(CURVIB)method in predicting near-wall velocity and pressure fluctuations in turbulent channel flows.Simulation results show that quantities including the time-averaged streamwise velocity,the rms(root-mean-square)of velocity fluctuations,the rms of vorticity fluctuations,the shear stresses,and the correlation coefficients of u'and v"computed from the CURVIB simulations are in good agreement with those from the body-fitted simulations.More importantly,it is found that the time-averaged pressure,the rms and wavenumber-frequency spectra of pressure fluctuations computed using the CURVIB method agree well with the body-fitted results.

Numerical Study of Turbulent Periodic Flow and Heat Transfer in a Square Channel with Different Ribs

A numerical investigation has been carried out to examine turbulent flow and heat transfer characteristics in a three-dimensional ribbed square channels. Fluent 6.3 CFD code has been used. The governing equations are discretized by the second order upwind differencing scheme, decoupling with the SIMPLE (semi-implicit method for pressure linked equations) algorithm and are solved using a finite volume approach. The fluid flow and heat transfer characteristics are presented for the Reynolds numbers based on the channel hydraulic diameter ranging from 104 to 4 ′ 104. The effects of rib shape and orientation on heat transfer and pressure drop in the channel are investigated for six different rib configurations. Rib arrays of 45° inclined and 45° V-shaped are mounted in inline and staggered arrangements on the lower and upper walls of the channel. In addition, the performance of these ribs is also compared with the 90° transverse ribs.

Error of large-eddy simulation in the wall pressure fluctuation of a turbulent channel flow

We analyze the error of large-eddy simulation(LES)in wall pressure fluctuation of a turbulent channel flow.To separate different sources of the error,we conduct both direct numerical simulations(DNS)and LES,and apply an explicit filter on DNS data to obtain filtered DNS(FDNS)data.The error of LES is consequently decomposed into two parts:The first part is the error of FDNS with respect to DNS,which quantifies the influence of the filter operation.The second part is the difference between LES and FDNS induced by the error of LES in velocity field.By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation,it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation.The present study provides a basis on future LES studies of the wall pressure fluctuation.