A new stable space–time formulation for two‐dimensional and three‐dimensional incompressible viscous flow

A space–time finite element method for the incompressible Navier–Stokes equations in a bounded domain in ?^d (with d=2 or 3) is presented. The method is based on the time‐discontinuous Galerkin method with the use of simplex‐type meshes together with the requirement that the space–time finite element discretization for the velocity and the pressure satisfy the inf–sup stability condition of Brezzi and Babu?ka. The finite element discretization for the pressure consists of piecewise linear functions, while piecewise linear functions enriched with a bubble function are used for the velocity. The stability proof and numerical results for some two‐dimensional problems are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

An extended κ–ϵ finite element model

An extended κ–? model (to include low-Reynolds-number regions) employing weighting functions is presented. Wall functions for the near-wall zones are developed giving correct boundary values for the Shear stress and κ–?. A finite element model using a penalty formulation for incompressible turbulent flow is applied to Solve a flow between two plates. Results with mesh boundaries situated in the near-wall region and a: the wall are compared with measured values.     

Simulating two-dimensional turbulent flow by using the κ–ϵ model and the vorticity–streamfunction formulation

A numerical procedure to solve turbulent flow which makes use of the κ–? model has been developed. The method is based on a control volume finite element method and an unstructured triangular domain discretization. The velocity-pressure coupling is addressed via the vorticity-streamfunction and special attention is given to the boundary conditions for the vorticity. Wall effects are taken into account via wail functions or a low-Reynolds-number model. The latter was found to perform better in recirculation regions. Source terms of the κ and ε transport equations have been linearized in a particular way to avoid non-realistic solutions. The vorticity and streamfunction discretized equations are solved in a coupled way to produce a faster and more stable computational procedure. Comparison between the numerical predictions and experimental data shows that the physics of the flow is correctly simulated.     

SOME VALIDATION OF STANDARD,MODIFIED AND NON-LINEAR k–ε TURBULENCE MODELS

Standard, modified and non-linear k–ε: turbulence models are validated against three axisymmetric flow problems—flow through a pipe expansion, flow through a pipe constriction and an impinging jet problem—to underpin knowledge about the solution quality obtained from two-equation turbulence models. The extended models improve the prediction of turbulence as a flow approaches a stagnation point and the non-linear model allows for the prediction of anisotropic turbulence. Significantly different values for the non-linear model coefficients are proposed in comparison with values found in the literature. Nevertheless, current turbulence models are still unable to accurately predict the spreading rate of shear layers. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids, 24: 965–986, 1997.     

AN ADAPTIVE FINITE ELEMENT METHOD FOR A TWO-EQUATION TURBULENCE MODEL IN WALL-BOUNDED FLOWS

This paper presents an adaptive finite element method for solving incompressible turbulent flows using a k–ϵ model of turbulence. Solutions are obtained in primitive variables using a highly accurate quadratic finite element on unstructured grids. A projection error estimator is presented that takes into account the relative importance of the errors in velocity, pressure and turbulence variables. The efficiency and convergence rate of the methodology are evaluated by solving problems with known analytical solutions. The method is then applied to turbulent flow over a backward-facing step and predictions are compared with experimental measurements. © 1997 John Wiley & Sons, Ltd.     

Upwind basis finite elements for convection-dominated problems

Finite elements using higher-order basis functions in the spirit of the QUICK method for convection-dominated fluid flow and transport problems are introduced and demonstrated. Instead of introducing new internal degrees of freedom, completeness is achieved by including functions based on nodal values exterior and upwind to the element domain. Applied with linear test functions to the weak statements for convection-dominated problems, a family of Petrov–Galerkin finite elements is developed. Quadratic and cubic versions are demonstrated for the one-dimensional convection–diffusion test problem. Elements of up to seventh degree are used for local solution refinement. The behaviour of these elements for one-dimensional linear and non-linear advection is investigated. A two-dimensional quadratic upwind element is demonstrated in a streamfunction–vorticity formulation of the Navier–Stokes equations for a driven cavity flow test problem. With some minor reservations, these elements are recommended for further study and application.     

Numerical implementation of the Crank–Nicolson/Adams–Bashforth scheme for the time‐dependent Navier–Stokes equations

This article considers numerical implementation of the Crank–Nicolson/Adams–Bashforth scheme for the two‐dimensional non‐stationary Navier–Stokes equations. A finite element method is applied for the spatial approximation of the velocity and pressure. The time discretization is based on the Crank–Nicolson scheme for the linear term and the explicit Adams–Bashforth scheme for the nonlinear term. Comparison with other methods, through a series of numerical experiments, shows that this method is almost unconditionally stable and convergent, i.e. stable and convergent when the time step is smaller than a given constant. Copyright © 2009 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.     

Lattice‐Boltzmann and meshless point collocation solvers for fluid flow and conjugate heat transfer

The applicability and performance of the lattice‐Boltzmann (LB) and meshless point collocation methods as CFD solvers in flow and conjugate heat transfer processes are investigated in this work. Lid‐driven cavity flow and flow in a slit with an obstacle including heat transfer are considered as case studies. A comparison of the computational efficiency accuracy of the two methods with that of a finite volume method as implemented in a commercial package (ANSYS CFX, ANSYS Inc., Canonsburg, PA) is made. Utilizing the analogy between heat and mass transfer, an advection–diffusion LB model was adopted to simulate the heat transfer part of the slit flow problem followed by a rigorous mapping of the mass transfer variables to the heat transfer quantities of interest, thus circumventing the need for a thermal LB model. Direct comparison among the results of the three methods revealed excellent agreement over a wide range of Reynolds and Prandtl number values. Furthermore, an integrated computational scheme is proposed, utilizing the rapid convergence of the LB model in the flow part of the conjugate heat transfer problem with that of the meshless collocation method for the heat transfer part. The meshless treatment remains sufficiently rapid even for conduction‐controlled processes in contrast to the LB method, which is very rapid in the convection‐controlled case only. A single, common computational grid, composed of regularly distributed nodes is used, saving significant computational and coding time and ensuring convergence of the discrete Laplacian operator in the heat transfer part of the computations. Copyright © 2012 John Wiley & Sons, Ltd.     

A New Model for the Dynamics of Dispersions in a Batch Reactor: Numerical Approach

In this paper, we develop some considerations concerning the structure of coalescence, breakage and volume scattering kernels appearing in the evolution equation related to a new model for the dynamics of liquid–liquid dispersions and show some numerical simulations. The mathematical model has been presented in [3, 4], where a proof of the existence and uniqueness for a classical solution to the integro–differential equation describing the physical phenomenon is provided as well as a complete analysis of the general characteristics of the integral kernels. Numerical simulations agree with experimental data and with the expected asymptotical behavior of the solution.     

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.     

A reappraisal of Taylor–Galerkin algorithm for drying–wetting areas in shallow water computations

Non‐linear shallow water equations are a useful approximation to phenomena such as estuary dynamics, tidal propagation, breaking of a dam, flood waves, etc. Quite frequently they involve propagation over dry beds and drying of wet zones, for which boundary changes. To solve this problem either special techniques such as remeshing and Arbitrary Lagrangian Eulerian formulations or algorithms initially developed for flows exhibiting shocks have been proposed in past years. The purpose of this paper is to show how classical finite elements formulations, such as Taylor–Galerkin can be applied to solve the problem to wetting–drying areas in a simple yet efficient manner. Copyright © 2002 John Wiley & Sons, Ltd.     

Asymptotic–Newton method for solving incompressible flows

In this paper we present a comparative study of three non-linear schemes for solving finite element systems of Navier–Stokes incompressible flows. The first scheme is the classical Newton–Raphson linearization, the second one is the modified Newton–Raphson linearization and the last one is a new scheme called the asymptotic–Newton method. The relative efficiency of these approaches is evaluated over a large number of examples. © 1997 John Wiley & Sons, Ltd.     

Development of a convection–diffusion‐reaction magnetohydrodynamic solver on non‐staggered grids

This paper presents a convection–diffusion‐reaction (CDR) model for solving magnetic induction equations and incompressible Navier–Stokes equations. For purposes of increasing the prediction accuracy, the general solution to the one‐dimensional constant‐coefficient CDR equation is employed. For purposes of extending this discrete formulation to two‐dimensional analysis, the alternating direction implicit solution algorithm is applied. Numerical tests that are amenable to analytic solutions were performed in order to validate the proposed scheme. Results show good agreement with the analytic solutions and high rate of convergence. Like many magnetohydrodynamic studies, the Hartmann–Poiseuille problem is considered as a benchmark test to validate the code. Copyright © 2004 John Wiley & Sons, Ltd.     

Utilization of bend–twist coupling for performance enhancement of composite marine propellers

Self-twisting composite marine propellers, when subject to hydrodynamic loading, will not only automatically bend but also twist due to passive bend–twist (BT) coupling characteristics of anisotropic composites. To exploit the BT coupling effects of self-twisting propellers, a two-level (material and geometry) design methodology is proposed, formulated, and implemented. The material design is formulated as a constrained, discrete, binary optimization problem, which is tackled using an enhanced genetic algorithm equipped with numerical and analytical tools as function evaluators. The geometry design is formulated as an inverse problem to determine the unloaded geometry, which is solved using an over-relaxed, nonlinear, iterative procedure. A sample design is provided to illustrate the design methodology, and the predicted performance is compared to that of a rigid propeller. The results show that the self-twisting propeller produced the same performance as the rigid propeller at the design flow condition, and it produced better performance than the rigid propeller at off-design flow conditions, including behind a spatially varying wake.     

Shallow-water flow solver with non-hydrostatic pressure: 2D vertical plane problems

A numerical solution for shallow-water flow is developed based on the unsteady Reynolds-averaged Navier–Stokes equations without the conventional assumption of hydrostatic pressure. Instead, the non-hydrostatic pressure component may be added in regions where its influence is significant, notably where bed slope is not small and separation in a vertical plane may occur or where the free-surface slope is not small. The equations are solved in the σ-co-ordinate system with semi-implicit time stepping and the eddy viscosity is calculated using the standard k–ϵ turbulence model. Conventionally, boundary conditions at the bed for shallow-water models only include vertical diffusion terms using wall functions, but here they are extended to include horizontal diffusion terms which can be significant when bed slope is not small. This is consistent with the inclusion of non-hydrostatic pressure. The model is applied to the 2D vertical plane flow of a current over a trench for which experimental data and other numerical results are available for comparison. Computations with and without non-hydrostatic pressure are compared for the same trench and for trenches with smaller side slopes, to test the range of validity of the conventional hydrostatic pressure assumption. The model is then applied to flow over a 2D mound and again the slope of the mound is reduced to assess the validity of the hydrostatic pressure assumption. © 1998 John Wiley & Sons, Ltd.     

Finite element simulation of turbulent Couette–Poiseuille flows using a low Reynolds number k–ϵ model

Developing Couette–Poiseuille flows at Re=5000 are studied using a low Reynolds number k–ϵ two‐equation model and a finite element formulation. Mesh‐independent solutions are obtained using a standard Galerkin formulation and a Galerkin/least‐squares stabilized method. The predictions for the velocity and turbulent kinetic energy are compared with available experimental results and to the DNS data. Second moment closure's solutions are also compared with those of the k–ϵ model. The deficiency of eddy viscosity models to predict dissymmetric low Reynolds number channel flows has been demonstrated. Copyright © 1999 John Wiley & Sons, Ltd.     

Computation of particle dispersion in turbulent liquid flows using an efficient Lagrangian trajectory model

The dispersion of solid particles in a turbulent liquid flow impinging on a centrebody through an axisymmetric sudden expansion was investigated numerically using a Eulerian–Lagrangian model. Detailed experimental measurements at the inlet were used to specify the inlet conditions for two-phase flow computations. The anisotropy of liquid turbulence was accounted for using a second-moment Reynold stress transport model. A recently developed stochastic–probabilistic model was used to enhance the computational efficiency of Lagrangian trajectory computations. Numerical results of the stochastic–probabilistic model using 650 particle trajectories were compared with those of the conventional stochastic discrete-delta-function model using 18 000 particle trajectories. In addition, results of the two models were compared with experimental measurements. © 1998 John Wiley & Sons, Ltd.     

Comment on ‘generalized potential flow theory and direct calculation of velocities applied to the numerical solution of the Navier–Stokes and the Boussinesq equations’

In a recent paper a generalized potential flow theory and its application to the solution of the Navier–Stokes equation are developed.¹ The purpose of this comment is to show that the analysis presented in that paper is in general not correct. We note that the theoretical development of Reference 1 is in fact an extension—although not cited—of some work first done by Hawthorne for steady inviscid flow.² Hawthorne's solution is correct, and his analysis, which we briefly describe, provides a useful introduction to this note.