Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

A spectral approach to inertial confined thin‐film flow

A spectral approach is proposed to determine the flow field of a thin film inside narrow channels of arbitrary shape. Although the method is easily extended to transient flow, only steady flow is considered here. The flow field is represented spectrally in the depthwise direction in terms of orthonormal shape functions, which together with the Galerkin projection lead to a system of ordinary differential equations that can be solved using standard methods. The method is particularly effective for nonlinear flow, including nonlinearities of geometrical or material origins. The validity of the proposed method is demonstrated for a flow with inertia, and, unlike the depth‐averaging method, is not limited to a flow at small Reynolds number. The problem is closely related to high‐speed lubrication flow. The validity of the spectral representation is assessed by examining the convergence of the method, and comparing it with the fully two‐dimensional finite‐element solution, and the widely used depth‐averaging method from shallow‐water theory. It is found that a low number of modes are usually sufficient to secure convergence and accuracy. The influence of inertia is examined on the velocity and pressure fields. The pressure distributions reflect excellent agreement between the low‐order spectral method and the finite‐element solution, even at moderately high Reynolds number. The depth‐averaging solution is unable to predict accurately (qualitatively and quantitatively) the high‐inertia flow. Comparison of the velocity field reflects the expected discrepancy in a boundary layer formulation. Copyright © 2012 John Wiley & Sons, Ltd.     

On the use of characteristic‐based split meshfree method for solving flow problems

This study presents characteristic‐based split (CBS) algorithm in the meshfree context. This algorithm is the extension of general CBS method which was initially introduced in finite element framework. In this work, the general equations of flow have been represented in the meshfree context. A new finite element and MFree code is developed for solving flow problems. This computational code is capable of solving both time‐dependent and steady‐state flow problems. Numerical simulation of some known benchmark flow problems has been studied. Computational results of MFree method have been compared to those of finite element method. The results obtained have been verified by known numerical, analytical and experimental data in the literature. A number of shape functions are used for field variable interpolation. The performance of each interpolation method is discussed. It is concluded that the MFree method is more accurate than FEM if the same numbers of nodes are used for each solver. Meshfree CBS algorithm is completely stable even at high Reynolds numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

A coupled continuous and discontinuous finite element method for the incompressible flows

In this paper, we develop a coupled continuous Galerkin and discontinuous Galerkin finite element method based on a split scheme to solve the incompressible Navier–Stokes equations. In order to use the equal order interpolation functions for velocity and pressure, we decouple the original Navier–Stokes equations and obtain three distinct equations through the split method, which are nonlinear hyperbolic, elliptic, and Helmholtz equations, respectively. The hybrid method combines the merits of discontinuous Galerkin (DG) and finite element method (FEM). Therefore, DG is concerned to accomplish the spatial discretization of the nonlinear hyperbolic equation to avoid using the stabilization approaches that appeared in FEM. Moreover, FEM is utilized to deal with the Poisson and Helmholtz equations to reduce the computational cost compared with DG. As for the temporal discretization, a second‐order stiffly stable approach is employed. Several typical benchmarks, namely, the Poiseuille flow, the backward‐facing step flow, and the flow around the cylinder with a wide range of Reynolds numbers, are considered to demonstrate and validate the feasibility, accuracy, and efficiency of this coupled method. Copyright © 2016 John Wiley & Sons, Ltd.     

Retracted and replaced: A flow‐condition‐based interpolation finite element procedure for triangular grids

A flow‐condition‐based interpolation finite element scheme is presented for use of triangular grids in the solution of the incompressible Navier–Stokes equations. The method provides spatially isotropic discretizations for low and high Reynolds number flows. Various example solutions are given to illustrate the capabilities of the procedure. This article and been retracted and replaced. See retraction and replacement notice DOI: 10.1002/fld.1247 . Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of mass transfer effects on the performance of journal bearings using micropolar lubricant

In this paper the static and dynamic characteristics of a finite hydrodynamic journal bearing with micropolar lubricant are analyzed. The effects of mass transfer of solid additives and contaminants in the lubricant oil, on the bearing characteristics are considered in this study. A generalized form of Reynolds equation is derived from the fluid flow equations with the effects of mass transfer across the fluid film considered. The generalized Reynolds equation is solved using Galerkin's weighted-residual finite element method to obtain the fluid pressure distribution in the bearing. The various static and dynamic characteristics are subsequently obtained and presented.     

Parallel domain decomposition method for finite element approximation of 3D steady state non‐Newtonian fluids

We introduce a stabilized finite element method for the 3D non‐Newtonian Navier–Stokes equations and a parallel domain decomposition method for solving the sparse system of nonlinear equations arising from the discretization. Non‐Newtonian flow problems are, generally speaking, more challenging than Newtonian flows because the nonlinearities are not only in the convection term but also in the viscosity term, which depends on the shear rate. Many good iterative methods and preconditioning techniques that work well for the Newtonian flows do not work well for the non‐Newtonian flows. We employ a Galerkin/least squares finite element method, with stabilization parameters adjusted to count the non‐Newtonian effect, to discretize the equations, and the resulting highly nonlinear system of equations is solved by a Newton–Krylov–Schwarz algorithm. In this study, we apply the proposed method to some inelastic power‐law fluid flows through the eccentric annuli with inner cylinder rotation and investigate the robustness of the method with respect to some physical parameters, including the power‐law index and the Reynolds number ratios. We then report the superlinear speedup achieved by the domain decomposition algorithm on a computer with up to 512 processors. Copyright © 2015 John Wiley & Sons, Ltd.     

Anisotropic adaptive stabilized finite element solver for RANS models

Aerodynamic characteristics of various geometries are predicted using a finite element formulation coupled with several numerical techniques to ensure stability and accuracy of the method. First, an edge‐based error estimator and anisotropic mesh adaptation are used to detect automatically all flow features under the constraint of a fixed number of elements, thus controlling the computational cost. A variational multiscale‐stabilized finite element method is used to solve the incompressible Navier‐Stokes equations. Finally, the Spalart‐Allmaras turbulence model is solved using the streamline upwind Petrov‐Galerkin method. This paper is meant to show that the combination of anisotropic unsteady mesh adaptation with stabilized finite element methods provides an adequate framework for solving turbulent flows at high Reynolds numbers. The proposed method was validated on several test cases by confrontation with literature of both numerical and experimental results, in terms of accuracy on the prediction of the drag and lift coefficients as well as their evolution in time for unsteady cases.     

Numerical simulations of viscous flows using a meshless method

This paper uses the element‐free Galerkin (EFG) method to simulate 2D, viscous, incompressible flows. The control equations are discretized with the standard Galerkin method in space and a fractional step finite element scheme in time. Regular background cells are used for the quadrature. Several classical fluid mechanics problems were analyzed including flow in a pipe, flow past a step and flow in a driven cavity. The flow field computed with the EFG method compared well with those calculated using the finite element method (FEM) and finite difference method. The simulations show that although EFG is more expensive computationally than FEM, it is capable of dealing with cases where the nodes are poorly distributed or even overlap with each other; hence, it may be used to resolve remeshing problems in direct numerical simulations. Flows around a cylinder for different Reynolds numbers are also simulated to study the flow patterns for various conditions and the drag and lift forces exerted by the fluid on the cylinder. These forces are calculated by integrating the pressure and shear forces over the cylinder surface. The results show how the drag and lift forces oscillate for high Reynolds numbers. The calculated Strouhal number agrees well with previous results. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of flow in the plate‐spiral of a reaction turbine using a streamline upwind Petrov–Galerkin method

The prediction of the flow field in a novel spiral casing has been accomplished. Hydraulic turbine manufacturers are considering the potential of using a special type of spiral casing because of the easier manufacturing process involved in its fabrication. These special spiral casings are known as plate‐spirals. Numerical simulation of complex three‐dimensional flow through such spiral casings has been accomplished using a finite element method (FEM). An explicit Eulerian velocity correction scheme has been deployed to solve the Reynolds‐average Navier–Stokes equations. The simulation has been performed to describe the flow in high Reynolds number (10⁶) regimes. For spatial discretization, a streamline upwind Petrov–Galerkin (SUPG) technique has been used. The velocity field and the pressure distribution inside the spiral casing reveal meaningful results. Copyright © 2000 John Wiley & Sons, Ltd.     

Application of a preconditioned high‐order accurate artificial compressibility‐based incompressible flow solver in wide range of Reynolds numbers

In the present study, the preconditioned incompressible Navier‐Stokes equations with the artificial compressibility method formulated in the generalized curvilinear coordinates are numerically solved by using a high‐order compact finite‐difference scheme for accurately and efficiently computing the incompressible flows in a wide range of Reynolds numbers. A fourth‐order compact finite‐difference scheme is utilized to accurately discretize the spatial derivative terms of the governing equations, and the time integration is carried out based on the dual time‐stepping method. The capability of the proposed solution methodology for the computations of the steady and unsteady incompressible viscous flows from very low to high Reynolds numbers is investigated through the simulation of different 2‐dimensional benchmark problems, and the results obtained are compared with the existing analytical, numerical, and experimental data. A sensitivity analysis is also performed to evaluate the effects of the size of the computational domain and other numerical parameters on the accuracy and performance of the solution algorithm. The present solution procedure is also extended to 3 dimensions and applied for computing the incompressible flow over a sphere. Indications are that the application of the preconditioning in the solution algorithm together with the high‐order discretization method in the generalized curvilinear coordinates provides an accurate and robust solution method for simulating the incompressible flows over practical geometries in a wide range of Reynolds numbers including the creeping flows.     

A finite element method for the two‐dimensional extended Boussinesq equations

A new numerical method for Nwogu's (ASCE Journal of Waterway, Port, Coastal and Ocean Engineering 1993; 119 :618)two‐dimensional extended Boussinesq equations is presented using a linear triangular finite element spatial discretization coupled with a sophisticated adaptive time integration package. The authors have previously presented a finite element method for the one‐dimensional form of these equations (M. Walkley and M. Berzins (International Journal for Numerical Methods in Fluids 1999; 29 (2):143)) and this paper describes the extension of these ideas to the two‐dimensional equations and the application of the method to complex geometries using unstructured triangular grids. Computational results are presented for two standard test problems and a realistic harbour model. Copyright © 2002 John Wiley & Sons, Ltd.     

A two‐phase adaptive finite element method for solid–fluid coupling in complex geometries

In this paper we present a method to solve the Navier–Stokes equations in complex geometries, such as porous sands, using a finite‐element solver but without the complexity of meshing the porous space. The method is based on treating the solid boundaries as a second fluid and solving a set of equations similar to those used for multi‐fluid flow. When combined with anisotropic mesh adaptivity, it is possible to resolve complex geometries starting with an arbitrary coarse mesh. The approach is validated by comparing simulation results with available data in three test cases. In the first we simulate the flow past a cylinder. The second test case compares the pressure drop in flow through random packs of spheres with the Ergun equation. In the last case simulation results are compared with experimental data on the flow past a simplified vehicle model (Ahmed body) at high Reynolds number using large‐eddy simulation (LES). Results are in good agreement with all three reference models. Copyright © 2010 John Wiley & Sons, Ltd.     

Die swell measurements of second‐order fluids: numerical experiments

An analysis of the flow of a second‐order fluid is presented. Reference values for some variables are defined, and with these a non‐dimensional formulation of the governing equations. From this formulation, three dimensionless numbers appear; one is the Reynolds number, and two numbers that are called the first‐ and second‐dimensionless normal stress (NSD) coefficients. The equations of motion are solved by a finite element method using a commercially available program (Fidap), and the steady state converged solution was used to measure the die swell. The factors that influence die swell and that are studied in this work include: the die geometry for circular cross sectional dies, including tubular, converging, diverging, half‐converging/half‐tubular shapes; fluid characteristics such as Reynolds number and first‐ and second‐DNS coefficients (both positive and negative values); and flow rates, as determined by the maximum velocity in a parabolic velocity profile at the entrance to the die. The results suggest that shear and deformation histories of the fluid directly influence not only swell characteristics, but also convergence characteristics of the numerical simulation. © 1999 John Wiley & Sons, Ltd.     

A comparison of preconditioners for incompressible Navier–Stokes solvers

We consider solution methods for large systems of linear equations that arise from the finite element discretization of the incompressible Navier–Stokes equations. These systems are of the so‐called saddle point type, which means that there is a large block of zeros on the main diagonal. To solve these types of systems efficiently, several block preconditioners have been published. These types of preconditioners require adaptation of standard finite element packages. The alternative is to apply a standard ILU preconditioner in combination with a suitable renumbering of unknowns. We introduce a reordering technique for the degrees of freedom that makes the application of ILU relatively fast. We compare the performance of this technique with some block preconditioners. The performance appears to depend on grid size, Reynolds number and quality of the mesh. For medium‐sized problems, which are of practical interest, we show that the reordering technique is competitive with the block preconditioners. Its simple implementation makes it worthwhile to implement it in the standard finite element method software. Copyright © 2007 John Wiley & Sons, Ltd.     

A comparative study of GLS finite elements with velocity and pressure equally interpolated for solving incompressible viscous flows

A comparative study of the bi‐linear and bi‐quadratic quadrilateral elements and the quadratic triangular element for solving incompressible viscous flows is presented. These elements make use of the stabilized finite element formulation of the Galerkin/least‐squares method to simulate the flows, with the pressure and velocity fields interpolated with equal orders. The tangent matrices are explicitly derived and the Newton–Raphson algorithm is employed to solve the resulting nonlinear equations. The numerical solutions of the classical lid‐driven cavity flow problem are obtained for Reynolds numbers between 1000 and 20 000 and the accuracy and converging rate of the different elements are compared. The influence on the numerical solution of the least square of incompressible condition is also studied. The numerical example shows that the quadratic triangular element exhibits a better compromise between accuracy and converging rate than the other two elements. Copyright © 2008 John Wiley & Sons, Ltd.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

Solving the fully nonlinear weakly dispersive Serre equations for flows over dry beds

We describe a numerical method for solving the Serre equations that can simulate flows over dry bathymetry. The method solves the Serre equations in conservation law form with a finite volume method. A finite element method is used to solve the auxiliary elliptic equation for the depth‐averaged horizontal velocity. The numerical method is validated against the lake at rest analytic solution, demonstrating that it is well‐balanced. Since there are currently no known nonstationary analytical solutions to the Serre equation that involve bathymetry, a nonstationary forced solution, involving bathymetry was developed. The method was further validated and its convergence rate established using the developed nonstationary forced solution containing the wetting and drying of bathymetry. Finally, the method is also validated against experimental results for the run‐up of a solitary wave on a sloped beach. The finite‐volume finite‐element approach to solving the Serre equation was found to be accurate and robust.     

Finite element solution of three‐dimensional turbulent flows applied to mold‐filling problems

This paper presents a finite element solution algorithm for three‐dimensional isothermal turbulent flows for mold‐filling applications. The problems of interest present unusual challenges for both the physical modelling and the solution algorithm. High‐Reynolds number transient turbulent flows with free surfaces have to be computed on complex three‐dimensional geometries. In this work, a segregated algorithm is used to solve the Navier–Stokes, turbulence and front‐tracking equations. The streamline–upwind/Petrov–Galerkin method is used to obtain stable solutions to convection‐dominated problems. Turbulence is modelled using either a one‐equation turbulence model or the κ–ε two‐equation model with wall functions. Turbulence equations are solved for the natural logarithm of the turbulence variables. The change of dependent variables allows for a robust solution algorithm and good predictions even on coarse meshes. This is very important in the case of large three‐dimensional applications for which highly refined meshes result in untreatable large numbers of elements. The position of the flow front in the mold cavity is computed using a level set approach. Finally, equations are integrated in time using an implicit Euler scheme. The methodology presents the robustness and cost effectiveness needed to tackle complex industrial applications. Copyright © 2000 John Wiley & Sons, Ltd.