首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 116 毫秒
1.
李勇  钱蔚旻  何录武 《力学季刊》2019,40(3):567-576
挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题.本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象,采用颜色模型模拟出口处粘弹流体和空气的两相流动,通过重新标色获得两种流体的界面,并最终获得胀大的形状.Navier-Stokes方程和本构方程的求解采用双分布函数模型.将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较,发现格子Boltzmann两相模型结果与解析解和实验结果相吻合,相比于单相模型,收敛速度更快,解的稳定性更高.研究了流道尺寸对胀大率的影响,并对挤出胀大的内在机理进行了分析.  相似文献   

2.
管内上随体Maxwell流体非定常流动   总被引:2,自引:1,他引:1  
韩式方  伍岳庆 《力学学报》1990,22(5):519-525
本文研究了上随体Maxwell流体在圆管内非定常流动规律,对于上随体Maxwell流体模型,导出了特殊的运动方程,分别应用隐式差分格式和Kantorovich变分法,求得数值解,对两类方法的结果进行比较,揭示了粘弹流效应对管内非定常流动规津的影响,根据上述研究认为,以上的特殊的变分方法适应于研究非定常流动。  相似文献   

3.
采用模拟黏弹流体挤出胀大的方法,计算了IUPAC-LDPE熔体经过4:1轴对称收缩流道的流动.计算的相对涡强度、入口校正和献中的结果基本一致,给出的流场也显示出计算结果是合理的.表明该方法能够适用于用积分型PSM模型表征的黏弹流体在收缩流道内的流动模拟.  相似文献   

4.
黏弹流体流动的数值模拟研究进展   总被引:5,自引:1,他引:4  
综述了黏弹流体流动数值模拟的研究进展,突出介绍近十年来有限元法在黏弹流体流动数值模拟研究中取得的成果,通过动量方程的适当变形和本构方程离散权函数的合理选择,可以显著增强数值计算的稳定性。得到较高Weissenberg数下的解,同时文中对黏弹流体流动数值模拟中本构方程的应用、非等温情况和三维空间下的研究进行了介绍。  相似文献   

5.
本文采用球形刚性颗粒悬浮于牛顿流体的血液模型,基于有相间滑移的两相流动层流Navier-Stokes方程,采用相间滑移算法(Inter-PhaseSlipAlgorithm)对圆管内定常轴对称血液入口两相流动问题进行了计算,获得了与实验结果吻合的计算结果,较好地模拟血液流动中红细胞的径向迁移现象,结果表明,采用二相流动模型研究血液流动是一种有前途的方法。  相似文献   

6.
非牛顿流体入口收敛流动分析   总被引:5,自引:0,他引:5  
梁基照 《力学学报》1990,22(1):79-85
本文讨论了非牛顿流体的入口收敛流动问题,考虑到粘弹性流体在流动中的粘滑行为,应用最小能原理,导出了扩展的入口收敛流边界流线方程和流体自然收敛锥角方程,并与前人的工作进行了比较和分析。  相似文献   

7.
粘弹性流体的入口收敛流动   总被引:2,自引:0,他引:2  
梁基照 《力学进展》1993,23(2):234-248
本文以聚合物流体为研究对象,对其在入口收敛流动中产生的粘弹效应及机理进行了初步的讨论和分析,并就近10年来国内外有关粘弹性流体入口流动研究及进展作了简要的评述。  相似文献   

8.
将单相格子Boltzmann方法(lattice Boltzmann method,LBM)引入到粘弹流体的瞬态挤出胀大的数值模拟中,建立了基于双分布函数的自由面粘弹性流动格子Boltzmann模型.分析得到的流道中流动速度分布和构型张量结果与理论解十分吻合.对粘弹流体瞬态挤出胀大过程进行了模拟,并分析了运动粘度比和剪切速率对挤出胀大率的影响,得到的胀大率结果与理论分析和其它模拟结果基本一致.表明给出的LBM可以捕捉挤出胀大的瞬态效应.  相似文献   

9.
应用宏观理论的共转模型研究纺丝拉伸过程中测粘-拉伸流动的流变学,以探讨液晶高分子流体纺丝拉伸初步机制。由共转Oldroyd B流体模型的本构方程导出了测粘-拉伸流动的法向应力差和拉伸粘度的解析表达式,利用边界条件及运动方程列出了主拉伸速度与纺程的微分方程。应用Matlab软件绘制出拉伸粘度和主拉伸速度随其它参数变化的曲线,从而得到这些参数对纺丝拉伸流动流变学性能的影响,这对纺丝过程的工艺控制有一定的指导意义。  相似文献   

10.
韩式方 《力学学报》1993,25(2):213-217
提出“准衰退记忆”新概念,发展了非牛顿流体扰动本构理论,并研究了粘弹流体拉伸流动的不稳定性规律  相似文献   

11.
A technique combining the features of parameter differentiation and finite differences is presented to compute the flow of viscoelastic fluids. Two flow problems are considered: (i) three-dimensional flow near a stagnation point and (ii) axisymmetric flow due to stretching of a sheet. Both flows are characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. The exact numerical solutions are obtained using the technique described in the paper. Also, the first-order perturbation solutions (in terms of the viscoelastic fluid parameter) are derived. A comparison of the results shows that the perturbation method is inadequate in predicting some of the vital characteristic features of the flows, which can possibly be revealed only by the exact numerical solution.  相似文献   

12.
Two-dimensional stagnation-point flow of viscoelastic fluids is studied theoretically assuming that the fluid obeys the upper-convected Maxwell (UCM) model. Boundary-layer theory is used to simplify the equations of motion which are further reduced to a single non-linear third-order ODE using the concept of stream function coupled with the technique of the similarity solution. The equation so obtained was solved using Chebyshev pseudo-spectral collocation-point method. Based on the results obtained in the present work, it is concluded that the well-established but controversial prediction that in stagnation-point flows of viscoelastic fluids the velocity inside the boundary layer may exceed that outside the layer may just be an artifact of the rheological model used in previous studies (namely, the second-grade model). No such peculiarity is predicted to exist for the Maxwell model. For a UCM fluid, a thickening of the boundary layer and a drop in wall skin friction coefficient is predicted to occur the higher the elasticity number. These predictions are in direct contradiction with those reported in the literature for a second-grade fluid.  相似文献   

13.
Two‐phase immiscible fluids in a two‐dimensional micro‐channels network are considered. The incompressible Stokes equations are used to describe the Newtonian fluid flow, while the Oldroyd‐B rheological model is used to capture the viscoelastic behavior. In order to perform numerical simulations in a complex geometry like a micro‐channels network, the volume penalization method is implemented. To follow the interface between the two fluids, the level‐set method is used, and the dynamics of the contact line is modeled by Cox law. Numerical results show the ability of the method to simulate two‐phase flows and to follow properly the contact line between the two immiscible fluids. Finally, simulations with realistic parameters are performed to show the difference when a Newtonian fluid is pushed by a viscoelastic fluid instead of a Newtonian one. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

14.
This paper presents a mathematical model for describing approximately the viscoelastic effects in non-Newtonian steady flows through a porous medium. The rheological behaviour of power law fluids is considered in the Maxwell model of elastic behaviour of the fluids. The equations governing the steady flow through porous media are derived and an analytical solution of these equations in the case of a simple flow system is obtained. The conditions for which the viscoelastic effects may become observable from the pressure distribution measurements are shown and expressed in terms of some dimensionless groups. These have been found to be relevant in the evaluation of viscoelastic effects in the steady flow through porous media.  相似文献   

15.
A low cost, low power laser-speckle photographic technique has been developed and is duscussed for the measurement of point velocities in slow laminar flows. The technique is particularly suitable for axisymmetric flows where the two velocity components can be easily measured. The accuracy of the technique is established by measurement of the velocity distribution for Poiseuille flow and from data obtained for acceleration of an inelastic Newtonian fluid through a four-to-one circular contraction. Preliminary results are also presented in the contracting flow field for a non shear-thinning highly elastic fluid. These data are particularly significant for verification of finite element numerical solutions currently being developed for viscoelastic fluids in circular entry flows.  相似文献   

16.
纤维悬浮液搅拌流动的数值模拟   总被引:2,自引:0,他引:2  
由于缺乏适当的本构方程,对纤维悬浮液流动的研究一直局限于纤维的牛顿流体悬浮液。本文采用MUCM模型对作者最近提出的纤维Oldroyd-B流体悬浮液的本构方程作了改进,并对锚式桨搅拌槽的二维Oldroyd-B流体和牛顿流体纤维悬浮液搅拌流动作了数值模拟。模拟的结果表明,本文所用的模型和方法能有效地抑制过大局部应力的影响并合理地处理流体的记忆效应。  相似文献   

17.
Research efforts made so far to simulate fibre suspension flows are limited to fibre suspensions in Newtonian fluids. Though short fibre composites are mostly made of polymers, lack of suitable constitutive equations for fibre suspensions in viscoelastic fluids frustrates attempts to simulate flows of these suspensions. A preliminary work done by the author led to a constitutive equation for semiconcentrated fibre suspensions in the Oldroyd-B fluid. This paper describes the mathematical formulation of the flow problems for such a suspension and the numerical procedure to solve them. Some numerical results of flow past a sphere in a tube are also presented. This paper is supported by the University of Melbourne, the National Natural Science Foundation of China and Zhejiang Province.  相似文献   

18.
For steady non-swirling axisymmetric flow () of an incompressible fluid two invariants of the rate of strain dyadic D are introduced, which directly enter into the expression for D. This being the case they - in conjunction with the vorticity - allow a flow classification into strong and weak flows. For a generalized Newtonian fluid an expression for the viscosity function is listed, which reduces for model fluids to correct results in shearing and, respectively, extensional flow. A possible modification of is proposed, which involves the relative vorticity as well (quasi-Newtonian fluid), since this allows to adjust itself to the local nature of the flow. As such it should prove useful for numerical calculation. Received April 23, 1998  相似文献   

19.
20.
In this work, the linear stability analysis of the viscoelastic Taylor-Couette flow against non-axisymmetric disturbances is investigated. A pseudospectrally generated, generalized algebraic eigenvalue problem is constructed from the linearized set of the three-dimensional governing equations around the steady-state azimuthal solution. Numerical evaluation of the critical eigenvalues shows that for an upper-convected Maxwell model and for the specific set of geometric and kinematic parameters examined in this work, the azimuthal Couette (base) flow becomes unstable against non-axisymmetric time periodic disturbances before it does so for axisymmetric ones, provided the elasticity number ε (De/Re) is larger than some non-zero but small value (ε 0.01). In addition, as ε increases, different families of eigensolutions become responsible for the onset of instability. In particular, the azimuthal wavenumber of the critical eigensolution has been found to change from 1 to 2 to 3 and then back to 2 as ε increases from 0.01 to infinity (inertialess flow).In an analogous fashion to the axisymmetric viscoelastic Taylor-Couette flow, two possible patterns of time-dependent solutions (limit cycles) can emerge after the onset of instability: ribbons and spirals, corresponding to azimuthal and traveling waves, respectively. These patterns are dictated solely by the symmetry of the primary flow and have already been observed in conjunction with experiments involving Newtonian fluids but with the two cylinders counter-rotatng instead of co-rotating as considered here. Inclusion of a non-zero solvent viscosity (Oldroyd-B model) has been found to affect the results quantitatively but not qualitatively. These theoretical predictions are of particular importance for the interpretation of the experimental data obtained in a Taylor-Couette flow using highly elastic viscoelastic fluids.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号