Three‐dimensional incompressible flow calculations using the characteristic based split (CBS) scheme

In this paper, the characteristic based split scheme is employed for the solution of three‐dimensional incompressible viscous flow problems on unstructured meshes. Many algorithm related issues are discussed. Fully explicit and semiimplicit forms of the scheme are explained and employed in the calculation of both isothermal and nonisothermal incompressible flows simulation. The extension of the scheme to porous medium flows is also demonstrated with relevant examples. Copyright © 2004 John Wiley & Sons, Ltd.     

An arbitrary Lagrangian Eulerian (ALE) formulation for free surface flows using the characteristic‐based split (CBS) scheme

An arbitrary Lagrangian Eulerian (ALE) method for non‐breaking free surface flow problems is presented. The characteristic‐based split (CBS) scheme has been employed to solve the ALE equations. A simple mesh smoothing procedure based on coordinate averaging (Laplacian smoothing) is employed in the calculations. The mesh velocity is calculated at each time step and incorporated as part of the scheme. Results presented show an excellent agreement with the available experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

The characteristic‐based split (CBS) meshfree method for free surface flow problems in ALE formulation

A semi‐implicit characteristic‐based split (CBS) meshfree algorithm in the arbitrary Lagrangian Eulerian (ALE) framework is proposed for the numerical solution of incompressible free surface flow problem in the paper. The algorithm is the extension of general CBS method which was initially introduced in finite element framework, this is due to the fact that CBS method not only can enhance the stability, but also avoid LBB condition when equal order basis function is used to approximate velocity and pressure variables. Meanwhile, a simple way for node update and node speed calculation is developed which is used to capture the free surface exactly. The numerical solutions are compared with available analytical and numerical solutions, which shows that the proposed method has better ability to simulate the free surface incompressible flow problem. Copyright © 2009 John Wiley & Sons, Ltd.     

An improvement for multidimensional characteristic‐based scheme by using different selected waves

In this study, an improved multidimensional characteristic‐based (MCB) scheme, used for calculation of incompressible flows via artificial compressibility method, is presented. The new idea uses four different waves instead of traditional waves, which were used in the original MCB scheme corresponding to coordinate axes. In the new selection of waves, the number of compatibility relations is increased to eight, and therefore, more information reaches to each face in order to estimate convective fluxes. Numerical tests show that the new improved MCB scheme is faster in terms of convergence and demonstrates the same accuracy by the older one. Also, it needs less CPU time in comparison with traditional MCB scheme for running the computer code at the same conditions. The new scheme is able to use higher CFL numbers compared with the old MCB scheme. Results obtained by the new proposed method are in good agreement with the benchmark solutions at the wide range of Reynolds numbers. Copyright © 2014 John Wiley & Sons, Ltd.     

Genuinely multidimensional characteristic‐based scheme for incompressible flows

This paper presents a novel multidimensional characteristic‐based (MCB) upwind method for the solution of incompressible Navier–Stokes equations. As opposed to the conventional characteristic‐based (CB) schemes, it is genuinely multidimensional in that the local characteristic paths, along which information is propagated, are used. For the first time, the multidimensional characteristic structure of incompressible flows modified by artificial compressibility is extracted and used to construct an inherent multidimensional upwind scheme. The new proposed MCB scheme in conjunction with the finite‐volume discretization is employed to model the convective fluxes. Using this formulation, the steady two‐dimensional incompressible flow in a lid‐driven cavity is solved for a wide range of Reynolds numbers. It was found that the new proposed scheme presents more accurate results than the conventional CB scheme in both their first‐ and second‐order counterparts in the case of cavity flow. Also, results obtained with second‐order MCB scheme in some cases are more accurate than the central scheme that in turn provides exact second‐order discretization in this grid. With this inherent upwinding technique for evaluating convective fluxes at cell interfaces, no artificial viscosity is required even at high Reynolds numbers. Another remarkable advantage of MCB scheme lies in its faster convergence rate with respect to the CB scheme that is found to exhibit substantial delays in convergence reported in the literature. The results obtained using new proposed scheme are in good agreement with the standard benchmark solutions in the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Implementing a high‐order accurate implicit operator scheme for solving steady incompressible viscous flows using artificial compressibility method

This paper uses a fourth‐order compact finite‐difference scheme for solving steady incompressible flows. The high‐order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two‐dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier–Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth‐order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block‐tridiagonal matrix inversions. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. In this study, the high‐order compact implicit operator scheme is also extended for computing three‐dimensional incompressible flows. The accuracy and efficiency of this high‐order compact method are demonstrated for different incompressible flow problems. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated. Test cases considered herein for validating the results are incompressible flows in a 2‐D backward facing step, a 2‐D cavity and a 3‐D cavity at different flow conditions. Results obtained for these cases are in good agreement with the available numerical and experimental results. The study shows that the scheme is robust, efficient and accurate for solving incompressible flow problems. Copyright © 2010 John Wiley & Sons, Ltd.     

A two‐step Taylor‐characteristic‐based Galerkin method for incompressible flows and its application to flow over triangular cylinder with different incidence angles

An alternative characteristic‐based scheme, the two‐step Taylor‐characteristic‐based Galerkin method is developed based on the introduction of multi‐step temporal Taylor series expansion up to second order along the characteristic of the momentum equation. Contrary to the classical characteristic‐based split (CBS) method, the current characteristic‐based method does not require splitting the momentum equation, and segregate the calculation of the pressure from that of the velocity by using the momentum–pressure Poisson equation method. Some benchmark problems are used to examine the effectiveness of the proposed algorithm and to compare with the original CBS method, and the results show that the proposed method has preferable accuracy with less numerical dissipation. We further applied the method to the numerical simulation of flow around equilateral triangular cylinder with different incidence angles in free stream. In this numerical investigation, the flow simulations are carried out in the low Reynolds number range. Instantaneous streamlines around the cylinder are used as a means to visualize the wake region behind, and they clearly show the flow pattern around the cylinder in time. The influence of incidence angle on flow characteristic parameters such as Strouhal number, Drag and Lift coefficients are discussed quantitatively. Copyright © 2009 John Wiley & Sons, Ltd.     

A cell‐based smoothed finite element method with semi‐implicit CBS procedures for incompressible laminar viscous flows

In this paper, the cell‐based smoothed finite element method (CS‐FEM) with the semi‐implicit characteristic‐based split (CBS) scheme (CBS/CS‐FEM) is proposed for computational fluid dynamics. The 3‐node triangular (T3) element and 4‐node quadrilateral (Q4) element are used for present CBS/CS‐FEM for two‐dimensional flows. The 8‐node hexahedral element (H8) is used for three‐dimensional flows. Two types of CS‐FEM are implemented in this paper. One is standard CS‐FEM with quadrilateral gradient smoothing cells for Q4 element and hexahedron cells for H8 element. Another is called as n‐sided CS‐FEM (nCS‐FEM) whose gradient smoothing cells are triangles for Q4 element and pyramids for H8 element. To verify the proposed methods, benchmarking problems are tested for two‐dimensional and three‐dimensional flows. The benchmarks show that CBS/CS‐FEM and CBS/nCS‐FEM are capable to solve incompressible laminar flow and can produce reliable results for both steady and unsteady flows. The proposed CBS/CS‐FEM method has merits on better robustness against distorted mesh with only slight more computation time and without losing accuracy, which is important for problems with heavy mesh distortion. The blood flow in carotid bifurcation is also simulated to show capabilities of proposed methods for realistic and complicated flow problems.     

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.     

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.     

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 stabilized finite element predictor–corrector scheme for the incompressible Navier–Stokes equations using a nodal‐based implementation

A finite element model to solve the incompressible Navier–Stokes equations based on the stabilization with orthogonal subscales, a predictor–corrector scheme to segregate the pressure and a nodal based implementation is presented in this paper. The stabilization consists of adding a least‐squares form of the component orthogonal to the finite element space of the convective and pressure gradient terms, which allows to deal with convection‐dominated flows and to use equal velocity–pressure interpolation. The pressure segregation is inspired in fractional step schemes, although the converged solution corresponds to that of a monolithic time integration. Finally, the nodal‐based implementation is based on an a priori calculation of the integrals appearing in the formulation and then the construction of the matrix and right‐hand side vector of the final algebraic system to be solved. After appropriate approximations, this matrix and this vector can be constructed directly for each nodal point, without the need to loop over the elements and thus making the calculations much faster. Some issues related to this implementation for fractional step and our predictor–corrector scheme, which is the main contribution of this paper, are discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

The multi‐stage centred‐scheme approach applied to a drift‐flux two‐phase flow model

For two‐phase flow models, upwind schemes are most often difficult do derive, and expensive to use. Centred schemes, on the other hand, are simple, but more dissipative. The recently proposed multi‐stage (MUSTA ) method is aimed at coming close to the accuracy of upwind schemes while retaining the simplicity of centred schemes. So far, the MUSTA approach has been shown to work well for the Euler equations of inviscid, compressible single‐phase flow. In this work, we explore the MUSTA scheme for a more complex system of equations: the drift‐flux model, which describes one‐dimensional two‐phase flow where the motions of the phases are strongly coupled. As the number of stages is increased, the results of the MUSTA scheme approach those of the Roe method. The good results of the MUSTA scheme are dependent on the use of a large‐enough local grid. Hence, the main benefit of the MUSTA scheme is its simplicity, rather than CPU ‐time savings. Copyright © 2006 John Wiley & Sons, Ltd.     

Similar solution for three‐dimensional flow in an Oldroyd‐B fluid over a stretching surface

The present investigation deals with the three‐dimensional flow of an Oldroyd‐B fluid over a stretching surface. The governing equations for the three‐dimensional flow are developed. Similarity transformations are invoked for the conversion of nonlinear partial differential equations into the coupled system of ordinary differential equations. Computations for the series solution are presented through implementation of homotopy analysis method. The salient features of the involved parameters have been pointed out. Copyright © 2011 John Wiley & Sons, Ltd.     

Characteristic‐based inlet and outlet boundary conditions for incompressible flows

Determining boundary conditions (BCs) for incompressible flows is such a delicate matter that affects the accuracy of the results. In this research, a new characteristic‐based BC for incompressible Navier‐Stokes equations is introduced. Discretization of equations has been done via finite volume. Additionally, artificial compressibility correction has been employed to deal with equations. Ordinary extrapolation from inner cells of a domain was used as a traditional way to estimate pressure and velocities on solid wall and inlet/outlet boundaries. Here, this method was substituted by the newly proposed BCs based on the characteristics of artificial compressibility equations. To follow this purpose, a computer code has been developed to carry out series of numerical tests for a flow over a backward‐facing step and was applied to a wide range of Reynolds numbers and grid combinations. Calculation of convective and viscous fluxes was done using Jameson's averaging scheme. Employing the characteristic‐based method for determining BCs has shown an improved convergence rate and reduced calculation time comparing with those of traditional ones. Furthermore, with the reduction of domain and computational cells, a similar accuracy was achieved for the results in comparison with the ones obtained from the traditional extrapolation method, and these results were in good agreement with the ones in the literature.     

Development of circular function‐based gas‐kinetic scheme (CGKS) on moving grids for unsteady flows through oscillating cascades

In this paper, the circular function‐based gas‐kinetic scheme (CGKS), which was originally developed for simulation of flows on stationary grids, is extended to solve moving boundary problems on moving grids. Particularly, the unsteady flows through oscillating cascades are our major interests. The main idea of the CGKS is to discretize the macroscopic equations by the finite volume method while the fluxes at the cell interface are evaluated by locally reconstructing the solution of the continuous Boltzmann Bhatnagar–Gross–Krook equation. The present solver is based on the fact that the modified Boltzmann equation, which is expressed in a moving frame of reference, can recover the corresponding macroscopic equations with Chapman–Enskog expansion analysis. Different from the original Maxwellian function‐based gas‐kinetic scheme, in improving the computational efficiency, a simple circular function is used to describe the equilibrium state of distribution function. Considering that the concerned cascade oscillating problems belong to cases that the motion of surface boundary is known a priori, the dynamic mesh method is suitable and is adopted in the present work. In achieving the mesh deformation with high quality and efficiency, a hybrid dynamic mesh method named radial basic functions‐transfinite interpolation is presented and applied for cascade geometries. For validation, several numerical test cases involving a wide range are investigated. Numerical results show that the developed CGKS on moving grids is well applied for cascade oscillating flows. And for some cases where nonlinear effects are strong, the solution accuracy could be effectively improved by using the present method. Copyright © 2017 John Wiley & Sons, Ltd.     

Efficient computation of mean drag for the subcritical flow past a circular cylinder using general Galerkin G2

General Galerkin (G2) is a new computational method for turbulent flow, where a stabilized Galerkin finite element method is used to compute approximate weak solutions to the Navier–Stokes equations directly, without any filtering of the equations as in a standard approach to turbulence simulation, such as large eddy simulation, and thus no Reynolds stresses are introduced, which need modelling. In this paper, G2 is used to compute the drag coefficient c_D for the flow past a circular cylinder at Reynolds number Re=3900, for which the flow is turbulent. It is found that it is possible to approximate c_D to an accuracy of a few percent, corresponding to the accuracy in experimental results for this problem, using less than 10⁵ mesh points, which makes the simulations possible using a standard PC. The mesh is adaptively refined until a stopping criterion is reached with respect to the error in a chosen output of interest, which in this paper is c_D. Both the stopping criterion and the mesh‐refinement strategy are based on a posteriori error estimates, in the form of a space–time integral of residuals times derivatives of the solution of a dual problem, linearized at the approximate solution, and with data coupling to the output of interest. Copyright © 2008 John Wiley & Sons, Ltd.     

A high‐order mass‐lumping procedure for B‐spline collocation method with application to incompressible flow simulations

This paper presents new developments of the staggered spline collocation method for cost‐effective solution to the incompressible Navier–Stokes equations. Maximal decoupling of the velocity and the pressure is obtained by using the fractional step method of Gresho and Chan, allowing the solution to sparse elliptic problems only. In order to preserve the high‐accuracy of the B‐spline method, this fractional step scheme is used in association with a sparse approximation to the inverse of the consistent mass matrix. Such an approximation is constructed from local spline interpolation method, and represents a high‐order generalization of the mass‐lumping technique of the finite‐element method. A numerical investigation of the accuracy and the computational efficiency of the resulting semi‐consistent spline collocation schemes is presented. These schemes generate a stable and accurate unsteady Navier–Stokes solver, as assessed by benchmark computations. Copyright © 2003 John Wiley & Sons, Ltd.     

Flux‐difference splitting‐based upwind compact schemes for the incompressible Navier–Stokes equations

Third‐order and fifth‐order upwind compact finite difference schemes based on flux‐difference splitting are proposed for solving the incompressible Navier–Stokes equations in conjunction with the artificial compressibility (AC) method. Since the governing equations in the AC method are hyperbolic, flux‐difference splitting (FDS) originally developed for the compressible Euler equations can be used. In the present upwind compact schemes, the split derivatives for the convective terms at grid points are linked to the differences of split fluxes between neighboring grid points, and these differences are computed by using FDS. The viscous terms are approximated with a sixth‐order central compact scheme. Comparisons with 2D benchmark solutions demonstrate that the present compact schemes are simple, efficient, and high‐order accurate. Copyright © 2008 John Wiley & Sons, Ltd.     

Three‐dimensional computation for flow‐induced vibrations in in‐line and cross‐flow directions of a circular cylinder

We present numerical results for in‐line and cross‐flow vibrations of a circular cylinder, which is immersed in a uniform flow and is elastically supported by damper‐spring systems to compute vibrations of a rigid cylinder. In the case of a circular cylinder with a low Scruton number, it is well‐known that two types of self‐excited vibrations appear in the in‐line direction in the range of low reduced velocities. On the other hand, a cross‐flow vibration of the circular cylinder can be excited in the range of high reduced velocities. Therefore, we compute the flow‐induced vibrations of the circular cylinder in the wide range of the reduced velocities at low and high Scruton numbers and discuss about excitation mechanisms in the in‐line and cross‐flow directions. Copyright © 2011 John Wiley & Sons, Ltd.