Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though loo...Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though look simple are numerically challenging to study, due to the loss of numerical stability. Over the years, various methodologies have been developed to overcome this numerical limitation. In spite of this, numerical solutions are considered distant from accuracy, as first-order upwind-differencing scheme （UDS） is often employed for improving the stability of algorithm. To elude this effect, some works been reported in the past, where high-resolution-schemes （HRS） were employed and Deborah number was varied. However, these works are limited to creeping flows and do not detail any information on the numerical stability of HRS. Hence, this article presents the numerical study of high shearing contraction flows, where stability of HRS are addressed in reference to fluid elasticity. Results suggest that all I-IRS show some order of undue oscillations in flow variable profiles, measured along vertical lines placed near contraction region in the upstream section of domain, at varied elasticity number E ~ 5. Furthermore, by E, a clear relationship between numerical stability of HRS and E was obtained, which states that the order of undue oscillations in flow variable profiles is directly proportional to E.展开更多
In this paper, the asymmetric laminar flow in a porous channel with expanding or contracting walls is investigated. The governing equations are reduced to ordinary ones by using suitable similar transformations. Homot...In this paper, the asymmetric laminar flow in a porous channel with expanding or contracting walls is investigated. The governing equations are reduced to ordinary ones by using suitable similar transformations. Homotopy analysis method （HAM） is employed to obtain the expres- sions for velocity fields. Graphs are sketched for values of parameters and associated dynamic characteristics, especially the expansion ratio, are analyzed in detail.展开更多
The simulation results on viscoelastic fluid flows in sudden expansion geometry with different expansion ratios are presented. Oldroyd-B, linear Phan-Thien-Tanner （L-PTT） and Finitely Extensible Nonlinear Elastic （...The simulation results on viscoelastic fluid flows in sudden expansion geometry with different expansion ratios are presented. Oldroyd-B, linear Phan-Thien-Tanner （L-PTT） and Finitely Extensible Nonlinear Elastic （FENE-P） based constitutive equations were applied in two-dimensional Cartesian coordinates. The governing equations in transient and fully developed regions were solved using open source software called OpenFOAM. The flow patterns, including velocity profiles, shear stresses and first normal stress differences in some horizontal and vertical sections are illustrated. In addition, effects of the fluid type, flow dynamics and expansion ratio on the flow and vortex patterns in transient and fully developed regions are presented and discussed. The presented results show that existences of vortices cause the inverse velocity and negative stresses in expansion regions of the channel which increase with increment of expansion ratio and Weissenberg number （We）. Furthermore, some dead spaces can be observed at channel expansion regions close to the wall which are also increased. The results also show that at low We numbers all fluids show close behavior while at high We numbers the FENE-P fluid behavior shows high divergence from that of the two other fluids.展开更多
The turbulent flows through the channels with abrupt cross-sectional changes are common and importantphysical process in nature.For a better prediction of the mean flow and turbulent characteristics for this problem,a...The turbulent flows through the channels with abrupt cross-sectional changes are common and importantphysical process in nature.For a better prediction of the mean flow and turbulent characteristics for this problem,atwo-dimensional depth-averaged numerical model is developed.The model is robust and accurate in reproducing therecirculation flow behind a groyne and turbulent flows in channels with abrupt cross-sectional changes,when com-pared to the available experimental data of mean velocities and turbulence kinetic energy.Our results reveal that theabrupt cross-sectional change of a channel can affect the flow pattern significantly and introduces the complex turbu-lence characteristics.In particular,when the channel has an abrupt expansion,the mean flow pattern is mainly in lon-gitudinal direction with rather small transverse component.Meanwhile,a recirculating region forms behind the expan-sion position and the turbulence has very strong intensity within this region.For the flow in the channel with an ab-rupt contraction,the longitudinal component of the flow is decreased by the obstruction on one side and accelerated onthe other side,whereas the transverse velocity is small.The turbulence is extraordinarily strong in the regions adja-cent to the contraction wall in the narrow channel.In both cases of abrupt cross-sectional changes,the TKE is genera-ted dominantly by the shear of the longitudinal velocities.展开更多
Tubing collars’influence on hydrodynamic behavior of annular duct flow has been investigated using Particle Image Velocity(PIV)technology.PIV has become an efficient method for complex transient flows visualization.A...Tubing collars’influence on hydrodynamic behavior of annular duct flow has been investigated using Particle Image Velocity(PIV)technology.PIV has become an efficient method for complex transient flows visualization.A water flow loop with two replaceable variable cross-sections(VCS),75-90 mm and 90-110 mm,in a 129 mm inner diameter(ID)pipe was used.The whole field of the variable cross-section annulus(VCSA)was visualized,including forward-facing step(FFS),narrow annulus(NA),and backwardfacing step(BFS)flow.The VCSA ratio and Reynolds(Re)number influence on streamline distribution,velocity distribution,and turbulence intensity were discussed.Results showed that the recirculation is easier to form in BFS than FFS flow under the same condition.The VCSA ratio affects the formation of recirculation zones and the location of the reattachment point.Reynolds number mainly affects BFS flow by influencing the main velocity component-axial velocity.The turbulence intensity is relatively high in BFS than FFS flow and is larger at y/h>1.0 than y/h<1.0.Furthermore,the streamwise cohere nt structures reveal that the first two modes are predominant and represent the main characteristics of the flow by proper orthogonal decomposition(POD)method.展开更多
Flow patterns in upstream and downstream straight tubes of sudden-changed areas in a horizontal straight pipe were experimentally examined. Both sudden-expansion cross-section (SECS) and sudden-contraction cross-secti...Flow patterns in upstream and downstream straight tubes of sudden-changed areas in a horizontal straight pipe were experimentally examined. Both sudden-expansion cross-section (SECS) and sudden-contraction cross-section (SCCS) were investigated. The flow pattern maps upstream and downstream were delineated and compared with those in straight tubes with uniform cross-sections. The effects of the SECS and SCCS on flow patterns were discussed and analyzed. Furthermore, flow pattern transition mechanisms resulting in occurrences of different flow patterns were simply discussed and some transition criteria for the flow pattern transitions were deduced by using the non-dimensionlized analysis method.展开更多
Deterministic simulation approach was applied to simulating viscoelastic complex flows,in which,the SIMPLE algorithm based on collocated grid was used to solve the conservation equations on the macroscopic level and t...Deterministic simulation approach was applied to simulating viscoelastic complex flows,in which,the SIMPLE algorithm based on collocated grid was used to solve the conservation equations on the macroscopic level and the spectral method was used to solve the Fokker-Planck equation on the mesoscopic level.Here,the transient 4∶1 planar contraction flow for dilute polymer solution was computed by using the coupled technique and the calculated polymeric stress distributions at steady state were identical with the results of continuum approach as well as the corresponding references.Therefore,the presented results indicated that the SIMPLE algorithm based on collocated grid coupled with the spectral method can be used to simulate viscoelastic complex flows effectively.Moreover,the influence of Deborah number and viscosity ratio on polymeric stress was also investigated.展开更多
文摘Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though look simple are numerically challenging to study, due to the loss of numerical stability. Over the years, various methodologies have been developed to overcome this numerical limitation. In spite of this, numerical solutions are considered distant from accuracy, as first-order upwind-differencing scheme （UDS） is often employed for improving the stability of algorithm. To elude this effect, some works been reported in the past, where high-resolution-schemes （HRS） were employed and Deborah number was varied. However, these works are limited to creeping flows and do not detail any information on the numerical stability of HRS. Hence, this article presents the numerical study of high shearing contraction flows, where stability of HRS are addressed in reference to fluid elasticity. Results suggest that all I-IRS show some order of undue oscillations in flow variable profiles, measured along vertical lines placed near contraction region in the upstream section of domain, at varied elasticity number E ~ 5. Furthermore, by E, a clear relationship between numerical stability of HRS and E was obtained, which states that the order of undue oscillations in flow variable profiles is directly proportional to E.
基金supported by the National Natural Science Foundations of China (50936003, 50905013)The Open Project of State Key Lab. for Adv. Matals and Materials (2009Z-02)Research Foundation of Engineering Research Institute of USTB
文摘In this paper, the asymmetric laminar flow in a porous channel with expanding or contracting walls is investigated. The governing equations are reduced to ordinary ones by using suitable similar transformations. Homotopy analysis method （HAM） is employed to obtain the expres- sions for velocity fields. Graphs are sketched for values of parameters and associated dynamic characteristics, especially the expansion ratio, are analyzed in detail.
文摘The simulation results on viscoelastic fluid flows in sudden expansion geometry with different expansion ratios are presented. Oldroyd-B, linear Phan-Thien-Tanner （L-PTT） and Finitely Extensible Nonlinear Elastic （FENE-P） based constitutive equations were applied in two-dimensional Cartesian coordinates. The governing equations in transient and fully developed regions were solved using open source software called OpenFOAM. The flow patterns, including velocity profiles, shear stresses and first normal stress differences in some horizontal and vertical sections are illustrated. In addition, effects of the fluid type, flow dynamics and expansion ratio on the flow and vortex patterns in transient and fully developed regions are presented and discussed. The presented results show that existences of vortices cause the inverse velocity and negative stresses in expansion regions of the channel which increase with increment of expansion ratio and Weissenberg number （We）. Furthermore, some dead spaces can be observed at channel expansion regions close to the wall which are also increased. The results also show that at low We numbers all fluids show close behavior while at high We numbers the FENE-P fluid behavior shows high divergence from that of the two other fluids.
基金supported,in part,by the National Natural Science Foundation of China(51061130547 and51279120)
文摘The turbulent flows through the channels with abrupt cross-sectional changes are common and importantphysical process in nature.For a better prediction of the mean flow and turbulent characteristics for this problem,atwo-dimensional depth-averaged numerical model is developed.The model is robust and accurate in reproducing therecirculation flow behind a groyne and turbulent flows in channels with abrupt cross-sectional changes,when com-pared to the available experimental data of mean velocities and turbulence kinetic energy.Our results reveal that theabrupt cross-sectional change of a channel can affect the flow pattern significantly and introduces the complex turbu-lence characteristics.In particular,when the channel has an abrupt expansion,the mean flow pattern is mainly in lon-gitudinal direction with rather small transverse component.Meanwhile,a recirculating region forms behind the expan-sion position and the turbulence has very strong intensity within this region.For the flow in the channel with an ab-rupt contraction,the longitudinal component of the flow is decreased by the obstruction on one side and accelerated onthe other side,whereas the transverse velocity is small.The turbulence is extraordinarily strong in the regions adja-cent to the contraction wall in the narrow channel.In both cases of abrupt cross-sectional changes,the TKE is genera-ted dominantly by the shear of the longitudinal velocities.
基金financial support of the National Key R&D Program of China(No.2019YFB1504102)the National Natural Science Foundation of China(52020105001)the National Science Fund for Distinguished Young Scholars(51725404)and their approval of publishing this paper。
文摘Tubing collars’influence on hydrodynamic behavior of annular duct flow has been investigated using Particle Image Velocity(PIV)technology.PIV has become an efficient method for complex transient flows visualization.A water flow loop with two replaceable variable cross-sections(VCS),75-90 mm and 90-110 mm,in a 129 mm inner diameter(ID)pipe was used.The whole field of the variable cross-section annulus(VCSA)was visualized,including forward-facing step(FFS),narrow annulus(NA),and backwardfacing step(BFS)flow.The VCSA ratio and Reynolds(Re)number influence on streamline distribution,velocity distribution,and turbulence intensity were discussed.Results showed that the recirculation is easier to form in BFS than FFS flow under the same condition.The VCSA ratio affects the formation of recirculation zones and the location of the reattachment point.Reynolds number mainly affects BFS flow by influencing the main velocity component-axial velocity.The turbulence intensity is relatively high in BFS than FFS flow and is larger at y/h>1.0 than y/h<1.0.Furthermore,the streamwise cohere nt structures reveal that the first two modes are predominant and represent the main characteristics of the flow by proper orthogonal decomposition(POD)method.
基金the National Natural Science Foundation (No.59995460)
文摘Flow patterns in upstream and downstream straight tubes of sudden-changed areas in a horizontal straight pipe were experimentally examined. Both sudden-expansion cross-section (SECS) and sudden-contraction cross-section (SCCS) were investigated. The flow pattern maps upstream and downstream were delineated and compared with those in straight tubes with uniform cross-sections. The effects of the SECS and SCCS on flow patterns were discussed and analyzed. Furthermore, flow pattern transition mechanisms resulting in occurrences of different flow patterns were simply discussed and some transition criteria for the flow pattern transitions were deduced by using the non-dimensionlized analysis method.
文摘Deterministic simulation approach was applied to simulating viscoelastic complex flows,in which,the SIMPLE algorithm based on collocated grid was used to solve the conservation equations on the macroscopic level and the spectral method was used to solve the Fokker-Planck equation on the mesoscopic level.Here,the transient 4∶1 planar contraction flow for dilute polymer solution was computed by using the coupled technique and the calculated polymeric stress distributions at steady state were identical with the results of continuum approach as well as the corresponding references.Therefore,the presented results indicated that the SIMPLE algorithm based on collocated grid coupled with the spectral method can be used to simulate viscoelastic complex flows effectively.Moreover,the influence of Deborah number and viscosity ratio on polymeric stress was also investigated.
