Simulation of planar converging flow of a Leonov viscoelastic fluid

An efficient finite element algorithm is presented to simulate the planar converging flow for the viscoelastic fluid of the Leonov model. The governing equation set, composed of the continuity, momentum and constitutive equations for the Leonov fluid flow, is conveniently decoupled and a two-stage cyclic iteration technique is employed to solve the velocity and elastic strain fields separately. Artificial viscosity terms are imposed on the momentum equations to relax the elastic force and data smoothing is performed on the iterative calculations for velocities to further stabilize the numerical computations. The calculated stresses agree qualitatively with the experimental measurements and other numerically simulated results available in the literature. Computations were successful to moderately high values of Deborah number of about 27·5.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.     

Numerical study around the corotational Maxwell model for the viscoelastic fluid flows

We present numerical results for the FEM (finite element method) presented in [Comput. Methods Appl. Mech. Engrg. 191 (2002) 5045–5065]. This method is devoted to the approximation of fluid flows obeying the Oldroyd model. A particularity of this method, is to take into account the purely viscoelastic case, the so-called Maxwell model, important in practice. Numerical results are given for a fluid flowing in an abrupt plane 4 to 1 contraction. We use the corotational Maxwell model as benchmark in the choice of our computations. Results are also given for the upper convected Maxwell model. Interesting effects appear on the velocity profile: a phenomenon of quasi slip at the downstream wall.     

A least squares finite element method for viscoelastic fluid flow problems

Viscoelastic flows remain a demanding class of problems for approximate analysis, particularly at increasing Weissenberg numbers. Part of the difficulty stems from the convective behavior and in the treatment of the stress field as a primary unknown. This latter aspect has led to the use of higher-order piecewise approximations for the stress approximation spaces in recent finite element research. The computational complexity of the discretized problem is increased significantly by this approach but at present it appears the most viable technique for solving these problems. Motivated by recent success in treating mixed systems and convective problems, we formulate here a least squares finite element method for the viscoelastic flow problem. Numerical experiments are conducted to test the method and examine its strengths and limitations. Some difficulties and open issues are identified through the numerical experiments. We consider the use of high degree elements (p refinement) to improve performance and accuracy.     

Using a segregated finite element scheme to solve the incompressible Navier-Stokes equations

In this paper, a segregated finite element scheme for the solution of the incompressible Navier-Stokes equations is proposed which is simpler in form than previously reported formulations. A pressure correction equation is derived from the momentum and continuity equations, and equal-order interpolation is used for both the velocity components and pressure. Algorithms such as this have been known to lead to checkerboard pressure oscillations; however, the pressure correction equation of this scheme should not produce these oscillations. The method is applied to several laminar flow situations, and details of the methods used to achieve converged solutions are given.     

Spectral-finite element method for compressible fluid flow

In this paper, a combined Fourier spectral-finite element method is proposed for solving n-dimensional (n=2, 3), semi-periodio compressible fiuid flow problems. The strict error estimation as well as the convergence rate, is presented.     

Finite element analysis of convective heat transfer in porous media

The finite element method is used to analyse convective heat transfer in a porous medium. Convection past a vertical surface embedded in the medium and convection in a confined porous medium enclosure are analysed using the above method. The results are compared with those available in the literature and the agreement is found to be good. The method is applicable for two-dimensional analysis in a porous body of any arbitrary shape. The restriction of the boundary layer assumption is relaxed.     

A new decoupled finite element algorithm for viscoelastic flow. Part 2: Convergence properties of the algorithm

We present the results of some numerical experiments which were carried out in order to investigate the general characteristics of the algorithm described in Part I of this paper.     

Finite element analysis of the axially symmetric motion of an incompressible viscous fluid in a spherical annulus

This paper presents results obtained by employing a modified Galerkin finite element method to analyse the steady state flow of a fluid contained between two concentric, rotating spheres. The spheres are assumed to be rigid and the cavity region between the spheres is filled with an incompressible, viscous, Newtonian fluid. The inner sphere is constrained to rotate about a vertical axis with a prescribed angular velocity, while the outer sphere is fixed. Results for the circumferential function Ω, streamfunction ψ, vorticity function ζ and inner boundary torque T₁ are presented for Reynolds numbers Re ? 2000 and radius ratios 0.1 ? α ? 0.9. The method proved effective for obtaining results for a wide range of radius ratios (0.1 ? α ? 0.9) and Reynolds numbers (0 ? Re ? 2000). Previous investigators who employed the finite difference method experienced difficulties in obtaining results for cases with radius ratios α ? 0.2, except for small Reynolds numbers (Re ? 100). Results for Ω, Ψ, ζ and T₁ obtained in this study for radius ratios 0.8 ≤ α ≤ 0.9 verified the development of Taylor vortices reported by other investigators. The research indicates that the method may be useful for analysing other non-linear fluid flow problems.     

A finite element formulation for laminar flow of a fluid with microstructure

A finite element formulation for the steady laminar flow of an incompressible fluid with microstructure has been developed. The particular fluids considered are commonly known as micropolar fluids, in which case suspended particulate microstructures are modelled by an ‘extended’ continuum formulation. The particle microspin is a new kinematic variable which is independent of the classical vorticity vector and thereby allows relative rotation between particles and the surrounding fluid. This formulation also gives rise to couple stresses in addition to classical force or traction stresses. The finite element formulation utilizes a variational approach and imposes conservation of mass through a penalty function. A general boundary condition for microspin has been incorporated whereby microspin at a solid boundary is constrained to be proportional to the fluid vorticity. The proportionality constant in this case can vary from zero to unity. Sample solutions are presented for fully developed flow through a straight tube and compared with an analytical solution. Results are also generated for flow through a constricted tube and compared with a Newtonian fluid solution.     

Finite element modelling of two-phase heat and fluid flow in deforming porous media

This paper details a finite element model which describes the flow of two-phase fluid and heat within a deforming porous medium. The coupled governing equations are derived in terms of displacements, pore pressures and temperatures, and details of the time-stepping algorithm and thermodynamic considerations are also presented. Two numerical examples are included for verification.     

An adaptive remeshing strategy for viscoelastic fluid flow simulations

In the last few years, we have developed a fairly general adaptive finite element solution procedure which can be applied to a large variety of problems. In this paper, this strategy is briefly recalled and applied to the solution of two-dimensional viscoelastic fluid flow problems. A log-conformation formulation recently introduced by Fattal and Kupferman [R. Fattal, R. Kupferman, Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation, J. Non-Newtonian Fluid Mech. 126 (2005) 23-37] was implemented in order to improve the convergence properties of the numerical scheme. We confirm some results obtained in Hulsen, Fattal and Kupferman [M. Hulsen, R. Fattal, R. Kupferman, Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithm, J. Non-Newtonian Fluid Mech. 127 (2005) 27-39] and in some instances, we show that mesh adaptation allows to almost automatically reproduce accurate results obtained on very fine structured meshes.     

A new decoupled finite element algorithm for viscoelastic flow. Part 1: Numerical algorithm and sample results

This paper presents an algorithm for two-dimensional Steady viscoelastic flow Simulation in which the Solution of the momentum and continuity equations is decoupled from that of the constitutive equations. The governing equations are discretized by the finite element method, with 3 × 3 element subdivision for the stress field approximation. Non-consistent Streamline upwinding is also used. Results are given for flow through a converging channel and through an abrupt planar 4:1 contraction.     

三晶体中三晶交线附近弹性应变能集中的有限单元分析

提出应变能集中这一概念,对三晶体中三晶交线附近的弹性应变能的集中特性进行了计算,发现当加载方向平行于三晶交线时,三晶交线附近有弹性应变能集中,但是集中程度不是很高,当加载方向垂直于三晶交线时,根据三晶体中组元晶体的相对位置的不同和取向的不同,三晶交线附近或者没有弹性应变能集中,或者有弹性应变能集中,而且在一些情况下与平行于三晶交线方向加载的情况相比,弹性应变能的集中程度要高很多。     

Numerical simulation method for viscoelastic flows with free surfaces—fringe element generation method

To simulate filling flow in injection moulding for viscoelastic fluids, a numerical method, based on a finite element method and a finite volume method, has been developed for incompressible isothermal viscoelastic flow with moving free surfaces. The advantages of this method are, first, good applicability to arbitrarily shaped mould geometries and, second, accurate treatment for boundary conditions on the free surface. Typical filling flows are simulated, namely filling flow into a 1:4 expansion cavity with and without an obstacle. Numerical results predict the position of weld lines and air-traps. The method also indicates the effects of elongational flow on molecular orientation.     

Finite element computations for unsteady fluid and elastic membrane interaction problems

The interaction between the hydrodynamic forces of a flow field and the elastic forces of adjacent deformable boundaries is described by elastohydrodynamics, a coupled fluid–elastic membrane problem. Direct numerical solution of the unsteady, highly non-linear equations requires that the dynamic evolution of both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. This paper describes a numerical algorithm based on the deformable spatial domain space–time (DSD/ST) finite element method for the unsteady motion of an incompressible, viscous fluid with elastic membrane interaction. The unsteady Navier–Stoke and elastic membrane equations are solved separately using an iterative procedure by the GMRES technique with an incomplete lower-upper (ILU) decomposition at every time instant. One-dimensional, two-dimensional and deformable domain model problems are used to demonstrate the capabilities and accuracy of the present algorithm. Both steady state and transient problems are studied. © 1997 John Wiley & Sons, Ltd.     

An orthogonal mapping technique for the computation of a viscous free-surface flow

In this paper we describe finite element computations of the free-surface flow of a viscous fluid down an undulating inclined plane. The technique developed here employs an orthogonal mapping that is computed along with the velocity and pressure. This is allied to a technique to compute symbolically the Jacobian and other derivatives required for numerical continuation methods. The solutions obtained are compared with laboratory experiments and finite element computations reported by Pritchard and co-workers. The finite element computational method used by these authors employs spines to represent the free surface. An excellent agreement is shown to exist between the new computations and the laboratory experiments, and with the numerical solutions of Pritchard and co-workers.     

A total linearization method for solving viscous free boundary flow problems by the finite element method

In this paper a total linearization method is derived for solving steady viscous free boundary flow problems (including capillary effects) by the finite element method. It is shown that the influence of the geometrical unknown in the totally linearized weak formulation can be expressed in terms of boundary integrals. This means that the implementation of the method is simple. Numerical experiments show that the iterative method gives accurate results and converges very fast.     

On the modelling of highly elastic flows of amorphous thermoplastics

Two approaches to the kinematic structuring of constitutive models for highly elastic flows of polymer melts have been examined systematically, assuming either: (1) additivity of elastic and viscous velocity gradients or (2) multiplicability of elastic and viscous deformation gradients. A series of constitutive models were compared, with differing kinematic structure but the same linear responses in elastic and viscous limits. They were solved numerically and their predictions compared, and they were also compared to those of the Giesekus model. Several variants, previously proposed as separate models, are shown to be equivalent and qualitatively in agreement with experiment, and therefore a sound basis for construction of models. But the assignment of viscous spin is critical: if it is assumed equal to the total spin with approach (1), or equal to zero with approach (2), then unphysical viscoelastic behaviour is predicted.     

Finite element analysis of lift build-up due to aerofoil indicial motion

In this paper the problem of impulsively started aerofoil or suden change of incidence of an aerofoil in incompressible potential flow is investigated. The essence of solution lies in the representation of a timely and spatially varying wake in a largely irrotational potential flow field. This is achieved by representing the wake through velocity potential difference, which seems to be the only way of imposing a velocity difference condition in the finite element context with velocity potentials as the basic unknowns. Superposition is employed to meet various boundary conditions, which is justified by the linearity of the problem. The finite element solutions are compared with those from singularity method.