A boundary element analysis of multiply connected three‐dimensional cavity mixing flow of polymer solutions

The emergence of non‐linear dynamics in cavity mixing is examined using the boundary element method (BEM). The method is implemented for the simulation of three‐dimensional transient creeping flow of Newtonian or linear viscoelastic fluids of the Jeffreys type. A boundary only formulation in the time domain is proposed for viscoelastic flow. Special emphasis is placed on cavity flow involving multiply connected moving domains. The BEM becomes particularly suited for this case, when part of the boundary (stirrer or rotor) is moving, and the remaining outer part (cavity) is at rest. In contrast to conventional volume methods, the BEM is shown to be much easier to implement since the kinematics of the elements bounding the fluid is known (imposed). It is found that, for a simple cavity flow induced by a rotating vane at constant angular velocity, the tractions at the vane tip and cavity face exhibit non‐linear periodic dynamical behaviour with time for fluids obeying linear constitutive equations. Copyright © 1999 John Wiley & Sons, Ltd.     

Three‐dimensional boundary element analysis of drop deformation in confined flow for Newtonian and viscoelastic systems

An adaptive (Lagrangian) boundary element approach is proposed for the general three‐dimensional drop deformation in confined flow. The adaptive method is stable as it includes remeshing capabilities of the deforming interface between drop and suspending fluid, and thus can handle large deformations. Both drop and surrounding fluid are viscous incompressible and can be Newtonian or viscoelastic. A boundary‐only formulation is implemented for fluids obeying the linear Jeffrey's constitutive equation. Similarly to the formulation for two‐dimensional Newtonian fluids (Khayat RE, Luciani A, Utracki LA. Boundary element analysis of planar drop deformation in confined flow. Part I. Newtonian fluids. Engineering Analysis of Boundary Elements 1997; 19 : 279), the method requires the solution of two simultaneous integral equations on the interface between the two fluids and the confining solid boundary. Although the problem is formulated for any confining geometry, the method is illustrated for a deforming drop as it is driven by the ambient flow inside a cylindrical tube. The accuracy of the method is assessed by comparison with the analytical solution for two‐phase radial spherical flow, leading to good agreement. The influence of mesh refinement is examined for a drop in simple shear flow. Copyright © 2000 John Wiley & Sons, Ltd.     

Boundary element analysis of three-dimensional mixing flow of Newtonian and viscoelastic fluids

The boundary element method (BEM) is implemented for the simulation of three-dimensional transient flows of typical relevance to mixing. Creeping Newtonian and viscoelastic fluids of the Maxwell type are examined. A boundary-only formulation in the time domain is proposed for linear viscoelastic flows. Special emphasis is placed on cavity flows involving simple- and multiple-connected moving domains. The BEM becomes particularly suited in multiple-connected flows, where part of the boundary (stirrer or rotor) is moving, and the remaining outer part (cavity or barrel) is at rest. In this case, conventional methods, such as the finite element method (FEM), generally require remeshing or mesh refinement of the three-dimensional fluid volume as the flow evolves and the domain of computation changes with time. The BEM is shown to be much easier to implement since the kinematics of the elements bounding the fluid is known (imposed). It is found that, for simple cavity flow induced by a rotating vane at constant angular velocity, the tractions at the vane tip and cavity face exhibit non-linear periodic dynamical behavior with time for fluids obeying linear constitutive equations. © 1998 John Wiley & Sons, Ltd.     

A 2D boundary element method for simulating the deformation of axisymmetric compound non‐Newtonian drops

The boundary integral formulation of the solution to the Stokes equations is used to describe the deformation of small compound non‐Newtonian axisymmetric drops suspended in a Newtonian fluid that is subjected to an axisymmetric flow field. The non‐Newtonian stress is treated as a source term in the Stokes equations, which yields an extra integral over the domains containing non‐Newtonian material. By transforming the integral representation for the velocity to cylindrical co‐ordinates and performing the integration over the azimuthal direction analytically, the dimension of the problem can be reduced from three to two. A boundary element method for the remaining two‐dimensional problem aimed at the simulation of the deformation of such axisymmetric compound non‐Newtonian drops is developed. Apart from a numerical validation of the method, simulation results for a drop consisting of an Oldroyd‐B fluid and a viscoelastic material are presented. Moreover, the method is extended to compound drops that are composed of a viscous inner core encapsulated by a viscoelastic material. The simulation results for these drops are verified against theoretical results from literature. Moreover, it is shown that the method can be used to identify the dominant break‐up mechanism of compound drops in relation to the specific non‐Newtonian character of the membrane. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

A boundary‐only approach to the deformation of a shear‐thinning drop in extensional Newtonian flow

The influence of shear thinning on drop deformation is examined through a numerical simulation. A two‐dimensional formulation within the scope of the boundary element method (BEM) is proposed for a drop driven by the ambient flow inside a channel of a general shape, with emphasis on a convergent–divergent channel. The drop is assumed to be shear thinning, obeying the Carreau–Bird model and the suspending fluid is Newtonian. The viscosity of the drop at any time is estimated on the basis of a rate‐of‐strain averaged over the region occupied by the drop. The viscosity thus changes from one time step to the next, and it is strongly influenced by drop deformation. It is found that small drops, flowing on the axis, elongate in the convergent part of the channel, then regain their spherical form in the divergent part; thus confirming experimental observations. Newtonian drops placed off‐axis are found to rotate during the flow with the period related to the initial extension, i.e. to the drop aspect ratio. This rotation is strongly prohibited by shear thinning. The formulation is validated by monitoring the local change of viscosity along the interface between the drop and the suspending fluid. It is found that the viscosity averaged over the drop compares, generally to within a few per cent, with the exact viscosity along the interface.     

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.     

Lattice‐Boltzmann and finite element simulations of fluid flow in a SMRX Static Mixer Reactor

A detailed comparison between the finite element method (FEM) and the lattice‐Boltzmann method (LBM) is presented. As a realistic test case, three‐dimensional fluid flow simulations in an SMRX static mixer were performed. The SMRX static mixer is a piece of equipment with excellent mixing performance and it is used as a highly efficient chemical reactor for viscous systems like polymers. The complex geometry of this mixer makes such three‐dimensional simulations non‐trivial. An excellent agreement between the results of the two simulation methods was found. Furthermore, the numerical results for the pressure drop as a function of the flow rate were close to experimental measurements. Results show that the relatively simple LBM is a good alternative to traditional methods. Copyright © 1999 John Wiley & Sons, Ltd.     

大规模有限元系统的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 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.     

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.     

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.