Due to the multiformity and complexity of chain conformation under external flow and the challenge of systematically investigating the transient conformation and dynamic evolution process of polymer chains at the mole...Due to the multiformity and complexity of chain conformation under external flow and the challenge of systematically investigating the transient conformation and dynamic evolution process of polymer chains at the molecular level by means of present experimental techniques,a universal description of both chain conformation and dynamics with respect to continuous volume extensional flow(CVEF)is still absent.Taking into account the temperature effect,we performed dissipative particle dynamics(DPD)simulations with the particles corresponding to the repeat units of polymers over a wide temperature range and analyzed the correlation with the conformational properties of ultra-high molecular weight polyethylene/polypropylene(UHMWPE/PP)blend in response to the CVEF.With time evolution,the polymer chains become highly oriented parallel to the flow direction instead of the initial random coiling and self-aggregation.It is found that a high temperature is necessary for more substantial compactness to take place than low temperature.The low-k plateau and low-k peak in structure factor S(k)curves suggest a low degree of conformational diversity and a high degree of chain stretching.It is also concluded that the intra-molecular C-C bond interaction is the main driving force for the dynamics process of the chain conformations undergoing CVEF,where the motion of the alkyl chains is seriously restricted owing to the increase in bond interaction potential,resulting in a reduction of the difference in diffusion rates among alkyl chains.展开更多
This work systematically simulates the external mass transfer from/to a spherical drop and solid particle suspended in a nonlinear uniaxial extensional creeping flow.The mass transfer problem is governed by three dime...This work systematically simulates the external mass transfer from/to a spherical drop and solid particle suspended in a nonlinear uniaxial extensional creeping flow.The mass transfer problem is governed by three dimensionless parameters:the viscosity ratio(λ),the Peclet number(Pe),and the nonlinear intensity of the flow(E).The existing mass transfer theory,valid for very large Peclet numbers only,is expanded,by numerical simulations,to include a much larger range of Peclet numbers(1≤Pe≤105).The simulation results show that the dimensionless mass transfer rate,expressed as the Sherwood number(5 h),agrees well with the theoretical results at the convection-dominated regime(Pe>103).Only when E>5/4,the simulated Sh for a solid sphere in the nonlinear uniaxial extensional flow is larger than theoretical results because the theory neglects the effect of the vortex formed outside the particle on the rate of mass transfer.Empirical correlations are proposed to predict the influence of the dimensionless governing parameters(λ,Pe,E)on the Sherwood number(Sh).The maximum deviations of all empirical correlations are less than 15%when compared to the numerical simulated results.展开更多
The internal flow of a droplet in the nonlinear extensional flow field will exhibit more than two internal circulations with the variation of nonlinear intensity(E).In this paper,the effect of positions and sizes of i...The internal flow of a droplet in the nonlinear extensional flow field will exhibit more than two internal circulations with the variation of nonlinear intensity(E).In this paper,the effect of positions and sizes of internal circulations on internal mass transfer rate of a single spherical droplet in a nonlinear extensional flow field is studied and compared with that in a linear extensional flow field.The simulation results show that when E≥0,there are two symmetrical internal circulations in the droplet,which is the same with that in a linear extensional flow.The limit value of mass transfer rate Sh is 15,which is equal to that in a linear extensional flow,no matter how large E is.When E≤-3/7,the number of internal flow circulation of a droplet increase to four and the transfer rate Sh increases.When E=-1,the maximum internal transfer rate Sh equals 30 which is twice of that in a linear extensional flow.The generation of new flow circulations in droplets and the circulation positions will enhance mass transfer when E≤-3/7,which provides a new idea for enhancing the internal mass transfer rate of droplets.展开更多
The asymmetric breakups of a droplet in an axisymmetric cross-like microfluidic device are investigated by using a three-dimensional volume of fluid(VOF) multiphase numerical model. Two kinds of asymmetries(droplet lo...The asymmetric breakups of a droplet in an axisymmetric cross-like microfluidic device are investigated by using a three-dimensional volume of fluid(VOF) multiphase numerical model. Two kinds of asymmetries(droplet location deviation from the symmetric geometry center and different flow rates at two symmetric outlets) generate asymmetric flow fields near the droplet, which results in the asymmetric breakup of the latter. Four typical breakup regimes(no breakup, one-side breakup, retraction breakup and direct breakup) have been observed.Two regime maps are plotted to describe the transition from one regime to another for the two types of different asymmetries, respectively. A power law model, which is based on the three critical factors(the capillary number,the asymmetry of flow fields and the initial volume ratio), is employed to predict the volume ratio of the two unequal daughter droplets generated in the direct breakup. The influences of capillary numbers and the asymmetries have been studied systematically in this paper. The larger the asymmetry is, the bigger the oneside breakup zone is. The larger the capillary number is, the more possible the breakup is in the direct breakup zone. When the radius of the initial droplet is 20 μm, the critical capillary numbers are 0.122, 0.128, 0.145,0.165, 0.192 and 0.226 for flow asymmetry factor AS= 0.05, 0.1, 0.2, 0.3, 0.4 and 0.5, respectively, in the flow system whose asymmetry is generated by location deviations. In the flow system whose asymmetry is generated by two different flow rates at two outlets, the critical capillary numbers are 0.121, 0.133, 0.145, 0.156 and 0.167 for AS= 1/21, 3/23, 1/5, 7/27 and 9/29, respectively.展开更多
A new continuum theory of the constitutive equation of co-rotational derivative type was developed by the author for anisotropic viscoelastic fluid-liquid crystalline (LC) polymers (S.F. Han, 2008, 2010) . This paper ...A new continuum theory of the constitutive equation of co-rotational derivative type was developed by the author for anisotropic viscoelastic fluid-liquid crystalline (LC) polymers (S.F. Han, 2008, 2010) . This paper is a continuation of the recent publication [1] to study extrusion-extensional flow of the fluid. A new concept of simple anisotropic fluid is introduced. On the basis of anisotropic simple fluid, stress behavior is described by velocity gradient tensor F and spin tensor W instead of the velocity gradient tensor D in the classic Leslie?Ericksen continuum theory. A special form of the constitutive equation of the co-rotational type is established for the fluid. Using the special form of the constitutive equation in components a computational analytical theory of the extrusion-extensional flow is developed for the LC polymer liquids - anisotropic viscoelastic fluid. Application of the constitutive theory to the flow is successful in predicting bifurcation of elongational viscosity and contraction of extrudate for LC polymer liquids–anisotropic viscoelastic fluid. The contraction of extrudate of LC polymer liquids may be associated with the stored elastic energy conversion into that necessary for bifurcation of elongational viscosity in extrusion extensional flow of the fluid.展开更多
The extensional viscosity of three polyethylenes(LDPE,LLDPE and HDPE)melts was measured in-uniaxialelongational flow at a constant extension rate,as well asby the process of fiber spinning"Rheotens"and the e...The extensional viscosity of three polyethylenes(LDPE,LLDPE and HDPE)melts was measured in-uniaxialelongational flow at a constant extension rate,as well asby the process of fiber spinning"Rheotens"and the entryflow methods.The melts showed the same melt flowingindex(MFI)but differed in their extensional viscositiesand attenuating forces,reflecting the differences in theirmolecular structures.The indirect methods could alsogive meaningful estimates for the extensional viscosity ofthe polymer melts,if the amount of the accumulated ex-tensional strain in each case was taken into consider-ation.展开更多
Recent numerical studies on pressure-drops in contraction flows have introduced a variety of constitutive models to compare and contrast the competing influences of extensional vis-cosity, normal stress and shear-thin...Recent numerical studies on pressure-drops in contraction flows have introduced a variety of constitutive models to compare and contrast the competing influences of extensional vis-cosity, normal stress and shear-thinning. Early work on pressure-drops employed the constant viscosity Oldroyd-B and Upper Convected Max- well (UCM) models to represent the behavior of so-called Boger fluids in axisymmetric contrac-tion flows, in (unsuccessful) attempts to predict the very large enhancements that were ob-served experimentally. In more recent studies, other constitutive models have been employed to interpret observed behavior and some pro-gress has been made, although finding a (re-spectable) model to describe observed contrac-tion-flow behavior, even qualitatively, has been frustratingly difficult. With this in mind, the present study discusses the ability of a well- known FENE type model (the so-called FENE- CR model) to describe observed behavior. For various reasons, an axisymmetric (4:1:4) con-traction/expansion geometry, with rounded corners, is singled out for special attention, and a new hybrid finite element/volume algo-rithm is utilized to conduct the modeling, which reflects an incremental pressure-correction time-stepping structure. New to this algo-rithmic formulation are techniques in time discretization, discrete treatment of pressure terms, and compatible stress/velocity-gradient representation. We shall argue that the current simulations for the FENE-CR model have re-sulted in a major improvement in the sort-for agreement between theory and experiment in this important bench-mark problem.展开更多
基金the National Key R&D Program of China(No.2016YFB0302301)the Guangdong YangFan Innovative&Ente preneurial Research TeamProgram(No.2016YT03C077)+1 种基金the Science and Technology Planning Project of Guangzhou(No.201704020008)the Open Foundation of Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanics(Wuhan University of Technology)(No.TAM202001)。
文摘Due to the multiformity and complexity of chain conformation under external flow and the challenge of systematically investigating the transient conformation and dynamic evolution process of polymer chains at the molecular level by means of present experimental techniques,a universal description of both chain conformation and dynamics with respect to continuous volume extensional flow(CVEF)is still absent.Taking into account the temperature effect,we performed dissipative particle dynamics(DPD)simulations with the particles corresponding to the repeat units of polymers over a wide temperature range and analyzed the correlation with the conformational properties of ultra-high molecular weight polyethylene/polypropylene(UHMWPE/PP)blend in response to the CVEF.With time evolution,the polymer chains become highly oriented parallel to the flow direction instead of the initial random coiling and self-aggregation.It is found that a high temperature is necessary for more substantial compactness to take place than low temperature.The low-k plateau and low-k peak in structure factor S(k)curves suggest a low degree of conformational diversity and a high degree of chain stretching.It is also concluded that the intra-molecular C-C bond interaction is the main driving force for the dynamics process of the chain conformations undergoing CVEF,where the motion of the alkyl chains is seriously restricted owing to the increase in bond interaction potential,resulting in a reduction of the difference in diffusion rates among alkyl chains.
基金support and helpful insight.This work was supported by the National Key Research and Development Program(2021YFC2902502)the National Natu-ral Science Foundation of China(21938009,91934301,22078320)+5 种基金the Major Scientific and Technological Innovation Projects in Shan-dong Province(2019JZZY010302)the Shandong Key Research and Development Program(International Cooperation Office)(2019GHZ018)the Shandong Province Postdoctoral Innovative Talents Support Plan(SDBX2020018)the External Cooperation Program of BIC,Chinese Academy of Sciences(122111KYSB20190032)Chemistry and Chemical Engineering Guangdong Laboratory(1922006)GHfund B(202107021062).
文摘This work systematically simulates the external mass transfer from/to a spherical drop and solid particle suspended in a nonlinear uniaxial extensional creeping flow.The mass transfer problem is governed by three dimensionless parameters:the viscosity ratio(λ),the Peclet number(Pe),and the nonlinear intensity of the flow(E).The existing mass transfer theory,valid for very large Peclet numbers only,is expanded,by numerical simulations,to include a much larger range of Peclet numbers(1≤Pe≤105).The simulation results show that the dimensionless mass transfer rate,expressed as the Sherwood number(5 h),agrees well with the theoretical results at the convection-dominated regime(Pe>103).Only when E>5/4,the simulated Sh for a solid sphere in the nonlinear uniaxial extensional flow is larger than theoretical results because the theory neglects the effect of the vortex formed outside the particle on the rate of mass transfer.Empirical correlations are proposed to predict the influence of the dimensionless governing parameters(λ,Pe,E)on the Sherwood number(Sh).The maximum deviations of all empirical correlations are less than 15%when compared to the numerical simulated results.
基金supported by the National Key Research and Development Program of China(2021YFC2902502)the National Natural Science Foundation of China(22078320,22035007)+4 种基金the NSFC-EU project(31961133018)the Shandong Provincial Key Research and Development Program(2022CXGC020106)the Shandong Key Research and Development Program(International Cooperation Office)(2019GHZ018)the Shandong Province Postdoctoral Innovative Talents Support Plan(SDBX2020018)the External Cooperation Program of BIC,Chinese Academy of Sciences(122111KYSB20190032)。
文摘The internal flow of a droplet in the nonlinear extensional flow field will exhibit more than two internal circulations with the variation of nonlinear intensity(E).In this paper,the effect of positions and sizes of internal circulations on internal mass transfer rate of a single spherical droplet in a nonlinear extensional flow field is studied and compared with that in a linear extensional flow field.The simulation results show that when E≥0,there are two symmetrical internal circulations in the droplet,which is the same with that in a linear extensional flow.The limit value of mass transfer rate Sh is 15,which is equal to that in a linear extensional flow,no matter how large E is.When E≤-3/7,the number of internal flow circulation of a droplet increase to four and the transfer rate Sh increases.When E=-1,the maximum internal transfer rate Sh equals 30 which is twice of that in a linear extensional flow.The generation of new flow circulations in droplets and the circulation positions will enhance mass transfer when E≤-3/7,which provides a new idea for enhancing the internal mass transfer rate of droplets.
基金Supported by Major State Basic Research Development Program of China(2012CB720305)the National Natural Science Foundation of China(21376162)
文摘The asymmetric breakups of a droplet in an axisymmetric cross-like microfluidic device are investigated by using a three-dimensional volume of fluid(VOF) multiphase numerical model. Two kinds of asymmetries(droplet location deviation from the symmetric geometry center and different flow rates at two symmetric outlets) generate asymmetric flow fields near the droplet, which results in the asymmetric breakup of the latter. Four typical breakup regimes(no breakup, one-side breakup, retraction breakup and direct breakup) have been observed.Two regime maps are plotted to describe the transition from one regime to another for the two types of different asymmetries, respectively. A power law model, which is based on the three critical factors(the capillary number,the asymmetry of flow fields and the initial volume ratio), is employed to predict the volume ratio of the two unequal daughter droplets generated in the direct breakup. The influences of capillary numbers and the asymmetries have been studied systematically in this paper. The larger the asymmetry is, the bigger the oneside breakup zone is. The larger the capillary number is, the more possible the breakup is in the direct breakup zone. When the radius of the initial droplet is 20 μm, the critical capillary numbers are 0.122, 0.128, 0.145,0.165, 0.192 and 0.226 for flow asymmetry factor AS= 0.05, 0.1, 0.2, 0.3, 0.4 and 0.5, respectively, in the flow system whose asymmetry is generated by location deviations. In the flow system whose asymmetry is generated by two different flow rates at two outlets, the critical capillary numbers are 0.121, 0.133, 0.145, 0.156 and 0.167 for AS= 1/21, 3/23, 1/5, 7/27 and 9/29, respectively.
文摘A new continuum theory of the constitutive equation of co-rotational derivative type was developed by the author for anisotropic viscoelastic fluid-liquid crystalline (LC) polymers (S.F. Han, 2008, 2010) . This paper is a continuation of the recent publication [1] to study extrusion-extensional flow of the fluid. A new concept of simple anisotropic fluid is introduced. On the basis of anisotropic simple fluid, stress behavior is described by velocity gradient tensor F and spin tensor W instead of the velocity gradient tensor D in the classic Leslie?Ericksen continuum theory. A special form of the constitutive equation of the co-rotational type is established for the fluid. Using the special form of the constitutive equation in components a computational analytical theory of the extrusion-extensional flow is developed for the LC polymer liquids - anisotropic viscoelastic fluid. Application of the constitutive theory to the flow is successful in predicting bifurcation of elongational viscosity and contraction of extrudate for LC polymer liquids–anisotropic viscoelastic fluid. The contraction of extrudate of LC polymer liquids may be associated with the stored elastic energy conversion into that necessary for bifurcation of elongational viscosity in extrusion extensional flow of the fluid.
文摘The extensional viscosity of three polyethylenes(LDPE,LLDPE and HDPE)melts was measured in-uniaxialelongational flow at a constant extension rate,as well asby the process of fiber spinning"Rheotens"and the entryflow methods.The melts showed the same melt flowingindex(MFI)but differed in their extensional viscositiesand attenuating forces,reflecting the differences in theirmolecular structures.The indirect methods could alsogive meaningful estimates for the extensional viscosity ofthe polymer melts,if the amount of the accumulated ex-tensional strain in each case was taken into consider-ation.
文摘Recent numerical studies on pressure-drops in contraction flows have introduced a variety of constitutive models to compare and contrast the competing influences of extensional vis-cosity, normal stress and shear-thinning. Early work on pressure-drops employed the constant viscosity Oldroyd-B and Upper Convected Max- well (UCM) models to represent the behavior of so-called Boger fluids in axisymmetric contrac-tion flows, in (unsuccessful) attempts to predict the very large enhancements that were ob-served experimentally. In more recent studies, other constitutive models have been employed to interpret observed behavior and some pro-gress has been made, although finding a (re-spectable) model to describe observed contrac-tion-flow behavior, even qualitatively, has been frustratingly difficult. With this in mind, the present study discusses the ability of a well- known FENE type model (the so-called FENE- CR model) to describe observed behavior. For various reasons, an axisymmetric (4:1:4) con-traction/expansion geometry, with rounded corners, is singled out for special attention, and a new hybrid finite element/volume algo-rithm is utilized to conduct the modeling, which reflects an incremental pressure-correction time-stepping structure. New to this algo-rithmic formulation are techniques in time discretization, discrete treatment of pressure terms, and compatible stress/velocity-gradient representation. We shall argue that the current simulations for the FENE-CR model have re-sulted in a major improvement in the sort-for agreement between theory and experiment in this important bench-mark problem.