Turbulent swirling flow inside a short pipe interacting with a conical bluff body was simulated using the commercial CFD code Fluent.The geometry used is a simplified version of a novel liquid/gas separator used in mu...Turbulent swirling flow inside a short pipe interacting with a conical bluff body was simulated using the commercial CFD code Fluent.The geometry used is a simplified version of a novel liquid/gas separator used in multiphase flow metering.Three turbulence models,belonging to the Reynolds averaged Navier-Stokes(RANS)equations framework,are used.These are,RNG k-ε,SST k-ωand the full Reynolds stress model(RSM)in their steady and unsteady versions.Steady and unsteady RSM simulations show similar behavior.Compared to other turbulence models,they yield the best predictions of the mean velocity profiles though they exhibit some discrepancies in the core region.The influence of the Reynolds number on velocity profiles,swirl decay,and wall pressure on the bluff body are also presented.For Reynolds numbers generating a Rankine-like velocity profile,the width and magnitude of flow reversal zone decreases along the pipe axis disappearing downstream for lower Reynolds numbers.The tangential velocity peaks increase with increasing Reynolds number.The swirl decay rate follows an exponential form in accordance with the existing literature.These flow features would affect the performance of the real separator and,thus,the multiphase flow meter,noticeably.展开更多
The motion of fibers in turbulent pipe flow was simulated by 3-D integral method based on the slender body theory and simplified model of turbulence.The orientation distribution of fibers in the computational area for...The motion of fibers in turbulent pipe flow was simulated by 3-D integral method based on the slender body theory and simplified model of turbulence.The orientation distribution of fibers in the computational area for different Re numbers was computed.The results which were consistent with the experimental ones show that the fluctuation velocity of turbulence cause fibers to orient randomly.The orientation distributions become broader as the Re number increases.Then the fluctuation velocity and angular velocity of fibers were obtained.Both are affected by the fluctuation velocity of turbulence.The fluctuation velocity intensity of fiber is stronger at longitudinal than at lateral,while it was opposite for the fluctuation angular velocity intensity of fibers.Finally,the spatial distribution of fiber was given.It is obvious that the fiber dispersion is strenghened with the increase of Re numbers.展开更多
To predict the characteristics of dense liquid-solid two-phase flow, K-Ε-T model is established, in which the turbulent flow of fluid phase was described with fluid turbulent kinetic energy Kf and its dissipation ra...To predict the characteristics of dense liquid-solid two-phase flow, K-Ε-T model is established, in which the turbulent flow of fluid phase was described with fluid turbulent kinetic energy Kf and its dissipation rate Εf, and the particles random motion was described with particle turbulent energy Kp and its dissipation rate Εp and pseudothermal temperature Tp. The governing equations were also derived. With K-Ε-T model, numerical study of dense liquid-solid two-phase turbulent up-flow in a pipe is performed. The calculated results are in good agreement with experimental data of Alajbegovic et al. (1994), and some flow features are captured.展开更多
In this paper we use the Green function method to solve the problem of steady one dimensional flow of an incompressible viscous, electrically conducting fluid through a pipe with partial circular ring cross sec- tion ...In this paper we use the Green function method to solve the problem of steady one dimensional flow of an incompressible viscous, electrically conducting fluid through a pipe with partial circular ring cross sec- tion and one with annular cross section, in the presence of an applied transverse uniform magnetic field, We ob- tain analytic solutions and carry out some numerical calculations of the velocity distribution and induced magnet- ic field.展开更多
The optimal transient growth process of perturbations driven by the pressure gradient is studied in a turbulent pipe flow. A new computational method is proposed, based on the projection operators which project the go...The optimal transient growth process of perturbations driven by the pressure gradient is studied in a turbulent pipe flow. A new computational method is proposed, based on the projection operators which project the governing equations onto the sub- space spanned by the radial vorticity and radial velocity. The method is validated by comparing with the previous studies. Two peaks of the maximum transient growth am- plification curve are found at different Reynolds numbers ranging from 20 000 to 250 000. The optimal flow structures are obtained and compared with the experiments and DNS results. The location of the outer peak is at the azimuthal wave number n = 1, while the location of the inner peak is varying with the Reynolds number. It is observed that the velocity streaks in the buffer layer with a spacing of 100δv are the most amplified flow structures. Finally, we consider the optimal transient growth time and its dependence on the azimuthal wave length. It shows a self-similar behavior for perturbations of different scales in the optimal transient growth process.展开更多
We address the flow of incompressible fluid with a pressure-dependent viscosity through a pipe with helical shape. The viscosity-pressure relation is defined by the Barus law. The thickness of the pipe and the helix s...We address the flow of incompressible fluid with a pressure-dependent viscosity through a pipe with helical shape. The viscosity-pressure relation is defined by the Barus law. The thickness of the pipe and the helix step are assumed to be of the same order and considered as the small parameter. After transforming the starting problem, we compute the asymptotic solution using curvilinear coordinates and standard perturbation technique. The solution is provided in the explicit form clearly showing the influence of viscosity-pressure dependence and pipe's geometry on the effective flow.展开更多
The nature of turbulent swirling and rotating flow in a straight pipe is investigated using a family of near-wall two-equation models. Specifically, the viability of three different near-wall two-equation models is as...The nature of turbulent swirling and rotating flow in a straight pipe is investigated using a family of near-wall two-equation models. Specifically, the viability of three different near-wall two-equation models is assessed. These models are asymptotically consistent near the wall. The first two models, one with isotropic and another with anisotropic eddy viscosity invoked, solved a dissipation rate equation with no explicit correction made to account for swirl and flow rotation. The third model assumes an isotropic eddy viscosity but solves an improved dissipation rate equation that has explicit corrections made to account for swirl and flow rotation. Calculations of turbulent flows in the swirl number range 0.25 - 1.3 with and without a central recirculation region reveal that, with the exception of the third model, neither one of the other two models can replicate the mean field at the swirl numbers tested. Furthermore, taking stress anisotropy into account also fails to model swirl effect correctly. Significant improvements can be realized from the third model, which is based on an improved dissipation rate equation and the assumption of isotropic eddy viscosity. The predicted mean flow and turbulence statistics correlate well with measurements at low swirl. At high swirl, the two-equation model with an improved dissipation rate equation can still be used to model swirling and rotating pipe flows with a central recirculation region. However, its simulation of flows without a central recirculation region is not as satisfactory.展开更多
The computed orientation distribution of fibers immersed in laminar pipe flows showed that the longitudinal distributions are wide for small Reynolds numbers and become narrower with increasing Re.For low Re number.th...The computed orientation distribution of fibers immersed in laminar pipe flows showed that the longitudinal distributions are wide for small Reynolds numbers and become narrower with increasing Re.For low Re number.the axial orientation distributions are broad with almost no preferred orientations.For high Re numer,the axial distribution becomes narrow.with sharp maxima..The mean values of the longitudinal orientation depend strongly on the Re number.The computed results are in qualitative agreemernt with relevant experimental results.展开更多
A gate valve is one of the main elements of a circular pipeline, but the flow characteristics around the gate valve are hardly known. In this study, clarification of the flow field in front of the gate valve model in ...A gate valve is one of the main elements of a circular pipeline, but the flow characteristics around the gate valve are hardly known. In this study, clarification of the flow field in front of the gate valve model in a pipe flow via flow visualization and PIV analysis was attempted. As a result, four kinds of steady necklace-type vortex systems, 2-vortex, 4-vortex, 6-vortex and 8-vortex systems, were clearly observed in a Reynolds number between 290 and 2130. In addition, the main vortex was observed in the Reynolds number range between 2130 and 4870 with difficulty. On this account, both the center position and vorticity in the main vortex are presented against Reynolds number.展开更多
The aim of this study is to evaluate a three-equation turbulence model applied to pipe flow. Uncertainty is approximated by comparing with published direct numerical simulation results for fully-developed average pipe...The aim of this study is to evaluate a three-equation turbulence model applied to pipe flow. Uncertainty is approximated by comparing with published direct numerical simulation results for fully-developed average pipe flow. The model is based on the Reynolds averaged Navier-Stokes equations. Boussinesq hypothesis is invoked for determining the Reynolds stresses. Three local length scales are solved, based on which the eddy viscosity is calculated. There are two parameters in the model;one accounts for surface roughness and the other is possibly attributed to the fluid. Error in the mean axial velocity and Reynolds stress is found to be negligible.展开更多
The Moody Diagram is widely used to determine the friction factor for fluid flow in pipes. The diagram combines the effects of Reynolds number and relative roughness to determine the friction factor. The relationship ...The Moody Diagram is widely used to determine the friction factor for fluid flow in pipes. The diagram combines the effects of Reynolds number and relative roughness to determine the friction factor. The relationship is highly non-linear and appears to have a complex interaction between viscous and boundary roughness effects. The Moody Diagram is based on predictions from an equation developed by Colebrook in 1939. The relationship requires an iteration process to make predictions. While empirical relationships have been developed that provide good predictions without an iteration process, no one has fully explained the cause for the observed results. The objective of this paper is to present a logical development for prediction of the friction factor. An equation has been developed that models the summed effect of both the laminar sublayer and the boundary roughness on the fluid profile and the resulting friction factor for pipes. The new equation does not require an iteration procedure to obtain values for the friction factor. Predicted results match well with values generated from Colebrook’s work that is expressed in the Moody Diagram. Predictions are within one percent of Colebrook values and generally less than 0.3 percent error from his values. The development provides insight to how processes operating at the boundary cause the friction factor to change.展开更多
In this study, the magnetohydrodynamic (MHD) flow through a circular pipe under the influence of a transverse mag- netic field when the outside medium is also electrically conducting is solved numerically by using FEM...In this study, the magnetohydrodynamic (MHD) flow through a circular pipe under the influence of a transverse mag- netic field when the outside medium is also electrically conducting is solved numerically by using FEM-BEM coupling approach. The coupled partial differential equations defined for the interior medium are transformed into homogenous modified Helmholtz equations. For the exterior medium on an infinite region, the Laplace equation is considered for the exterior magnetic field. Unknowns in the equations are also related with the corresponding Dirichlet and Neumann type coupled boundary conditions. Unknown values of the magnetic field on the boundary and for the exterior region are obtained by using BEM, and the unknown velocity and magnetic field inside the pipe are obtained by using SUPG type stabilized FEM. Computations are carried for very high values of magnetic Reynolds numbers Rm1, Reynolds number Re and magnetic pressure Rh of the fluid. The results show that using stabilized method enables us to get stable and accurate numerical approximations consistent with the physical configuration of the problem over rough mesh which also results a cheap computational cost.展开更多
The steady laminar pipe flow of a suspension with a gas volume fraction ∅≤0.5 and small or intermediate bubble deformations in long, and straight sections of a circular pipe is calculated. The calcula...The steady laminar pipe flow of a suspension with a gas volume fraction ∅≤0.5 and small or intermediate bubble deformations in long, and straight sections of a circular pipe is calculated. The calculations are based on the constitutive equation that was originally derived for dilute emulsions and further developed for concentrated suspensions containing bubbles. In contrast to the literature, an analytical procedure is used to determine the solution of a pipe flow more accurately. The results are presented and discussed with respect to the Reynolds number Re and capillary number Ca. If Ca 1, a bubble suspension has a parabolic velocity profile indicating a Newtonian rheology. If Ca ≈1, two regimes of flow are observed in agreement with the literature;that is, an inner plug flow where deformation rates are low and an outer flow where deformation rates are high. These results imply that, if Ca ∅?and that, if Ca ≥1, the opposite effect occurs;that is, the Reynolds number Re increases with increasing gas volume fraction.展开更多
The calculation of the mechanical energy loss is one of the fundamental problems in the field of Hydraulics and Engineering Fluid Mechanics. However, for a non-uniform flow the relation between the mechanical energy l...The calculation of the mechanical energy loss is one of the fundamental problems in the field of Hydraulics and Engineering Fluid Mechanics. However, for a non-uniform flow the relation between the mechanical energy loss in a volume of fluid and the kinematical and dynamical characteristics of the flow field is not clearly established. In this paper a new mechanical energy equation for the incompressible steady non-uniform pipe flow of homogeneous fluid is derived, which includes the variation of the mean turbulent kinetic energy, and the formula for the calculation of the mechanical energy transformation loss for the non-uniform flow between two cross sections is obtained based on this equation. This formula can be simplified to the Darcy-Weisbach formula for the uniform flow as widely used in Hydraulics. Furthermore, the contributions of the mechanical energy loss relative to the time averaged velocity gradient and the dissipation of the turbulent kinetic energy in the turbulent uniform pipe flow are discussed, and the contributions of the mechanical energy loss in the viscous sublayer, the buffer layer and the region above the buffer layer for the turbulent uniform flow are also analyzed.展开更多
A linear stability analysis on the circular pipe flow of fiber suspensions ispresented. The fiber orientation tensors are used to describe the characteristics of theorientation dislribution of fibers. The constitutive...A linear stability analysis on the circular pipe flow of fiber suspensions ispresented. The fiber orientation tensors are used to describe the characteristics of theorientation dislribution of fibers. The constitutive equation for the fiber suspensions is set upand the modified Orr-Sommerfeld stability equation is derived. An improved finite difference methodwith high order accuracy is employed to solve the e-quation. The Newtonian pipe Poiseuille flowcorresponding to H = 0 is also analyzed for comparison. The results reveal that the fiber additiveswill enhance the flow stability, the degree of enhancement becomes high with the increases of theparameter H which accounts for the fiber resistance to the stretching along its axis. Fibersuspensions with large H can suppress the influence of high Reynolds number on the flow stability.Particularly, fibers give a higher attenuation of the short waves of disturbance.展开更多
The deposition of particles in turbulent pipe flow was investigated in terms of two mechanisms, turbulent and thermophoretic. A general equation incorporating these two mechanisms was formulated to calculate the depos...The deposition of particles in turbulent pipe flow was investigated in terms of two mechanisms, turbulent and thermophoretic. A general equation incorporating these two mechanisms was formulated to calculate the deposition efficiency of aerosol particles in turbulent pipe flow together with thermophoretic deposition. The validity of the equation was confirmed by good agreement between calculated and measured results.展开更多
In the non-spherical particulate turbulent flows, a set of new fluidfluctuating velocity equations with the non-spherical particle source term were derived, then a newmethod, which treats the slowly varying functions ...In the non-spherical particulate turbulent flows, a set of new fluidfluctuating velocity equations with the non-spherical particle source term were derived, then a newmethod, which treats the slowly varying functions and rapidly varying functions separately, wasproposed to solve the equations, and finally the turbulent intensity and the Reynolds stress of theflu-id were obtained by calculating the fluctuating velocity statlsti-cally. The equations andmethod were used to a paniculate tur-bulent pipe flow. The results show that the turbulent intensityand the Reynolds stress are decreased almost inverse proportion-ally to the fluctuating velocityratio of particle to fluid. Non-spherical particles have a greater suppressing effect on thetur-bulence than the spherical particles. The particles with short re-laxation time reduce theturbulence intensity of fluid, while the particles with long relaxation time increase the turbulenceinten-sity of fluid. For fixed particle and fluid, the small particles sup-press the turbulence andthe large particles increase the turbu-ience.展开更多
Fully developed turbulence measurements in pipe flow were made in theReynolds number range from 10 X 10~3 to 350 X 10~3 with hot-wire anemometer and a Pilot tube.Comparisons were made with the experimental results of ...Fully developed turbulence measurements in pipe flow were made in theReynolds number range from 10 X 10~3 to 350 X 10~3 with hot-wire anemometer and a Pilot tube.Comparisons were made with the experimental results of previous researchers. The mean velocityprofile and the turbulent intensity in the experiments indicate that for the mean velocity profile,in the fully developed turbulent pipe flow, von Karman's constant κ is a function of the Reynoldsnumber, i. e. κ increases slowly with the Reynolds number. For turbulent pipe flow, the outer limitdepends on whether the Kdrmdn number R^+ is greater or less than 850 in the centerline velocityprofile: a log law exists for 850 < R^+< 1750 in the experiment, and von Karman' s constant κ isshown to be 0. 408. Under the effects of the test trip at the inlet, fully developed turbulence wasobtained in pipe flow at lower Reynolds number when the entrance length (x/D) was larger than 40. Inthe experiment it was also found that turbulence quantities in pipe flow remain independent of theupstream conditions when the trip blockage ratio is higher than 20%, and the comparison with channelwater flow was also performed.展开更多
In this article,we employ a fully-resolved numerical simulation method(the fictitious domain method)to investigate the effects of large neutrally-buoyant particles on the turbulent flow in a pipe at low Reynolds num...In this article,we employ a fully-resolved numerical simulation method(the fictitious domain method)to investigate the effects of large neutrally-buoyant particles on the turbulent flow in a pipe at low Reynolds number and non-dilute regimes.The tube Reynolds number is fixed to be 4 900,the particle-pipe diameter ratio is 0.1,and the particle volume fraction ranges from 0.33%to 10%.Our results indicate that the presence of large particles decreases the maximum root-of-mean-square(rms)of the streamwise velocity fluctuation near the wall by weakening the intensity of large-scale streamwise vortices,although in the region very close to the wall the particles increase the rms of streamwise velocity fluctuation.On the other hand,the particles induce small-scale vortices in the near-wall region,resulting in the enhancement of the rms of radial and circumferential velocity fluctuations there.展开更多
基金ADNOC Onshore Company(ADCO)for the financial support of this research project.
文摘Turbulent swirling flow inside a short pipe interacting with a conical bluff body was simulated using the commercial CFD code Fluent.The geometry used is a simplified version of a novel liquid/gas separator used in multiphase flow metering.Three turbulence models,belonging to the Reynolds averaged Navier-Stokes(RANS)equations framework,are used.These are,RNG k-ε,SST k-ωand the full Reynolds stress model(RSM)in their steady and unsteady versions.Steady and unsteady RSM simulations show similar behavior.Compared to other turbulence models,they yield the best predictions of the mean velocity profiles though they exhibit some discrepancies in the core region.The influence of the Reynolds number on velocity profiles,swirl decay,and wall pressure on the bluff body are also presented.For Reynolds numbers generating a Rankine-like velocity profile,the width and magnitude of flow reversal zone decreases along the pipe axis disappearing downstream for lower Reynolds numbers.The tangential velocity peaks increase with increasing Reynolds number.The swirl decay rate follows an exponential form in accordance with the existing literature.These flow features would affect the performance of the real separator and,thus,the multiphase flow meter,noticeably.
文摘The motion of fibers in turbulent pipe flow was simulated by 3-D integral method based on the slender body theory and simplified model of turbulence.The orientation distribution of fibers in the computational area for different Re numbers was computed.The results which were consistent with the experimental ones show that the fluctuation velocity of turbulence cause fibers to orient randomly.The orientation distributions become broader as the Re number increases.Then the fluctuation velocity and angular velocity of fibers were obtained.Both are affected by the fluctuation velocity of turbulence.The fluctuation velocity intensity of fiber is stronger at longitudinal than at lateral,while it was opposite for the fluctuation angular velocity intensity of fibers.Finally,the spatial distribution of fiber was given.It is obvious that the fiber dispersion is strenghened with the increase of Re numbers.
文摘To predict the characteristics of dense liquid-solid two-phase flow, K-Ε-T model is established, in which the turbulent flow of fluid phase was described with fluid turbulent kinetic energy Kf and its dissipation rate Εf, and the particles random motion was described with particle turbulent energy Kp and its dissipation rate Εp and pseudothermal temperature Tp. The governing equations were also derived. With K-Ε-T model, numerical study of dense liquid-solid two-phase turbulent up-flow in a pipe is performed. The calculated results are in good agreement with experimental data of Alajbegovic et al. (1994), and some flow features are captured.
文摘In this paper we use the Green function method to solve the problem of steady one dimensional flow of an incompressible viscous, electrically conducting fluid through a pipe with partial circular ring cross sec- tion and one with annular cross section, in the presence of an applied transverse uniform magnetic field, We ob- tain analytic solutions and carry out some numerical calculations of the velocity distribution and induced magnet- ic field.
基金Project supported by the National Natural Science Foundation of China(Nos.11322221,11132005,and 11490551)
文摘The optimal transient growth process of perturbations driven by the pressure gradient is studied in a turbulent pipe flow. A new computational method is proposed, based on the projection operators which project the governing equations onto the sub- space spanned by the radial vorticity and radial velocity. The method is validated by comparing with the previous studies. Two peaks of the maximum transient growth am- plification curve are found at different Reynolds numbers ranging from 20 000 to 250 000. The optimal flow structures are obtained and compared with the experiments and DNS results. The location of the outer peak is at the azimuthal wave number n = 1, while the location of the inner peak is varying with the Reynolds number. It is observed that the velocity streaks in the buffer layer with a spacing of 100δv are the most amplified flow structures. Finally, we consider the optimal transient growth time and its dependence on the azimuthal wave length. It shows a self-similar behavior for perturbations of different scales in the optimal transient growth process.
基金supported by the Croatian Science Foundation(scientific project 3955:Mathematical modeling and numerical simulations of processes in thin or porous domains)
文摘We address the flow of incompressible fluid with a pressure-dependent viscosity through a pipe with helical shape. The viscosity-pressure relation is defined by the Barus law. The thickness of the pipe and the helix step are assumed to be of the same order and considered as the small parameter. After transforming the starting problem, we compute the asymptotic solution using curvilinear coordinates and standard perturbation technique. The solution is provided in the explicit form clearly showing the influence of viscosity-pressure dependence and pipe's geometry on the effective flow.
文摘The nature of turbulent swirling and rotating flow in a straight pipe is investigated using a family of near-wall two-equation models. Specifically, the viability of three different near-wall two-equation models is assessed. These models are asymptotically consistent near the wall. The first two models, one with isotropic and another with anisotropic eddy viscosity invoked, solved a dissipation rate equation with no explicit correction made to account for swirl and flow rotation. The third model assumes an isotropic eddy viscosity but solves an improved dissipation rate equation that has explicit corrections made to account for swirl and flow rotation. Calculations of turbulent flows in the swirl number range 0.25 - 1.3 with and without a central recirculation region reveal that, with the exception of the third model, neither one of the other two models can replicate the mean field at the swirl numbers tested. Furthermore, taking stress anisotropy into account also fails to model swirl effect correctly. Significant improvements can be realized from the third model, which is based on an improved dissipation rate equation and the assumption of isotropic eddy viscosity. The predicted mean flow and turbulence statistics correlate well with measurements at low swirl. At high swirl, the two-equation model with an improved dissipation rate equation can still be used to model swirling and rotating pipe flows with a central recirculation region. However, its simulation of flows without a central recirculation region is not as satisfactory.
文摘The computed orientation distribution of fibers immersed in laminar pipe flows showed that the longitudinal distributions are wide for small Reynolds numbers and become narrower with increasing Re.For low Re number.the axial orientation distributions are broad with almost no preferred orientations.For high Re numer,the axial distribution becomes narrow.with sharp maxima..The mean values of the longitudinal orientation depend strongly on the Re number.The computed results are in qualitative agreemernt with relevant experimental results.
文摘A gate valve is one of the main elements of a circular pipeline, but the flow characteristics around the gate valve are hardly known. In this study, clarification of the flow field in front of the gate valve model in a pipe flow via flow visualization and PIV analysis was attempted. As a result, four kinds of steady necklace-type vortex systems, 2-vortex, 4-vortex, 6-vortex and 8-vortex systems, were clearly observed in a Reynolds number between 290 and 2130. In addition, the main vortex was observed in the Reynolds number range between 2130 and 4870 with difficulty. On this account, both the center position and vorticity in the main vortex are presented against Reynolds number.
文摘The aim of this study is to evaluate a three-equation turbulence model applied to pipe flow. Uncertainty is approximated by comparing with published direct numerical simulation results for fully-developed average pipe flow. The model is based on the Reynolds averaged Navier-Stokes equations. Boussinesq hypothesis is invoked for determining the Reynolds stresses. Three local length scales are solved, based on which the eddy viscosity is calculated. There are two parameters in the model;one accounts for surface roughness and the other is possibly attributed to the fluid. Error in the mean axial velocity and Reynolds stress is found to be negligible.
文摘The Moody Diagram is widely used to determine the friction factor for fluid flow in pipes. The diagram combines the effects of Reynolds number and relative roughness to determine the friction factor. The relationship is highly non-linear and appears to have a complex interaction between viscous and boundary roughness effects. The Moody Diagram is based on predictions from an equation developed by Colebrook in 1939. The relationship requires an iteration process to make predictions. While empirical relationships have been developed that provide good predictions without an iteration process, no one has fully explained the cause for the observed results. The objective of this paper is to present a logical development for prediction of the friction factor. An equation has been developed that models the summed effect of both the laminar sublayer and the boundary roughness on the fluid profile and the resulting friction factor for pipes. The new equation does not require an iteration procedure to obtain values for the friction factor. Predicted results match well with values generated from Colebrook’s work that is expressed in the Moody Diagram. Predictions are within one percent of Colebrook values and generally less than 0.3 percent error from his values. The development provides insight to how processes operating at the boundary cause the friction factor to change.
文摘In this study, the magnetohydrodynamic (MHD) flow through a circular pipe under the influence of a transverse mag- netic field when the outside medium is also electrically conducting is solved numerically by using FEM-BEM coupling approach. The coupled partial differential equations defined for the interior medium are transformed into homogenous modified Helmholtz equations. For the exterior medium on an infinite region, the Laplace equation is considered for the exterior magnetic field. Unknowns in the equations are also related with the corresponding Dirichlet and Neumann type coupled boundary conditions. Unknown values of the magnetic field on the boundary and for the exterior region are obtained by using BEM, and the unknown velocity and magnetic field inside the pipe are obtained by using SUPG type stabilized FEM. Computations are carried for very high values of magnetic Reynolds numbers Rm1, Reynolds number Re and magnetic pressure Rh of the fluid. The results show that using stabilized method enables us to get stable and accurate numerical approximations consistent with the physical configuration of the problem over rough mesh which also results a cheap computational cost.
文摘The steady laminar pipe flow of a suspension with a gas volume fraction ∅≤0.5 and small or intermediate bubble deformations in long, and straight sections of a circular pipe is calculated. The calculations are based on the constitutive equation that was originally derived for dilute emulsions and further developed for concentrated suspensions containing bubbles. In contrast to the literature, an analytical procedure is used to determine the solution of a pipe flow more accurately. The results are presented and discussed with respect to the Reynolds number Re and capillary number Ca. If Ca 1, a bubble suspension has a parabolic velocity profile indicating a Newtonian rheology. If Ca ≈1, two regimes of flow are observed in agreement with the literature;that is, an inner plug flow where deformation rates are low and an outer flow where deformation rates are high. These results imply that, if Ca ∅?and that, if Ca ≥1, the opposite effect occurs;that is, the Reynolds number Re increases with increasing gas volume fraction.
文摘The calculation of the mechanical energy loss is one of the fundamental problems in the field of Hydraulics and Engineering Fluid Mechanics. However, for a non-uniform flow the relation between the mechanical energy loss in a volume of fluid and the kinematical and dynamical characteristics of the flow field is not clearly established. In this paper a new mechanical energy equation for the incompressible steady non-uniform pipe flow of homogeneous fluid is derived, which includes the variation of the mean turbulent kinetic energy, and the formula for the calculation of the mechanical energy transformation loss for the non-uniform flow between two cross sections is obtained based on this equation. This formula can be simplified to the Darcy-Weisbach formula for the uniform flow as widely used in Hydraulics. Furthermore, the contributions of the mechanical energy loss relative to the time averaged velocity gradient and the dissipation of the turbulent kinetic energy in the turbulent uniform pipe flow are discussed, and the contributions of the mechanical energy loss in the viscous sublayer, the buffer layer and the region above the buffer layer for the turbulent uniform flow are also analyzed.
文摘A linear stability analysis on the circular pipe flow of fiber suspensions ispresented. The fiber orientation tensors are used to describe the characteristics of theorientation dislribution of fibers. The constitutive equation for the fiber suspensions is set upand the modified Orr-Sommerfeld stability equation is derived. An improved finite difference methodwith high order accuracy is employed to solve the e-quation. The Newtonian pipe Poiseuille flowcorresponding to H = 0 is also analyzed for comparison. The results reveal that the fiber additiveswill enhance the flow stability, the degree of enhancement becomes high with the increases of theparameter H which accounts for the fiber resistance to the stretching along its axis. Fibersuspensions with large H can suppress the influence of high Reynolds number on the flow stability.Particularly, fibers give a higher attenuation of the short waves of disturbance.
基金support from the National Nature Science Foundation of China(50006005 and 50476010)the INET Research Fund and the 0utstanding Professorship of Chinese National Education Commission.
文摘The deposition of particles in turbulent pipe flow was investigated in terms of two mechanisms, turbulent and thermophoretic. A general equation incorporating these two mechanisms was formulated to calculate the deposition efficiency of aerosol particles in turbulent pipe flow together with thermophoretic deposition. The validity of the equation was confirmed by good agreement between calculated and measured results.
文摘In the non-spherical particulate turbulent flows, a set of new fluidfluctuating velocity equations with the non-spherical particle source term were derived, then a newmethod, which treats the slowly varying functions and rapidly varying functions separately, wasproposed to solve the equations, and finally the turbulent intensity and the Reynolds stress of theflu-id were obtained by calculating the fluctuating velocity statlsti-cally. The equations andmethod were used to a paniculate tur-bulent pipe flow. The results show that the turbulent intensityand the Reynolds stress are decreased almost inverse proportion-ally to the fluctuating velocityratio of particle to fluid. Non-spherical particles have a greater suppressing effect on thetur-bulence than the spherical particles. The particles with short re-laxation time reduce theturbulence intensity of fluid, while the particles with long relaxation time increase the turbulenceinten-sity of fluid. For fixed particle and fluid, the small particles sup-press the turbulence andthe large particles increase the turbu-ience.
文摘Fully developed turbulence measurements in pipe flow were made in theReynolds number range from 10 X 10~3 to 350 X 10~3 with hot-wire anemometer and a Pilot tube.Comparisons were made with the experimental results of previous researchers. The mean velocityprofile and the turbulent intensity in the experiments indicate that for the mean velocity profile,in the fully developed turbulent pipe flow, von Karman's constant κ is a function of the Reynoldsnumber, i. e. κ increases slowly with the Reynolds number. For turbulent pipe flow, the outer limitdepends on whether the Kdrmdn number R^+ is greater or less than 850 in the centerline velocityprofile: a log law exists for 850 < R^+< 1750 in the experiment, and von Karman' s constant κ isshown to be 0. 408. Under the effects of the test trip at the inlet, fully developed turbulence wasobtained in pipe flow at lower Reynolds number when the entrance length (x/D) was larger than 40. Inthe experiment it was also found that turbulence quantities in pipe flow remain independent of theupstream conditions when the trip blockage ratio is higher than 20%, and the comparison with channelwater flow was also performed.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11072217, 10872181)the Fundamental Research Funds for the Central Universities (Grant Nos. 2009QNA4036, 2010QNA4015)
文摘In this article,we employ a fully-resolved numerical simulation method(the fictitious domain method)to investigate the effects of large neutrally-buoyant particles on the turbulent flow in a pipe at low Reynolds number and non-dilute regimes.The tube Reynolds number is fixed to be 4 900,the particle-pipe diameter ratio is 0.1,and the particle volume fraction ranges from 0.33%to 10%.Our results indicate that the presence of large particles decreases the maximum root-of-mean-square(rms)of the streamwise velocity fluctuation near the wall by weakening the intensity of large-scale streamwise vortices,although in the region very close to the wall the particles increase the rms of streamwise velocity fluctuation.On the other hand,the particles induce small-scale vortices in the near-wall region,resulting in the enhancement of the rms of radial and circumferential velocity fluctuations there.