Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

Bubble collapse in compressible fluids using a spectral element marker particle method. Part 2. Viscoelastic fluids

This paper is concerned with the development of a high‐order numerical scheme for two‐phase viscoelastic flows. In the companion paper, herein referred to as Part 1, the scheme is applied to the modelling of two‐phase Newtonian flows. The particular problem of the collapse of a 2D bubble in the vicinity of a rigid boundary is considered. Attention is given to the construction of the most general form of the compressible Oldroyd B model that is consistent with the compressible Newtonian and upper‐convected Maxwell models in the appropriate limits. The governing equations are discretized using the spectral element method, and the two phases are modelled using a marker particle method. A comprehensive set of results is presented for the problem of bubble collapse near a rigid wall, and qualitative agreement is obtained with other numerical studies and experimental observations. Viscoelastic effects that are predicted include increased bubble oscillation with increasing Weissenberg number and considerable bubble deformation and cusping near the wall. Most importantly, it has been shown that viscoelasticity has the ability to prevent jet formation and therefore is likely to have a mitigating effect on cavitation damage. Copyright © 2012 John Wiley & Sons, Ltd.     

Time‐dependent finite element simulations of a shear‐thinning viscoelastic fluid with application to blood flow

A new finite element method is developed to simulate time‐dependent viscoelastic shear‐thinning flows characterized by the generalized Oldroyd‐B model. The focus of the algorithm is improved stability through a free‐energy dissipative scheme by using low‐order piecewise‐constant finite element approximations for stress. The algorithm is further modified by incorporating a pressure‐projection method, a DG‐upwinding scheme, a symmetric interior penalty DG method to solve the elliptic pressure‐update equation and a geometric multigrid preconditioner. The improved stability and cost to accuracy is compared when using higher order discontinuous bilinear approximation, where in addition, we consider the influence of a slope limiter for these elements. The algorithm is applied to the 2D start‐up‐driven cavity problem, and the stability of the free energy is illustrated and compared between element choices. An application of the model to modelling blood in small arterioles and channels is considered by simulating pulsatile blood flow through a stenotic arteriole. The individual influences of viscoelasticity and shear‐thinning within the generalized Oldroyd‐B model are investigated by comparing results to the Newtonian, generalized Newtonian and Oldroyd‐B models. Copyright © 2014 John Wiley & Sons, Ltd.     

Development and implementation of VOF-PROST for 3D viscoelastic liquid–liquid simulations

We implement a volume-of-fluid algorithm with a parabolic re-construction of the interface for the calculation of the surface tension force (VOF-PROST). This achieves higher accuracy for drop deformation simulations in comparison with existing VOF methods based on a piecewise linear interface re-construction. The algorithm is formulated for the Giesekus constitutive law. The evolution of a drop suspended in a second liquid and undergoing simple shear is simulated. Numerical results are first checked against two cases in the literature: the small deformation theory for second-order liquids, and an Oldroyd-B extensional flow simulation. We then address the experimental data of Guido et al. (2003) for a Newtonian drop in a viscoelastic matrix liquid. The data deviate from existing theories as the capillary number increases, and reasons for this are explored here with the Oldroyd-B and Giesekus models.     

大规模有限元系统的GPU加速计算研究

研究了GPU(Graphics Processing Units)计算应用于有限元方法中的总刚计算和组装、稀疏矩阵与向量乘积运算、线性方程组求解问题,并基于CUDA(Compute Unified Device Architecture)平台利用GTX295GPU进行程序实现和测试。系统总刚采用CSR(Compressed Sparse Row)压缩格式存放于GPU显存中,用单元染色方法实现总刚并行计算组装,用共轭梯度迭代法求解大规模线性方程组。对300万自由度以内的空间桁架和平面问题算例,GPU有限元计算分别获得最高9.5倍和6.5倍的计算加速比,并且加速比随系统自由度的增加而近似线性增加,GFLOP/s峰值也有近10倍的增加。     

Effects of viscoelasticity on the retraction of a sheared drop

Effects of drop and matrix viscoelasticity on the retraction of a sheared drop are numerically investigated. Retraction of an Oldroyd-B drop in a Newtonian matrix is initially faster and later slower with increasing drop Deborah number. The observed behavior is explained using an ordinary differential equation model representing the dominant balance between various forces during retraction. The initial faster relaxation of viscoelastic drops is due to viscoelastic stresses pulling the drop interface at the tips inward. The later slower retraction is due to the slowly-relaxing viscoelastic forces at the equator, where they act against the capillary force. The drop inclination decreases substantially during retraction unlike in a Newtonian case. Matrix viscoelasticity slows the relaxation of a Newtonian drop because of the increasingly slow relaxation of highly stretched polymers near the drop tip with increasing Deborah number. Increasing the ratio of polymeric to total viscosity further accentuates the viscoelastic effects in both cases. For an Oldroyd-B drop in an Oldroyd-B matrix, a competition between the dispersed and the continuous phase elasticities, represented by their ratio, determines the dynamics; larger values of the ratio leads again to initial faster and later slower retraction.     

Treatment of the small time instability in the finite element analysis of fluid structure interaction problems

In this paper, the fluid–structure interaction problem in mechanical systems in which a high frequency vibrating solid structure interacts with the surrounding fluid flow is considered. Such a situation normally appears in many microelectromechanical systems like a wide variety of microfluidic devices. A different implementation of the residual‐based variational multiscale flow method is employed within the arbitrary Lagrangian–Eulerian formulation. The combination of the variational multiscale method with appropriate stabilization parameters is used to handle the so‐called small time step instability in the finite element analysis of the fluid part in the coupled fluid–structure interaction problem. The capability of the employed approach has been demonstrated through finite element study of a benchmark example and FEM simulation of two different mechanical micropumping devices. High frequency vibrations of the solid membrane are used to derive the fluid flow in these micropumps. Results of FEM simulations are shown to be in good agreement with available experimental data.Copyright © 2012 John Wiley & Sons, Ltd.     

Spectral element method for viscoelastic flows in a planar contraction channel

A new algorithm, which combines the spectral element method with elastic viscous splitting stress (EVSS) method, has been developed for viscoelastic fluid flows in a planar contraction channel. The system of spectral element approximations to the velocity, pressure, extra stress and the rate of deformation variables is solved by a preconditioned conjugate gradient method based on the Uzawa iteration procedure. The numerical approach is implemented on a planar four‐to‐one contraction channel for a fluid governed by an Oldroyd‐B constitutive equation. The behaviour of the Oldroyd‐B fluids in the contraction channel is investigated with various Weissenberg numbers. It is shown that numerical solutions obtained here agree well with experimental measurements and other numerical predictions. Copyright © 2003 John Wiley & Sons, Ltd.     

Rheology of an emulsion of viscoelastic drops in steady shear

Steady shear rheology of a dilute emulsion with viscoelastic inclusions is numerically investigated using direct numerical simulations. Batchelor's formulation for rheology of a viscous emulsion is extended for a viscoelastic system. Viscoelasticity is modeled using the Oldroyd-B constitutive equation. A front-tracking finite difference code is used to numerically determine the drop shape, and solve for the velocity and stress fields. The effective stress of the viscoelastic emulsion has three different components due to interfacial tension, viscosity difference (not considered here) and the drop phase viscoelasticity. The interfacial contributions – first and second normal stress differences and shear stresses – vary with Capillary number in a manner similar to those of a Newtonian system. However the shear viscosity decreases with viscoelasticity at low Capillary numbers, and increases at high Capillary numbers. The first normal stress difference due to interfacial contribution decreases with increasing drop phase viscoelasticity. The first normal stress difference due to the drop phase viscoelasticity is found to have a complex dependence on Capillary and Deborah numbers, in contrast with the linear mixing rule. Drop phase viscoelasticity does not contribute significantly to effective shear viscosity of the emulsion. The total first normal stress difference shows an increase with drop phase viscoelasticity at high Capillary numbers. However at low Capillary numbers, a non-monotonic behavior is observed. The results are explained by examining the stress field and the drop shape.     

Concentration dependence of the linear viscoelastic properties of particle suspensions

The response under small amplitude oscillatory deformations of a suspension of non-Brownian spheres dispersed in a viscoelastic fluid is investigated. The correspondence principle of linear viscoelasticity is used to derive a simple constitutive model from a model for a suspension in a Newtonian liquid. The theory predicts that for a specific particulate system the concentration dependence of the viscoelastic properties should collapse to a single master curve when the values are normalized with those of the carrier fluid alone. Measurements with the micro-Fourier rheometer using oscillatory squeeze flow are carried out on two suspensions of 60 and 80 μm sized particles dispersed in polymeric fluid and in silicon oil, and the master curve is verified. Received: 27 April 1999/Accepted: 15 October 1999     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

A DOMAIN DECOMPOSITION ALGORITHM WITH FINITE ELEMENT-BOUNDARY ELEMENT COUPLING

A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.     

Boundary element simulations of spheres settling in circular,square and triangular conduits

Numerical simulations of a spherical particle sedimenting in circular, triangular and square conduits containing a viscous, inertialess, Newtonian fluid were investigated using the Boundary Element Method (BEM). Settling velocities and pressure drops for spheres falling along the centre-lines of the conduits were computed for a definitive range of sphere sizes. The numerical simulations for the settling velocities showed good agreement with existing experimental data. The most accurate analytic solution for a sphere settling along the axis of a circular conduit produced results which were almost indistinguishable from the present BEM calculations. For a sphere falling along the centre-line of a square conduit, the BEM calculations for small spheres agreed well with analytic results. No analytic results for a sphere falling along the axis of a triangular conduit were available for comparison. Extrapolation of the BEM predictions for the pressure drops, to infinitely small spheres, showed remarkable agreement with analytic results. For the circular conduit, the sphere's settling velocity and angular velocity were computed as a function of drop position for small, medium and large spheres. Excellent agreement with a reflection solution was achieved for the small sphere. In addition, end effects were investigated for centre-line drops and compared where possible with available experimental data and analytic results.Los Alamos National Laboratory, Los Alamos, New Mexico, USA.     

Instability of stationary and uniformly moving cylindrical fluid bodies—II. Viscoelastic threads and experimental observations

The instability analysis of Part I is extended to the breakup of viscoelastic threads in fluid media (also possibly viscoelastic). Critical Growth rates and wave-numbers are calculated in terms of the viscosity ratio, the Ohnesorge numbers (continuous and dispersed phases), and elasticity numbers for each of the respective phases. Comparisons with results for Newtonian systems indicate viscoelastic threads to be less stable than Newtonian threads under similar conditions. Also, the critical wave-numbers observed with viscoelastic threads can differ significantly from those observed with Newtonian systems, particularly if the relative magnitudes of elasticity of the dispersed and continuous phases are quite different. Systems with similar magnitudes of elasticity in each phase exhibit wave-numbers similar to Newtonian systems of similar viscosities.Experimental results obtained from observations of fluid thread breakup in a Taylor four-roller device provide a basis for checking the predictions of the lineararized theory for both Newtonian and viscoelastic systems. In general, the agreement is good and the theoretical predictions of Parts I and II seem to be reasonable representations of experimental fact.     

Drop shape under slow steady shear flow and during relaxation. Experimental results and comparison with theory

The three-dimensional small deformation of a single Newtonian drop immersed in an immiscible Newtonian liquid was investigated, both in slow steady shear and during retraction after cessation of shear. The experiments were performed in a parallel plate apparatus equipped with video- enhanced microscopy. The drop was observed from two perpendicular directions, and accurate measurements were obtained in each view by image analysis. The results were compared to existing theoretical predictions from the exact fluid dynamic problem, obtained perturbatively for the case of small deformations of the drop. Excellent agreement between data and theory was found, thus providing the first assessment of drop shape predictions in steady flow and relaxation for Newtonian fluids. Received: 6 February 2000 Accepted: 24 May 2000     

可压缩流场中气泡脉动数值模拟

在应用边界元方法对气泡动力学的研究中, 绝大多数模型是建立在不压缩势流理论基础之上, 针对可压缩流场中气泡运动特性的研究很少. 从波动方程出发, 分别在气泡运动前期和后期对波动方程进行简化, 得到气泡运动局部和全局简化方程, 采用双渐进方法对简化方程进行匹配, 提出了考虑流场可压缩性的非球状气泡运动模型. 该模型的计算结果与Prospertti 等的解析结果吻合很好, 气泡脉动最大半径和内部最大压力随气泡脉动逐渐减小. 基于该模型对比了自由场中药包爆炸考虑可压缩性与不考虑可压缩性的计算结果, 发现考虑可压缩性气泡射流速度较小, 随后基于该模型计算了刚性边界下气泡的运动特性.     

Front‐tracking by the level‐set and the volume penalization methods in a two‐phase microfluidic network

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.     

Experimental investigation of viscoelastic drop deformation in Newtonian matrix at high capillary number under simple shear flow

The steady state deformation of a viscoelastic drop (Boger fluid) in a Newtonian liquid at high capillary number under simple shear flow is investigated by direct visualization using a specially designed Couette apparatus which enables visualization from two perpendicular directions. Two drop deformation modes are found: (1) Mode I – drop deformation in the flow direction and (2) Mode II – drop deformation in the vorticity direction. The drop deformation mode depends on the relative strength of the elastic contribution to viscous contribution. If the elastic contribution is weak compared to the viscous contribution, the drop elongates in the flow direction via Mode I. If the elastic contribution is strong, the drop elongates in the vorticity direction via Mode II. The drop size also affects the drop deformation. At the same capillary number, bigger drops have larger deformations than smaller drops.     

Bubble collapse in compressible fluids using a spectral element marker particle method. Part 1. Newtonian fluids

This paper is concerned with the development of a high‐order numerical scheme for the modelling of two‐phase Newtonian flows. The companion paper, herein referred to as Part 2, extends the scheme to two‐phase viscoelastic flows. The particular problem of the collapse of a two‐dimensional bubble in the vicinity of a rigid boundary is considered. The governing equations are discretized using the spectral element method, and the two phases are modelled using a marker particle method. The marker particle scheme is validated using the Zalesak slotted disk rotation test problem. A comprehensive set of results is presented for the problem of bubble collapse near a rigid wall, and qualitative agreement is obtained with other numerical studies and experimental observations. Viscous effects are shown to inhibit bubble collapse and prevent jet formation and are therefore likely to have a mitigating effect on cavitation damage.Copyright © 2011 John Wiley & Sons, Ltd.