ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged u...ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.展开更多
The objective of the present investigation is to study the vortex-induced vibrations (VIV) for flow past a circular cylinder. The turbulent flow is simulated by using a 2-D standard k-ε model incorporating the fini...The objective of the present investigation is to study the vortex-induced vibrations (VIV) for flow past a circular cylinder. The turbulent flow is simulated by using a 2-D standard k-ε model incorporating the finite volume method (FVM) and the Semi-Implicit Method for the Pressure Linked Equations (SIMPLE) algorithm on non-orthogonal boundary-fitted collocated grids. The wall boundaries are approximated with wall functions. In the numerical cases, the turbulent wake patterns are studied by plotting the streamlines and the turbulent kinetic energy contours. The pressure distributions are investigated. Analyses of the vortex-induced force coefficients and the structural vibrations are carried out. The variations of the Strouhal number with the Reynolds number and of the vortex-induced force coefficients with the reduced velocity are obtained. The results show that this numerical approach is feasible and efficient in investigating the VIV problem for a circular cylinder.展开更多
We present the simulation of the dynamics of fluid-cylinder interactions in a narrow three-dimensional channel filled with a Newtonian fluid, using a Lagrange multiplier based fictitious domain methodology combined wi...We present the simulation of the dynamics of fluid-cylinder interactions in a narrow three-dimensional channel filled with a Newtonian fluid, using a Lagrange multiplier based fictitious domain methodology combined with a finite element method and an operator splitting technique. As expected, a settling truncated cylinder turns its broadside perpendicular to the main stream direction and the center of mass moves to the central axis of the channel. In the case of two truncated cylinders, they first move around each other for a while and then stay together in a "T" shape. After the "T" shape has been formed for a long enough time, we found no vortex shedding behind the cylinders. When simulating the fluidization of 60 truncated cylinders, we captured the features of interactions among fluidized cylinders as observed in experiments.展开更多
A new approach that models lift and drag hydrodynamic force signals operating over cylindrical structures was developed and validated. This approach is based on stochastic auto regressive moving average with exogenous...A new approach that models lift and drag hydrodynamic force signals operating over cylindrical structures was developed and validated. This approach is based on stochastic auto regressive moving average with exogenous (ARMAX) input and its time-varying form, TARMAX. Model structure selection and parameter estimation were discussed while considering the validation stage. In this papel; the cylindrical structure was considered as a dynamic system with an incoming water wave and resulting forces as the input and outputs, respectively. The experimental data, used in this study, were collected from a full-scale rough vertical cylinder at the Delft Hydraulics Laboratory. The practicality of the proposed method and also its efficiency in structural modeling were demonstrated through applying two hydrodynamic force components. For this purpose, an ARMAX model is first used to capture the dynamics of the process, relating in-line forces provided by water waves; secondly, the TARMAX model was applied to modeling and analysis of the lift forces on the cylinder. The evaluation of the lift force by the TARMAX model shows the model is successful in modeling the force from the surface elevation.展开更多
The numerical simulation of extrudate swell is significant in extrusion processing.Precise prediction of extrudate swell is propitious to the control of melt flow and the quality of final products.A mathematical model...The numerical simulation of extrudate swell is significant in extrusion processing.Precise prediction of extrudate swell is propitious to the control of melt flow and the quality of final products.A mathematical model of three-dimensional(3D)viscoelastic flow through elliptical ring die for polymer extrusion was investigated.The penalty function formulation of viscoelastic incompressible fluid was introduced to the finite element model to analyze 3D extrusion problem.The discrete elastic viscous split stress(DEVSS)and streamline-upwind PetrovGalerkin(SUPG)technology were used to obtain stable simulation results.Free surface was updated by updating the streamlines which needs less memory space.According to numerical simulation results,the effect of zero-shear viscosity and elongation parameter on extrudate swell was slight,but with the increase of volumetric flow rate and relax time the extrudate swell ratio increased markedly.Finally,the numerical simulation of extrudate swell flow for low-density polyethylene(LDPE)melts was investigated and the results agreed well with others’work.These conclusions provided quantitative basis for the forecasting extrudate swell ratio and the controlling of extrusion productivity shape.展开更多
The flow past a primary cylinder with one tandem control rod and one staggered control rod is simulated in this paper through solving the Navier-Stokes equations. Two examples are simulated to validate the model, and ...The flow past a primary cylinder with one tandem control rod and one staggered control rod is simulated in this paper through solving the Navier-Stokes equations. Two examples are simulated to validate the model, and the results matched well with those of previous researches. The Reynolds number based on the diameter of the primary cylinder is 500. The diameter ratio between the control rod and the primary cylinder (d/D) is 0.25. It was found that the effect of the combination of one upstream tandem control rod and one staggered control rod on the hydrodynamics of the primary cylinder is a linear superposition of the effect of a corresponding single control rod, and the effect of the upstream tandem control rod is dominant at larger spacing ratios such as G/D=2. For the combination of a downstream tandem control rod and a staggered control rod, the effect of the control rods is different from that of the corresponding single control rod in the region of 0.2〈G/D〈0.5 & 30°〈a〈120° and 0.9〈G/D〈1.4 & 30°〈a〈50°, where the additional effect is obvious. In this case, the effect of the downstream tandem control rod is dominant at small spacing ratios (such as G/D=0.1). At moderate spacing ratios such as G/D=0.4, the effects of the tandem control rod and the staggered control rod are comparable in both cases.展开更多
The optimal design of heating and cooling systems must take into account heat radiation which is a non-linear process.In this study,the mixed convection in a radiative magnetohydrodynamic Eyring-Powell copperwater nan...The optimal design of heating and cooling systems must take into account heat radiation which is a non-linear process.In this study,the mixed convection in a radiative magnetohydrodynamic Eyring-Powell copperwater nanofluid over a stretching cylinder was investigated.The energy balance is modeled,taking into account the non-linear thermal radiation and a thermal slip condition.The effects of the embedded flow parameters on the fluid properties,as well as on the skin friction coefficient and heat transfer rate,are analyzed.Unlike in many existing studies,the recent spectral quasi-linearization method is used to solve the coupled nonlinear boundary-value problem.The computational result shows that increasing the nanoparticle volume fraction,thermal radiation parameter and heat generation parameter enhances temperature profile.We found that the velocity slip parameter and the fluid material parameter enhance the skin friction.A comparison of the current numerical results with existing literature for some limiting cases shows excellent agreement.展开更多
Nonlinear Rossby waves in a Boussinesq fluid model which includes both the vertical and horizontalcomponents of Coriolis force are studied by using the semi-geostrophic approximation and the method of travelling-waves...Nonlinear Rossby waves in a Boussinesq fluid model which includes both the vertical and horizontalcomponents of Coriolis force are studied by using the semi-geostrophic approximation and the method of travelling-wavesolution.Taylor series expansion has been employed to isolate the characteristics of the linear Rossby waves and identifythe Rossby cnoidal and solitary waves.Qualitative analysis indicates that if the disturbances are independent of latitude,the effect of horizontal components of Coriolis force disappears.展开更多
The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding freque...The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...展开更多
The two-dimensional flows around a cylinder between two parallel walls at Re=40 and Re=100 are simulated with computational fluid dynamics(CFD). The governing equations are Navier-Stokes equations. They are discretize...The two-dimensional flows around a cylinder between two parallel walls at Re=40 and Re=100 are simulated with computational fluid dynamics(CFD). The governing equations are Navier-Stokes equations. They are discretized with finite volume method(FVM) and the solution is iterated with PISO Algorithm. Then, the calculating results are compared with the numerical results in literature, and good agreements are obtained. After that, the mechanism of the formation of Karman vortex street is investigated and the instability of the entire flow field is analyzed with the energy gradient theory. It is found that the two eddies attached at the rear of the cylinder have no effect on the flow instability for steady flow, i.e., they don't contribute to the formation of Karman vortex street. The formation of Karman vortex street originates from the combinations of the interaction of two shear layers at two lateral sides of the cylinder and the absolute instability in the cylinder wake. For the flow with Karman vortex street, the initial instability occurs at the region in a vortex downstream of the wake and the center of a vortex firstly loses its stability in a vortex. For pressure driven flow, it is confirmed that the inflection point on the time-averaged velocity profile leads to the instability. It is concluded that the energy gradient theory is potentially applicable to study the flow stability and to reveal the mechanism of turbulent transition.展开更多
Fresh cement mortar is a type of workable paste, which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering. In this paper, Papanastasiou's model for Bingham fluids...Fresh cement mortar is a type of workable paste, which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering. In this paper, Papanastasiou's model for Bingham fluids is solved by using the multiple- relaxation-time lattice Boltzmann model (MRT-LB). Analysis of the stress growth exponent m in Bingham fluid flow simulations shows that Papanastasiou's model provides a good approximation of realistic Bingham plastics for values of m 〉 108. For lower values of m, Papanastasiou's model is valid for fluids between Bingham and Newtonian fluids. The MRT-LB model is validated by two benchmark problems: 2D steady Poiseuille flows and lid-driven cavity flows. Comparing the numerical results of the velocity distributions with corresponding analytical solutions shows that the MRT-LB model is appropriate for studying Bingham fluids while also providing better numerical stability. We further apply the MRT-LB model to simulate flow through a sudden expansion channel and the flow surrounding a round particle. Besides the rich flow structures obtained in this work, the dynamics fhi d force on the round particle is calculated. Results show that both the Reynolds number Re and the Bingham number Bn affect the drag coefficients Co, and a drag coefficient with Re and Bn being taken into account is proposed. The relationship of Bn and the ratio of unyielded zone thickness to particle diameter is also analyzed. Finally, the Bingham fluid flowing around a set of randomly dispersed particles is simulated to obtain the apparent viscosity and velocity fields. These results help simulation of fresh concrete flowing in porous media.展开更多
A numerical method for solving isentropic compressible now problems is presented. This method uses the vorticity and the density as variables. The crux of the method lies in the numerical simulation of the process of ...A numerical method for solving isentropic compressible now problems is presented. This method uses the vorticity and the density as variables. The crux of the method lies in the numerical simulation of the process of vorticity and density variation. The fundamental equation of a compressible discrete vortex method is derived. Unlike the BiotSavart law in the incompressible fluids, the modified Blot-Savart law in compressible fluid should contain the change of the non-solenoidal. The change of density induces the variation of the fluid velocity. The effects of compressibility on flow past a circular cylinder have been investigated by using the compressible discrete vortex method at a Reynolds number of l.E+6. The Mach number is 0.5. The results show that the form of vortex shedding is different from the incompressible now. The separation positions shift upstream, the wake more wide and the street is not clear like the incompressible.展开更多
A three-dimensional large eddy simulation (LES) of a spatially developing round jet is carried out in cylindrical coordinates using a dynamic subgrid model with strong inflow instability. Evolutions of large-scale v...A three-dimensional large eddy simulation (LES) of a spatially developing round jet is carried out in cylindrical coordinates using a dynamic subgrid model with strong inflow instability. Evolutions of large-scale vortex structures represented by tangential vortices are obtained and compared with flow visualization. Also presented are three-dimensional spatial evolutions of coherent structure, which are of quasi two-dimensional Kelvin-Helmholtz instability and vortex rings as well as breaking up of the vortex rings with fully three-dimensional characteristics. Predicted results of mean velocity and turbulent intensity agree well with experiments. They are also compared with the results predicted by LES using standard Smagorinsky model and show good self-similarity. Turbulence spectrum of the predicted velocity shows the -5/3 decay for higher wave number, as expected for turbulent round jet flows. In addition, fl-test and y-test are carded out for the instantaneous velocity, showing that the present LES method can successfully predict the hierarchical structure of round jet.展开更多
The main purpose of this article is to present a mathematical model of ciliary motion in an annulus. In this analysis, the peristaltic motion of non-Newtonian Jeffrey six constant fluid is observed in an annulus with ...The main purpose of this article is to present a mathematical model of ciliary motion in an annulus. In this analysis, the peristaltic motion of non-Newtonian Jeffrey six constant fluid is observed in an annulus with ciliated tips in the presence of heat and mass transfer. The effects of viscous dissipation are also considered. The flow equations of non-Newtonian fluid for the two-dimensional tube in cylindrical coordinates are simplified using the low Reynolds number and long wave-length approximations. The main equations for Jeffrey six constant fluid are considered in cylindrical coordinates system. The resulting nonlinear problem is solved using the regular perturbation technique in terms of a variant of small dimensionless parameter α. The results of the solutions for velocity, temperature and concentration field are presented graphically. B_k is Brinkman number, ST is soret number, and SH is the Schmidth number. Outcome for the longitudinal velocity, pressure rise, pressure gradient and stream lines are represented through graphs. In the history, the viscous-dissipation effect is usually represented by the Brinkman number.展开更多
基金Project(2006BAJ04B04)supported by the National Science and Technology Pillar Program in the Eleventh Five-year Plan PeriodProject(2006AA05Z229)supported by the National High Technology Research and Development Program of China+1 种基金Project supportedby the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education MinistryProject(06wk3023)supported by Hunan Science and Technology Office
文摘ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.
基金The present research is financially supported by the National Natural Science Foundation of China(No.50479027)the Natural Science Foundation of Qingdao,China(No.04-3-jj-01).The support is gratefully appreciated.
文摘The objective of the present investigation is to study the vortex-induced vibrations (VIV) for flow past a circular cylinder. The turbulent flow is simulated by using a 2-D standard k-ε model incorporating the finite volume method (FVM) and the Semi-Implicit Method for the Pressure Linked Equations (SIMPLE) algorithm on non-orthogonal boundary-fitted collocated grids. The wall boundaries are approximated with wall functions. In the numerical cases, the turbulent wake patterns are studied by plotting the streamlines and the turbulent kinetic energy contours. The pressure distributions are investigated. Analyses of the vortex-induced force coefficients and the structural vibrations are carried out. The variations of the Strouhal number with the Reynolds number and of the vortex-induced force coefficients with the reduced velocity are obtained. The results show that this numerical approach is feasible and efficient in investigating the VIV problem for a circular cylinder.
文摘We present the simulation of the dynamics of fluid-cylinder interactions in a narrow three-dimensional channel filled with a Newtonian fluid, using a Lagrange multiplier based fictitious domain methodology combined with a finite element method and an operator splitting technique. As expected, a settling truncated cylinder turns its broadside perpendicular to the main stream direction and the center of mass moves to the central axis of the channel. In the case of two truncated cylinders, they first move around each other for a while and then stay together in a "T" shape. After the "T" shape has been formed for a long enough time, we found no vortex shedding behind the cylinders. When simulating the fluidization of 60 truncated cylinders, we captured the features of interactions among fluidized cylinders as observed in experiments.
文摘A new approach that models lift and drag hydrodynamic force signals operating over cylindrical structures was developed and validated. This approach is based on stochastic auto regressive moving average with exogenous (ARMAX) input and its time-varying form, TARMAX. Model structure selection and parameter estimation were discussed while considering the validation stage. In this papel; the cylindrical structure was considered as a dynamic system with an incoming water wave and resulting forces as the input and outputs, respectively. The experimental data, used in this study, were collected from a full-scale rough vertical cylinder at the Delft Hydraulics Laboratory. The practicality of the proposed method and also its efficiency in structural modeling were demonstrated through applying two hydrodynamic force components. For this purpose, an ARMAX model is first used to capture the dynamics of the process, relating in-line forces provided by water waves; secondly, the TARMAX model was applied to modeling and analysis of the lift forces on the cylinder. The evaluation of the lift force by the TARMAX model shows the model is successful in modeling the force from the surface elevation.
基金Supported by the National Science Foundation for Distinguished Young Scholars of China(50425517) the Shandong Province Natural Science Foundation(Y2007F59)
文摘The numerical simulation of extrudate swell is significant in extrusion processing.Precise prediction of extrudate swell is propitious to the control of melt flow and the quality of final products.A mathematical model of three-dimensional(3D)viscoelastic flow through elliptical ring die for polymer extrusion was investigated.The penalty function formulation of viscoelastic incompressible fluid was introduced to the finite element model to analyze 3D extrusion problem.The discrete elastic viscous split stress(DEVSS)and streamline-upwind PetrovGalerkin(SUPG)technology were used to obtain stable simulation results.Free surface was updated by updating the streamlines which needs less memory space.According to numerical simulation results,the effect of zero-shear viscosity and elongation parameter on extrudate swell was slight,but with the increase of volumetric flow rate and relax time the extrudate swell ratio increased markedly.Finally,the numerical simulation of extrudate swell flow for low-density polyethylene(LDPE)melts was investigated and the results agreed well with others’work.These conclusions provided quantitative basis for the forecasting extrudate swell ratio and the controlling of extrusion productivity shape.
基金the support from the National Natural Science Foundation of China(Nos.11372188,and 51490674)the National Basic Research Program of China(973 Program)(No.2015CB251203)
文摘The flow past a primary cylinder with one tandem control rod and one staggered control rod is simulated in this paper through solving the Navier-Stokes equations. Two examples are simulated to validate the model, and the results matched well with those of previous researches. The Reynolds number based on the diameter of the primary cylinder is 500. The diameter ratio between the control rod and the primary cylinder (d/D) is 0.25. It was found that the effect of the combination of one upstream tandem control rod and one staggered control rod on the hydrodynamics of the primary cylinder is a linear superposition of the effect of a corresponding single control rod, and the effect of the upstream tandem control rod is dominant at larger spacing ratios such as G/D=2. For the combination of a downstream tandem control rod and a staggered control rod, the effect of the control rods is different from that of the corresponding single control rod in the region of 0.2〈G/D〈0.5 & 30°〈a〈120° and 0.9〈G/D〈1.4 & 30°〈a〈50°, where the additional effect is obvious. In this case, the effect of the downstream tandem control rod is dominant at small spacing ratios (such as G/D=0.1). At moderate spacing ratios such as G/D=0.4, the effects of the tandem control rod and the staggered control rod are comparable in both cases.
文摘The optimal design of heating and cooling systems must take into account heat radiation which is a non-linear process.In this study,the mixed convection in a radiative magnetohydrodynamic Eyring-Powell copperwater nanofluid over a stretching cylinder was investigated.The energy balance is modeled,taking into account the non-linear thermal radiation and a thermal slip condition.The effects of the embedded flow parameters on the fluid properties,as well as on the skin friction coefficient and heat transfer rate,are analyzed.Unlike in many existing studies,the recent spectral quasi-linearization method is used to solve the coupled nonlinear boundary-value problem.The computational result shows that increasing the nanoparticle volume fraction,thermal radiation parameter and heat generation parameter enhances temperature profile.We found that the velocity slip parameter and the fluid material parameter enhance the skin friction.A comparison of the current numerical results with existing literature for some limiting cases shows excellent agreement.
基金Supported by National Natural Science Foundation of China under Grant No.40475023
文摘Nonlinear Rossby waves in a Boussinesq fluid model which includes both the vertical and horizontalcomponents of Coriolis force are studied by using the semi-geostrophic approximation and the method of travelling-wavesolution.Taylor series expansion has been employed to isolate the characteristics of the linear Rossby waves and identifythe Rossby cnoidal and solitary waves.Qualitative analysis indicates that if the disturbances are independent of latitude,the effect of horizontal components of Coriolis force disappears.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20060056036)
文摘The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...
基金supported by the Natural Science Foundation of Zhejiang Province LY14E060003 )the Special Major Project of Science and Technology of Zhejiang Province (No.2013C01139)+1 种基金Zhejiang Province Key Science and Technology Innovation Team (2 013TD18)the Science Foundation of Zhejiang Sci-Tech University (No.11130032661215)
文摘The two-dimensional flows around a cylinder between two parallel walls at Re=40 and Re=100 are simulated with computational fluid dynamics(CFD). The governing equations are Navier-Stokes equations. They are discretized with finite volume method(FVM) and the solution is iterated with PISO Algorithm. Then, the calculating results are compared with the numerical results in literature, and good agreements are obtained. After that, the mechanism of the formation of Karman vortex street is investigated and the instability of the entire flow field is analyzed with the energy gradient theory. It is found that the two eddies attached at the rear of the cylinder have no effect on the flow instability for steady flow, i.e., they don't contribute to the formation of Karman vortex street. The formation of Karman vortex street originates from the combinations of the interaction of two shear layers at two lateral sides of the cylinder and the absolute instability in the cylinder wake. For the flow with Karman vortex street, the initial instability occurs at the region in a vortex downstream of the wake and the center of a vortex firstly loses its stability in a vortex. For pressure driven flow, it is confirmed that the inflection point on the time-averaged velocity profile leads to the instability. It is concluded that the energy gradient theory is potentially applicable to study the flow stability and to reveal the mechanism of turbulent transition.
基金supported by the National Key Basic Research Program of China(Grant No.2010CB731504)the Natural Science Foundation of China(Grant Nos.11034010,11272048 and 51239006)+1 种基金European Commission Marie Curie Actions(Grant No.IRSES-294976)the State Key Laboratory of Hydroscience and Engineering(Grant No.2013-KY-2)
文摘Fresh cement mortar is a type of workable paste, which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering. In this paper, Papanastasiou's model for Bingham fluids is solved by using the multiple- relaxation-time lattice Boltzmann model (MRT-LB). Analysis of the stress growth exponent m in Bingham fluid flow simulations shows that Papanastasiou's model provides a good approximation of realistic Bingham plastics for values of m 〉 108. For lower values of m, Papanastasiou's model is valid for fluids between Bingham and Newtonian fluids. The MRT-LB model is validated by two benchmark problems: 2D steady Poiseuille flows and lid-driven cavity flows. Comparing the numerical results of the velocity distributions with corresponding analytical solutions shows that the MRT-LB model is appropriate for studying Bingham fluids while also providing better numerical stability. We further apply the MRT-LB model to simulate flow through a sudden expansion channel and the flow surrounding a round particle. Besides the rich flow structures obtained in this work, the dynamics fhi d force on the round particle is calculated. Results show that both the Reynolds number Re and the Bingham number Bn affect the drag coefficients Co, and a drag coefficient with Re and Bn being taken into account is proposed. The relationship of Bn and the ratio of unyielded zone thickness to particle diameter is also analyzed. Finally, the Bingham fluid flowing around a set of randomly dispersed particles is simulated to obtain the apparent viscosity and velocity fields. These results help simulation of fresh concrete flowing in porous media.
文摘A numerical method for solving isentropic compressible now problems is presented. This method uses the vorticity and the density as variables. The crux of the method lies in the numerical simulation of the process of vorticity and density variation. The fundamental equation of a compressible discrete vortex method is derived. Unlike the BiotSavart law in the incompressible fluids, the modified Blot-Savart law in compressible fluid should contain the change of the non-solenoidal. The change of density induces the variation of the fluid velocity. The effects of compressibility on flow past a circular cylinder have been investigated by using the compressible discrete vortex method at a Reynolds number of l.E+6. The Mach number is 0.5. The results show that the form of vortex shedding is different from the incompressible now. The separation positions shift upstream, the wake more wide and the street is not clear like the incompressible.
基金supported by the National Natural Science Foundation of China (Grant Nos. 50176027 and 50706021)a grant from the Research Committee of The Hong Kong Polytechnic University (Grant No.G-U294)
文摘A three-dimensional large eddy simulation (LES) of a spatially developing round jet is carried out in cylindrical coordinates using a dynamic subgrid model with strong inflow instability. Evolutions of large-scale vortex structures represented by tangential vortices are obtained and compared with flow visualization. Also presented are three-dimensional spatial evolutions of coherent structure, which are of quasi two-dimensional Kelvin-Helmholtz instability and vortex rings as well as breaking up of the vortex rings with fully three-dimensional characteristics. Predicted results of mean velocity and turbulent intensity agree well with experiments. They are also compared with the results predicted by LES using standard Smagorinsky model and show good self-similarity. Turbulence spectrum of the predicted velocity shows the -5/3 decay for higher wave number, as expected for turbulent round jet flows. In addition, fl-test and y-test are carded out for the instantaneous velocity, showing that the present LES method can successfully predict the hierarchical structure of round jet.
文摘The main purpose of this article is to present a mathematical model of ciliary motion in an annulus. In this analysis, the peristaltic motion of non-Newtonian Jeffrey six constant fluid is observed in an annulus with ciliated tips in the presence of heat and mass transfer. The effects of viscous dissipation are also considered. The flow equations of non-Newtonian fluid for the two-dimensional tube in cylindrical coordinates are simplified using the low Reynolds number and long wave-length approximations. The main equations for Jeffrey six constant fluid are considered in cylindrical coordinates system. The resulting nonlinear problem is solved using the regular perturbation technique in terms of a variant of small dimensionless parameter α. The results of the solutions for velocity, temperature and concentration field are presented graphically. B_k is Brinkman number, ST is soret number, and SH is the Schmidth number. Outcome for the longitudinal velocity, pressure rise, pressure gradient and stream lines are represented through graphs. In the history, the viscous-dissipation effect is usually represented by the Brinkman number.