非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

弹性湍流及其对微混合效率影响研究的若干进展

介绍和评述弹性湍流的产生及其对于微混合效率的影响等问题上的若干研究进展. 弹性 湍流和惯性湍流具有类似的流场特征,但引发机理有所不同. 惯性湍流产生的原因 是惯性引起的Reynolds应力,而弹性湍流则是由弹性应力所引起的. 鉴于在微流动 中,惯性力可忽略不计,因此牛顿流体的混合变得十分困难. 此时可在流体中加入 微量高分子聚合物以生成黏弹性流体. 由黏弹性流体所引发的弹性湍流在提高流体 微混合的效率上可发挥重要作用.     

基于分形导数对非牛顿流体层流的数值研究

非牛顿流体具有复杂的流变特性,揭示该流变特性可以更加合理地指导非牛顿流体在工农业生产中的应用.经典的非牛顿流体本构模型往往形式复杂,仅能应用于某些特定的情况.分数阶导数模型具有参数少和形式简单的特点,己成功地应用于描述非牛顿流体的运动.Hausdorff分形导数作为一个备选的建模方法,相比分数阶导数具有更简单的形式以及更高的计算效率.本文基于Hausdorff分形导数改进现有牛顿黏性模型,提出分形黏壶模型.通过研究分形黏壶在常应变率下表观黏度的变化情况,以及在加、卸载条件下的蠕变及恢复特性,发现分形黏壶模型适合于描述具有黏弹性的非牛顿流体(本文称之为分形流体).结合连续性方程及运动微分方程,推导出分形流体在平行板间层流的基本方程.按是否拖动上板和是否存在水平的压力梯度分为3种工况,分别用数值方法计算这3种工况下流速在板间的分布及其随时间变化的情况.通过分析不同工况下的流速分布,发现水平的压力梯度会改变流速随时间变化的形状,且会推迟流速到达稳定的时间.在水平压力梯度不存在的情况下,不同阶数的分形流体具有相同的流速分布或是演变过程.另外,在水平压力梯度存在的情况下,上板速度不影响不同阶数分形流体间稳定速度的差值.     

颗粒在黏弹性流体扩张和收缩流中沉降的格子Boltzmann模拟分析

对于Oldroyd-B型黏弹性流体,本文应用格子Boltzmann方法(LBM),实现了流体在二维1:3扩张流道及3:1收缩流道中流动的数值模拟,获得了黏弹性流体在扩张和收缩流道中的流场分布．结合颗粒的受力和运动规则,基于点源颗粒模型,数值分析了颗粒在扩张流和收缩流中的沉降过程和特征,讨论了颗粒相对质量和起始位置以及雷诺数Re和威森伯格数Wi对颗粒沉降特征的影响．结果表明,颗粒相对质量和起始位置以及Re对颗粒沉降轨迹和落点影响较大,而Wi的影响则较小．     

复杂流动中的智能颗粒游动策略研究进展

智能颗粒定义为可以根据环境的变化而主动调整自身在流场中运动的颗粒, 一般用于描述微小的游动物体, 如微生物、浮游生物和微纳机器人. 由于颗粒运动学特性和流场环境的复杂性, 探索智能颗粒的运动策略是一项具有现实意义与挑战性的研究. 近期强化学习算法被应用于智能颗粒的运动策略研究, 并取得一定进展. 本文将讨论强化学习在智能颗粒研究中的应用, 并介绍浮游生物运动策略的相关研究进展, 包括适用于海洋浮游生物的游动颗粒模型, 以及强化学习的基本原理和此类问题的研究框架. 海洋浮游生物的昼夜垂直迁移对其生存和繁衍至关重要, 生物学研究表明某些浮游生物可感知周围的流场信息, 但能否利用这些信息提高垂直方向游动效率仍是一个未知的问题. 基于这一背景, 相关研究考虑了重力沉降和颗粒形状对浮游生物垂直游动策略的影响. 研究发现细长的颗粒能够更加高效地垂直向上运动, 而重力沉降则导致智能游动策略发生较大变化. 在此基础上, 后续研究进一步考虑了局部流体信号的作用, 并讨论了基于局部信息在全局坐标系中定向运动的可能性. 当颗粒只能感受到局部流体信号时, 必须破坏动力学系统的对称性, 否则颗粒无法学到有效的垂向游动策略. 此外, 研究还发现强化学习能够利用流体信号背后的机制, 得到颗粒在二维定常流动和三维非定常湍流中高效垂直运动的策略. 这些智能游动策略依赖于本质的物理规律, 因此这些策略在更加复杂或真实的流动中也可能有优秀的表现.      

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.      

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型,介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展.通过线性稳定性理论,分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, DarcyBrinkman-Oldroyd以及Darcy-Brinkman-Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响.利用弱非线性分析方法,揭示失稳临界点附近热对流流动的分叉情况,以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式.采用数值模拟方法,研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的,而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞,最后发展为混沌状态.     

Oldroyd-B黏弹性液滴碰撞过程的数值模拟

复杂的流变特性使凝胶推进剂的雾化过程存在一定困难,这制约了它的发展.聚合物胶凝剂的加入使凝胶推进剂具有黏弹性,从而在雾化时会产生黏弹性液滴,因此为了进一步认识凝胶推进剂的雾化机理、提高凝胶推进剂的雾化性能,对黏弹性液滴的碰撞行为进行数值模拟研究.针对凝胶推进剂雾化过程中出现的液滴撞击现象,考虑流体具有的黏弹性效应,采用...     

表面活性剂流体弹性应力对颗粒沉降的影响

颗粒在粘弹性表面活性剂流体中的沉降问题是复杂固液两相流研究的基础问题。采用实验和数值模拟的方法研究了单颗粒在粘弹性表面活性剂流体中的沉降过程。实验结果表明,粘弹性表面活性剂流体的粘弹性会引起颗粒的终端速度随时间产生波动。Giesekus本构方程可以描述颗粒在表面活性剂流体中沉降的非线性剪切变稀行为和弹性导致的拉伸变形。颗粒在粘弹性表面活性剂流体中沉降时其尾部会出现负尾迹现象,负尾迹区随着松弛时间的增大向下游延伸,表征流体弹性的松弛时间越大,负尾迹区越长。颗粒表面剪切速率大于临界剪切速率时,负尾迹的出现使颗粒在粘弹性表面活性剂流体中沉降所受曳力减小,沉降速度增加;颗粒表面剪切速率小于临界剪切速率时,粘弹性表面活性剂流体负尾迹区的拉伸效应与颗粒表面的弹性应力的共同作用,导致颗粒终端速度随着松弛时间的增加先增大后减小。     

微流道内椭圆粒子黏弹性聚焦数值模拟研究

黏弹性聚焦技术借助微尺度黏弹性流体的惯性和弹性耦合效应,能够实现生物粒子在流道中心的单一位置聚焦排列,被认为是未来生物粒子计数以及检测的理想预处理单元,因而引起了广泛的关注．自然界中的生物粒子往往是非球形的,故而研究不同形状粒子在黏弹性流体中的迁移特性具有十分重要的价值．本文通过格子玻尔兹曼方法耦合浸入边界法,对椭球粒子在直流道内黏弹性流体中的聚焦行为进行了系统的数值模拟研究．结果表明,面积相同但长径比不同的椭圆粒子在黏弹性流体中有不同的旋转周期与迁移速度．长径比更大的粒子旋转周期更长,且长径比大于3.5 的粒子甚至不再有明显的旋转．长径比更大的粒子上下两侧的黏弹性力分布更加平缓,受到指向流道中心的弹性力更小,使得粒子横向迁移速度更慢从而导致了长径比不同的椭圆粒子聚焦至流道中心所需时间的差异．此外,Weissenberg 数Wi 的增加同样能够减弱粒子的旋转,使得长径比稍小的粒子也能和长径比为1.0 的圆形粒子产生明显的分离．上述数值模拟的结论,为不同长径比粒子在黏弹性流体中的聚焦与分选应用提供了重要的理论指导．     

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall considering the interaction between squirmer and wall are numerically studied with an immersed boundary-lattice Boltzmann method. The power-law index, Reynolds number, initial orientation angle of squirmer, and initial distance of squirmer from the wall are all taken into account to investigate the swimming characteristics for pusher (β?<?0), neutral squirmer (β?=?0), and puller (β?>?0) (three kinds of swimmer types) near the no-slip boundary. Four new kinds of swimming modes are found. Results show that, for the pushers and pullers, the wall displays an increasing attraction with increasing power-law index n, which differs from the neutral squirmer who always departs from the wall after the first collision with the wall. Both the initial orientation angle and initial distance from the wall only affect the moving situations rather than the moving modes of the squirmers. However, the squirmers depart from the wall as the Reynolds number increases and chaotic orbits appear for some squirmers at Re?=?5. Several typical flow fields are analyzed and the power consumption and torque for different kinds of flows are also studied. It is found that, as the absolute value of β increases, the power consumption generally increases in shear-thinning (n?=?0.4), Newtonian (n?=?1), and shear-thickening (n?=?1.6) fluids. Moreover, the pushers (β?<?0) and the pullers (β?>?0) expend almost the same power if the absolute value of β remains the same. In addition, the power consumption of the squirmers is highly dependent on the power-law index n.     

Viscoelastic flow analysis using the software OpenFOAM and differential constitutive equations

Viscoelastic fluids are of great importance in many industrial sectors, such as in food and synthetic polymers industries. The rheological response of viscoelastic fluids is quite complex, including combination of viscous and elastic effects and non-linear phenomena. This work presents a numerical methodology based on the split-stress tensor approach and the concept of equilibrium stress tensor to treat high Weissenberg number problems using any differential constitutive equations. The proposed methodology was implemented in a new computational fluid dynamics (CFD) tool and consists of a viscoelastic fluid module included in the OpenFOAM, a flexible open source CFD package. Oldroyd-B/UCM, Giesekus, Phan-Thien–Tanner (PTT), Finitely Extensible Nonlinear Elastic (FENE-P and FENE-CR), and Pom–Pom based constitutive equations were implemented, in single and multimode forms. The proposed methodology was evaluated by comparing its predictions with experimental and numerical data from the literature for the analysis of a planar 4:1 contraction flow, showing to be stable and efficient.     

A numerical study on nonlinear dynamics of three-dimensional time-depended viscoelastic Taylor-Couette flow

In this paper, three-dimensional viscoelastic Taylor-Couette instability between concentric rotating cylinders is studied numerically. The aim is to investigate and provide additional insight about the formation of time-dependent secondary flows in viscoelastic fluids between rotating cylinders. Here, the Giesekus model is used as the constitutive equation. The governing equations are solved using the finite volume method (FVM) and the PISO algorithm is employed for pressure correction. The effects of elasticity number, viscosity ratio, and mobility factor on various instability modes (especially high order ones) are investigated numerically and the origin of Taylor-Couette instability in Giesekus fluids is studied using the order of magnitude technique. The created instability is simulated for large values of fluid elasticity and high orders of nonlinearity. Also, the effect of elastic properties of fluid on the time-dependent secondary flows such as wave family and traveling wave and also on the critical conditions are studied in detail.     

Numerical study of Saffman–Taylor instability in immiscible nonlinear viscoelastic flows

In this paper, a numerical solution for Saffman–Taylor instability of immiscible nonlinear viscoelastic-Newtonian displacement in a Hele–Shaw cell is presented. Here, a nonlinear viscoelastic fluid pushes a Newtonian fluid and the volume of fluid method is applied to predict the formation of two phases. The Giesekus model is considered as the constitutive equation to describe the nonlinear viscoelastic behavior. The simulation is performed by a parallelized finite volume method (FVM) using second order in both the spatial and the temporal discretization. The effect of rheological properties and surface tension on the immiscible Saffman–Taylor instability are studied in detail. The destabilizing effect of shear-thinning behavior of nonlinear viscoelastic fluid on the instability is studied by changing the mobility factor of Giesekus model. Results indicate that the fluid elasticity and capillary number decrease the intensity of Saffman–Taylor instability.     

Numerical simulation of delayed die swell

In a recent paper, Joseph et al. showed that, for a number of viscoelastic fluids, one can observe the phenomenon of delayed die swell beyond a critical extrusion velocity, or beyond a critical value of the viscoelastic Mach number. Giesekus had also observed that delayed die swell is a critical phenomenon.In the present paper, we find a set of material and flow parameters under which it is possible to simulate delayed die swell. For the viscoelastic flow calculation, we use the finite element algorithm with sub-elements for the stresses and streamline upwinding in the discretized constitutive equations. For the free surface, we use an implicit technique which allows us to implement Newton's method for solving the non-linear system of equations. The fluid is Oldroyd-B which, in the present problem, is a singular perturbation of the Maxwell fluid. The results show very little sensitivity to the size of the retardation time. We also show delayed die swell for a Giesekus fluid.This paper is dedicated to Professor Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.     

Comparison of a quasi-Newtonian fluid with a viscoelastic fluid in planar contraction flow

The flow of a 5.0 wt.% solution of polyisobutylene in tetradecane through a planar 4 : 1 contraction exhibiting a shear thinning viscosity is simulated using the flow-type sensitive quasi-Newtonian fluid model. The shear viscosity is fitted by the Giesekus model, which, with the chosen parameters, leads to an extension thickening elongational viscosity. The stress and velocity fields of the numerical simulations are compared with the experimental results of Quinzani et al. [J. Non-Newtonian Fluid Mech. 52 (1994) 1–36] and the numerical results of the viscoelastic simulation using the Giesekus model of Azaiez et al. [J. Non-Newtonian Fluid Mech. 62 (1996) 253–277]. It can be shown that the quasi-Newtonian fluid qualitatively predicts the essential features of the flow in the vicinity of the contraction.     

Viscoelastic effects in double‐pipe single‐pass counterflow heat exchangers

The conducto‐convective heat loss from a viscoelastic liquid, in the core of a double‐pipe heat exchanger arrangement, to a cooler Newtonian fluid flowing in the outer annulus is investigated with direct numerical simulations. A numerical algorithm based on the finite difference method is implemented in time and space with the Giesekus constitutive model for the viscoelastic liquids. The flow of both the annulus and core‐fluids is considered to be Poiseuille flow, driven by respective pressure gradients. In general, the results show that a viscoelastic core‐fluid leads to slightly lower (albeit comparable) attainable temperatures in the core‐fluid stream as compared with a corresponding Newtonian fluid. Copyright © 2008 John Wiley & Sons, Ltd.     

Criteria of negative wake generation behind a cylinder

It was demonstrated by simulation in our previous study that both the normal stress and its gradient are responsible for the negative wake generation (overshoot in the axial velocity) and streamline shifting. Extensional properties of the fluids dominate the generation of the negative wake, while other factors strengthen or weaken the formation of velocity overshoot. In this study, the criteria for the negative wake generation are discussed in detail for various fluid models, including the PTT, the FENE-CR, the FENE-P, and the Giesekus models. With the FENE-CR fluid, it is easier to generate negative wake than with the FENE-P fluid. This confirms that the constant shear viscosity FENE-CR fluid enhances the velocity overshoot, and that the shear-thinning viscosity FENE-P fluid delays the negative wake generation. The Giesekus fluid has a similar behaviour to the PTT fluid with regarding to the critical conditions of negative wake generation when appropriate fluid parameters are selected. The mechanism of wall proximity in enhancing the negative wake generation is also demonstrated with the analysis for the first time.     

Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method

An efficient immersed boundary-lattice Boltzmann method (IB-LBM) is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows. Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor. A boundary thickening-based direct forcing IB method is adopted to solve the particle–fluid interactions with high accuracy for non-slip boundary conditions. A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computational parameters. An asynchronous calculation strategy is adopted to further improve the computing efficiency. The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods. The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed. Subsequently, 504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh–Taylor instabilities was observed. At last, simulations of a dense flow involving 11001 particles, the largest number of particles to date, were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid. The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time. These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.     

An extended finite element method for the simulation of particulate viscoelastic flows

We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.