A non‐conforming least‐squares finite element method for incompressible fluid flow problems

In this paper, we develop least‐squares finite element methods (LSFEMs) for incompressible fluid flows with improved mass conservation. Specifically, we formulate a new locally conservative LSFEM for the velocity–vorticity–pressure Stokes system, which uses a piecewise divergence‐free basis for the velocity and standard C⁰ elements for the vorticity and the pressure. The new method, which we term dV‐VP improves upon our previous discontinuous stream‐function formulation in several ways. The use of a velocity basis, instead of a stream function, simplifies the imposition and implementation of the velocity boundary condition, and eliminates second‐order terms from the least‐squares functional. Moreover, the size of the resulting discrete problem is reduced because the piecewise solenoidal velocity element is approximately one‐half of the dimension of a stream‐function element of equal accuracy. In two dimensions, the discontinuous stream‐function LSFEM [1] motivates modification of our functional, which further improves the conservation of mass. We briefly discuss the extension of this modification to three dimensions. Computational studies demonstrate that the new formulation achieves optimal convergence rates and yields high conservation of mass. We also propose a simple diagonal preconditioner for the dV‐VP formulation, which significantly reduces the condition number of the LSFEM problem. Published 2012. This article is a US Government work and is in the public domain in the USA.     

Parallelization of a vorticity formulation for the analysis of incompressible viscous fluid flows

A parallel computer implementation of a vorticity formulation for the analysis of incompressible viscous fluid flow problems is presented. The vorticity formulation involves a three‐step process, two kinematic steps followed by a kinetic step. The first kinematic step determines vortex sheet strengths along the boundary of the domain from a Galerkin implementation of the generalized Helmholtz decomposition. The vortex sheet strengths are related to the vorticity flux boundary conditions. The second kinematic step determines the interior velocity field from the regular form of the generalized Helmholtz decomposition. The third kinetic step solves the vorticity equation using a Galerkin finite element method with boundary conditions determined in the first step and velocities determined in the second step. The accuracy of the numerical algorithm is demonstrated through the driven‐cavity problem and the 2‐D cylinder in a free‐stream problem, which represent both internal and external flows. Each of the three steps requires a unique parallelization effort, which are evaluated in terms of parallel efficiency. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of boundary vorticity discretization on explicit stream‐function vorticity calculations

The numerical solution of the time‐dependent Navier–Stokes equations in terms of the vorticity and a stream function is a well tested process to describe two‐dimensional incompressible flows, both for fluid mixing applications and for studies in theoretical fluid mechanics. In this paper, we consider the interaction between the unsteady advection–diffusion equation for the vorticity, the Poisson equation linking vorticity and stream function and the approximation of the boundary vorticity, examining from a practical viewpoint, global iteration stability and error. Our results show that most schemes have very similar global stability constraints although there may be small stability gains from the choice of method to determine boundary vorticity. Concerning accuracy, for one model problem we observe that there were cases where the boundary vorticity discretization did not propagate to the interior, but for the usual cavity flow all the schemes tested had error close to second order. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite volume solution of the Navier–Stokes equations in velocity–vorticity formulation

For the incompressible Navier–Stokes equations, vorticity‐based formulations have many attractive features over primitive‐variable velocity–pressure formulations. However, some features interfere with the use of the numerical methods based on the vorticity formulations, one of them being the lack of a boundary conditions on vorticity. In this paper, a novel approach is presented to solve the velocity–vorticity integro‐differential formulations. The general numerical method is based on standard finite volume scheme. The velocities needed at the vertexes of each control volume are calculated by a so‐called generalized Biot–Savart formula combined with a fast summation algorithm, which makes the velocity boundary conditions implicitly satisfied by maintaining the kinematic compatibility of the velocity and vorticity fields. The well‐known fractional step approaches are used to solve the vorticity transport equation. The paper describes in detail how we accurately impose no normal‐flow and no tangential‐flow boundary conditions. We impose a no‐flux boundary condition on solid objects by the introduction of a proper amount of vorticity at wall. The diffusion term in the transport equation is treated implicitly using a conservative finite update. The diffusive fluxes of vorticity into flow domain from solid boundaries are determined by an iterative process in order to satisfy the no tangential‐flow boundary condition. As application examples, the impulsively started flows through a flat plate and a circular cylinder are computed using the method. The present results are compared with the analytical solution and other numerical results and show good agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

Fourth‐order finite difference scheme for the incompressible Navier–Stokes equations in a disk

We develop an efficient fourth‐order finite difference method for solving the incompressible Navier–Stokes equations in the vorticity‐stream function formulation on a disk. We use the fourth‐order Runge–Kutta method for the time integration and treat both the convection and diffusion terms explicitly. Using a uniform grid with shifting a half mesh away from the origin, we avoid placing the grid point directly at the origin; thus, no pole approximation is needed. Besides, on such grid, a fourth‐order fast direct method is used to solve the Poisson equation of the stream function. By Fourier filtering the vorticity in the azimuthal direction at each time stage, we are able to increase the time step to a reasonable size. The numerical results of the accuracy test and the simulation of a vortex dipole colliding with circular wall are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

A high‐order accurate method for two‐dimensional incompressible viscous flows

A high‐order accurate solution method for complex geometries is developed for two‐dimensional flows using the stream function–vorticity formulation. High‐order accurate spectrally optimized compact schemes along with appropriate boundary schemes are used for spatial discretization while a two‐level backward Euler implicit scheme is used for the time integration. The linear system of equations for stream function and vorticity are solved by an inner iteration while contravariant velocities constitute outer iterations. The effect of curvilinear grids on the solution accuracy is studied. The method is used to compute Cartesian and inclined driven cavity, flow in a triangular cavity and viscous flow in constricted channel. Benchmark‐like accuracy is obtained in all the problems with fewer grid points compared to reported studies. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical study of the atmospheric flow over a coastal cliff

This work presents a numerical method suitable for the study of the development of internal boundary layers (IBL) and their characteristics for flows over various types of coastal cliffs. The IBL is an important meteorological occurrence for flows with surface roughness and topographical step changes. A two‐dimensional flow program was used for this study. The governing equations were written using the vorticity–velocity formulation. The spatial derivatives were discretized by high‐order compact finite differences schemes. The time integration was performed with a low storage fourth‐order Runge–Kutta scheme. The coastal cliff (step) was specified through an immersed boundary method. The validation of the code was done by comparison of the results with experimental and observational data. The numerical simulations were carried out for different coastal cliff heights and inclinations. The results show that the predominant factors for the height of the IBL and its characteristics are the upstream velocity, and the height and form (inclination) of the coastal cliff. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element solution of the Navier–Stokes equations by a velocity–vorticity method

A velocity–vorticity formulation of the Navier–Stokes equations is presented as an alternative to the primitive variables approach. The velocity components and the vorticity are solved for in a fully coupled manner using a Newton method. No artificial viscosity is required in this formulation. The pressure is updated by a method allowing natural imposition of boundary conditions. Incompressible and subsonic results are presented for two-dimensional laminar internal flows up to high Reynolds numbers.     

Dynamic vorticity condition: Theoretical analysis and numerical implementation

The dynamic boundary conditions for vorticity, derived from the incompressible Navier-Stokes equations, are examined from both theoretical and computational points of view. It is found that these conditions can be either local (Neumann type) or global (Dirichlet type), both containing coupling with the boundary pressure, which is the main difficulty in applying vorticity-based methods. An integral formulation is presented to analyse the structure of vorticity and pressure solutions, especially the strength of the coupling. We find that for high-Reynolds-number flows the coupling is weak and, if necessary, can be effectively bypassed by simple iteration. In fact, even a fully decoupled approximation is well applicable for most Reynolds numbers of practical interest. The fractional step method turns out to be especially appropriate for implementing the decoupled approximation. Both integral and finite difference methods are tested for some simple cases with known exact solutions. In the integral approach smoothed heat kernels are used to increase the accuracy of numerical quadrature. For the more complicated problem of impulsively started flow over a circular cylinder at Re = 9500 the finite difference method is used. The results are compared against numerical solutions and fine experiments with good agreement. These numerical experiments confirm our thoeretical analysis and show the advantages of the dynamic condition in computing high-Reynolds-number flows.     

The DRBEM solution of incompressible MHD flow equations

This paper presents a dual reciprocity boundary element method (DRBEM) formulation coupled with an implicit backward difference time integration scheme for the solution of the incompressible magnetohydrodynamic (MHD) flow equations. The governing equations are the coupled system of Navier‐Stokes equations and Maxwell's equations of electromagnetics through Ohm's law. We are concerned with a stream function‐vorticity‐magnetic induction‐current density formulation of the full MHD equations in 2D. The stream function and magnetic induction equations which are poisson‐type, are solved by using DRBEM with the fundamental solution of Laplace equation. In the DRBEM solution of the time‐dependent vorticity and current density equations all the terms apart from the Laplace term are treated as nonhomogeneities. The time derivatives are approximated by an implicit backward difference whereas the convective terms are approximated by radial basis functions. The applications are given for the MHD flow, in a square cavity and in a backward‐facing step. The numerical results for the square cavity problem in the presence of a magnetic field are visualized for several values of Reynolds, Hartmann and magnetic Reynolds numbers. The effect of each parameter is analyzed with the graphs presented in terms of stream function, vorticity, current density and magnetic induction contours. Then, we provide the solution of the step flow problem in terms of velocity field, vorticity, current density and magnetic field for increasing values of Hartmann number. Copyright © 2010 John Wiley & Sons, Ltd.     

Large eddy simulation (2D) using diffusion–velocity method and vortex‐in‐cell

A large eddy Simulation based on the diffusion‐velocity method and the discrete vortex method is presented. The vorticity‐based and eddy viscosity type subgrid scale model simulating the enstrophy transfer between the large and small scale appears as a convective term in the diffusion‐velocity formulation. The methodology has been tested on a spatially growing mixing layer using the two‐dimensional vortex‐in‐cell method and the Smagorinsky subgrid scale model. The effects on the vorticity contours, momemtum thickness, mean streamwise velocity profiles, root‐mean‐square velocity and vorticity fluctuations and negative cross‐stream correlation are discussed. Comparison is made with experiment and numerical work where diffusion is simulated using random walk. Copyright © 2004 John Wiley & Sons, Ltd.     

A viscous vortex particle method for deforming bodies with application to biolocomotion

Bio‐inspired mechanics of locomotion generally consist of the interaction of flexible structures with the surrounding fluid to generate propulsive forces. In this work, we extend, for the first time, the viscous vortex particle method (VVPM) to continuously deforming two‐dimensional bodies. The VVPM is a high‐fidelity Navier–Stokes computational method that captures the fluid motion through evolution of vorticity‐bearing computational particles. The kinematics of the deforming body surface are accounted for via a surface integral in the Biot–Savart velocity. The spurious slip velocity in each time step is removed by computing an equivalent vortex sheet and allowing it to flux to adjacent particles; hence, no‐slip boundary conditions are enforced. Particles of both uniform and variable size are utilized, and their relative merits are considered. The placement of this method in the larger class of immersed boundary methods is explored. Validation of the method is carried out on the problem of a periodically deforming circular cylinder immersed in a stagnant fluid, for which an analytical solution exists when the deformations are small. We show that the computed vorticity and velocity of this motion are both in excellent agreement with the analytical solution. Finally, we explore the fluid dynamics of a simple fish‐like shape undergoing undulatory motion when immersed in a uniform free stream, to demonstrate the application of the method to investigations of biomorphic locomotion. Copyright © 2008 John Wiley & Sons, Ltd.     

A fluid solver based on vorticity–helical density equations with application to a natural convection in a cubic cavity

We study numerically a recently introduced formulation of incompressible Newtonian fluid equations in vorticity–helical density and velocity–Bernoulli pressure variables. Unlike most numerical methods based on vorticity equations, the current approach provides discrete solutions with mass conservation, divergence‐free vorticity, and accurate kinetic energy balance in a simple and natural way. The method is applied to compute buoyancy‐driven flows in a differentially heated cubic enclosure in the Boussinesq approximation for Ra ∈ {10⁴,10⁵,10⁶}. The numerical solutions on a finer grid are of benchmark quality. The computed helical density allows quantification of the three‐dimensional nature of the flow. Copyright © 2011 John Wiley & Sons, Ltd.     

Large eddy simulation (2D) of spatially developing mixing layer using vortex‐in‐cell for flow field and filtered probability density function for scalar field

A large eddy simulation based on filtered vorticity transport equation has been coupled with filtered probability density function transport equation for scalar field, to predict the velocity and passive scalar fields. The filtered vorticity transport has been formulated using diffusion‐velocity method and then solved using the vortex method. The methodology has been tested on a spatially growing mixing layer using the two‐dimensional vortex‐in‐cell method in conjunction with both Smagorinsky and dynamic eddy viscosity subgrid scale models for an anisotropic flow. The transport equation for filtered probability density function is solved using the Lagrangian Monte‐Carlo method. The unresolved subgrid scale convective term in filtered density function transport is modelled using the gradient diffusion model. The unresolved subgrid scale mixing term is modelled using the modified Curl model. The effects of subgrid scale models on the vorticity contours, mean streamwise velocity profiles, root‐mean‐square velocity and vorticity fluctuations profiles and negative cross‐stream correlations are discussed. Also the characteristics of the passive scalar, i.e. mean concentration profiles, root‐mean‐square concentration fluctuations profiles and filtered probability density function are presented and compared with previous experimental and numerical works. The sensitivity of the results to the Schmidt number, constant in mixing frequency and inflow boundary conditions are discussed. Copyright © 2005 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.     

Numerical simulations of inviscid and viscous free‐surface flows with application to nonlinear water waves

The purpose of the present study is to establish a numerical model appropriate for solving inviscid/viscous free‐surface flows related to nonlinear water wave propagation. The viscous model presented herein is based on the Navier–Stokes equations, and the free‐surface is calculated through an arbitrary Lagrangian–Eulerian streamfunction‐vorticity formulation. The streamfunction field is governed by the Poisson equation, and the vorticity is obtained on the basis of the vorticity transport equation. For computing the inviscid flow the Laplace streamfunction equation is used. These equations together with the respective (appropriate) fully nonlinear free‐surface boundary conditions are solved using a finite difference method. To demonstrate the model feasibility, in the present study we first simulate collision processes of two solitary waves of different amplitudes, and compute the phenomenon of overtaking of such solitary waves. The developed model is subsequently applied to calculate (both inviscid and the viscous) flow field, as induced by passing of a solitary wave over submerged rectangular structures and rigid ripple beds. Our study provides a reasonably good understanding of the behavior of (inviscid/viscous) free‐surface flows, within the framework of streamfunction‐vorticity formulation. The successful simulation of the above‐mentioned test cases seems to suggest that the arbitrary Lagrangian–Eulerian/streamfunction‐vorticity formulation is a potentially powerful approach, capable of effectively solving the fully nonlinear inviscid/viscous free‐surface flow interactions. Copyright © 2009 John Wiley & Sons, Ltd.     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

An augmented velocity–vorticity–pressure formulation for the Brinkman equations

This paper deals with the analysis of a new augmented mixed finite element method in terms of vorticity, velocity, and pressure, for the Brinkman problem with nonstandard boundary conditions. The approach is based on the introduction of Galerkin least‐squares terms arising from the constitutive equation relating the aforementioned unknowns and from the incompressibility condition. We show that the resulting augmented bilinear form is continuous and elliptic, which, thanks to the Lax–Milgram theorem, and besides proving the well‐posedness of the continuous formulation, ensures the solvability and stability of the Galerkin scheme with any finite element subspace of the continuous space. In particular, Raviart–Thomas elements of any order for the velocity field, and piecewise continuous polynomials of degree k + 1 for both the vorticity and the pressure, can be utilized. A priori error estimates and the corresponding rates of convergence are also given here. Next, we derive two reliable and efficient residual‐based a posteriori error estimators for this problem. The ellipticity of the bilinear form together with the local approximation properties of the Clément interpolation operator are the main tools for showing the reliability. In turn, inverse inequalities and the localization technique based on triangle‐bubble and edge‐bubble functions are utilized to show the efficiency. Finally, several numerical results illustrating the good performance of the method, confirming the properties of the estimators and showing the behavior of the associated adaptive algorithms, are reported. Copyright © 2015 John Wiley & Sons, Ltd.     

An improved immersed boundary‐lattice Boltzmann method based on force correction technique

In this paper, an improved immersed boundary‐lattice Boltzmann method based on the force correction technique is presented for fluid‐structure interaction problems including the moving boundary interfaces. By introducing a force correction coefficient, the non‐slip boundary conditions are much better enforced compared with the conventional immersed boundary‐lattice Boltzmann methods. In addition, the implicit and iterative calculations are avoided; thus, the computational cost is reduced dramatically. Several numerical experiments are carried out to test the efficiency of the method. It is found that the method has the second‐order accuracy, and the non‐slip boundary conditions are enforced indeed. The numerical results also show that the present method is a suitable tool for fluid‐structure interaction problems involving complex moving boundaries.     

Solution to transient Navier–Stokes equations by the coupling of differential quadrature time integration scheme with dual reciprocity boundary element method

The two‐dimensional time‐dependent Navier–Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equation, which links stream function and vorticity with an initial vorticity guess, produces velocity components in turn for the solution to vorticity transport equation. The DRBEM formulation of the vorticity transport equation results in an initial value problem represented by a system of first‐order ordinary differential equations in time. When the DQM discretizes this system in time direction, we obtain a system of linear algebraic equations, which gives the solution vector for vorticity at any required time level. The procedure outlined here is also applied to solve the problem of two‐dimensional natural convection in a cavity by utilizing an iteration among the stream function, the vorticity transport and the energy equations as well. The test problems include two‐dimensional flow in a cavity when a force is present, the lid‐driven cavity and the natural convection in a square cavity. The numerical results are visualized in terms of stream function, vorticity and temperature contours for several values of Reynolds (Re) and Rayleigh (Ra) numbers. Copyright © 2008 John Wiley & Sons, Ltd.