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.     

Orthogonality of modal bases in hp finite element models

In this paper, we exploit orthogonality of modal bases (SIAM J. Sci. Comput. 1999; 20 :1671–1695) used in hp finite element models. We calculate entries of coefficient matrix analytically without using any numerical integration, which can be computationally very expensive. We use properties of Jacobi polynomials and recast the entries of the coefficient matrix so that they can be evaluated analytically. We implement this in the context of the least‐squares finite element model although this procedure can be used in other finite element formulations. In this paper, we only develop analytical expressions for rectangular elements. Spectral convergence of the L₂ least‐squares functional is verified using exact solution of Kovasznay flow. Numerical results for transient flow over a backward‐facing step are also presented. We also solve steady flow past a circular cylinder and show the reduction in computational cost using expressions developed herein. Copyright © 2007 John Wiley & Sons, Ltd.     

A stress‐based least‐squares finite‐element model for incompressible Navier–Stokes equations

In this paper we present a stress‐based least‐squares finite‐element formulation for the solution of the Navier–Stokes equations governing flows of viscous incompressible fluids. Stress components are introduced as independent variables to make the system first order. Continuity equation becomes an algebraic equation and is eliminated from the system with suitable modifications. The h and p convergence are verified using the exact solution of Kovasznay flow. Steady flow past a large circular cylinder in a channel is solved to test mass conservation. Transient flow over a backward‐facing step problem is solved on several meshes. Results are compared with that obtained using vorticity‐based first‐order formulation for both benchmark problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Nonisothermal two‐ and three‐dimensional flow simulations of inelastic and viscoelastic fluids by a finite‐volume method

This paper applies the finite‐volume method to computations of steady flows of viscous and viscoelastic incompressible fluids in complex two and three‐dimensional geometries. The materials adopted in the study obey different constitutive laws: Newtonian, purely viscous Carreau–Yasuda as also Upper‐Convected Maxwell and Phan‐Thien/Tanner differential models, with a Williams–Landel–Ferry (WLF) equation for temperature dependence. Specific analyses are made depending on the rheological model. A staggered grid is used for discretizing the equations and unknowns. Stockage possibilities allow us to solve problems involving a great number of degrees of freedom, up to 1 500 000 unknowns with a desk computer. In relation to the fluid properties, our numerical simulations provide flow characteristics for various 2D and 3D configurations and demonstrate the possibilities of the code to solve problems involving complex nonlinear constitutive equations with thermal effects. Copyright © 2009 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.     

An hp‐adaptive discontinuous Galerkin method for shallow water flows

An adaptive spectral/hp discontinuous Galerkin method for the two‐dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non‐conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p?1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h‐type refinement, the parent element is subdivided into four similar sibling elements. The time‐stepping is performed using a third‐order Runge–Kutta scheme. The performance of the hp‐adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p‐adaptivity is more efficient than h‐adaptivity with respect to degrees of freedom and computational time. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple mixture theory for ν Newtonian and generalized Newtonian constituents

This work presents the development of mathematical models based on conservation laws for a saturated mixture of ν homogeneous, isotropic, and incompressible constituents for isothermal flows. The constituents and the mixture are assumed to be Newtonian or generalized Newtonian fluids. Power law and Carreau–Yasuda models are considered for generalized Newtonian shear thinning fluids. The mathematical model is derived for a ν constituent mixture with volume fractions ${\phi_\alpha}$ using principles of continuum mechanics: conservation of mass, balance of momenta, first and second laws of thermodynamics, and principles of mixture theory yielding continuity equations, momentum equations, energy equation, and constitutive theories for mechanical pressures and deviatoric Cauchy stress tensors in terms of the dependent variables related to the constituents. It is shown that for Newtonian fluids with constant transport properties, the mathematical models for constituents are decoupled. In this case, one could use individual constituent models to obtain constituent deformation fields, and then use mixture theory to obtain the deformation field for the mixture. In the case of generalized Newtonian fluids, the dependence of viscosities on deformation field does not permit decoupling. Numerical studies are also presented to demonstrate this aspect. Using fully developed flow of Newtonian and generalized Newtonian fluids between parallel plates as a model problem, it is shown that partial pressures p _α of the constituents must be expressed in terms of the mixture pressure p. In this work, we propose ${p_\alpha=\phi_\alpha p}$ and ${\sum_\alpha^\nu p_\alpha = p}$ which implies ${\sum_\alpha^\nu \phi_\alpha = 1}$ which obviously holds. This rule for partial pressure is shown to be valid for a mixture of Newtonian and generalized Newtonian constituents yielding Newtonian and generalized Newtonian mixture. Modifications of the currently used constitutive theories for deviatoric Cauchy stress tensor are proposed. These modifications are demonstrated to be essential in order for the mixture theory for ν constituents to yield a valid mathematical model when the constituents are the same. Dimensionless form of the mathematical models is derived and used to present numerical studies for boundary value problems using finite element processes based on a residual functional, that is, least squares finite element processes in which local approximations are considered in ${H^{k,p}\left(\bar{\Omega}^e\right)}$ scalar product spaces. Fully developed flow between parallel plates and 1:2 asymmetric backward facing step is used as model problems for a mixture of two constituents.     

An hp‐adaptive unstructured finite volume solver for compressible flows

The idea of hp‐adaptation, which has originally been developed for compact schemes (such as finite element methods), suggests an adaptation scheme using a mixture of mesh refinement and order enrichment based on the smoothness of the solution to obtain an accurate solution efficiently. In this paper, we develop an hp‐adaptation framework for unstructured finite volume methods using residual‐based and adjoint‐based error indicators. For the residual‐based error indicator, we use a higher‐order discrete operator to estimate the truncation error, whereas this estimate is weighted by the solution of the discrete adjoint problem for an output of interest to form the adaptation indicator for adjoint‐based adaptations. We perform our adaptation by local subdivision of cells with nonconforming interfaces allowed and local reconstruction of higher‐order polynomials for solution approximations. We present our results for two‐dimensional compressible flow problems including subsonic inviscid, transonic inviscid, and subsonic laminar flow around the NACA 0012 airfoil and also turbulent flow over a flat plate. Our numerical results suggest the efficiency and accuracy advantages of adjoint‐based hp‐adaptations over uniform refinement and also over residual‐based adaptation for flows with and without singularities.     

Augmented mixed finite element methods for a vorticity‐based velocity–pressure–stress formulation of the Stokes problem in 2D

In this paper, we consider an augmented velocity–pressure–stress formulation of the 2D Stokes problem, in which the stress is defined in terms of the vorticity and the pressure, and then we introduce and analyze stable mixed finite element methods to solve the associated Galerkin scheme. In this way, we further extend similar procedures applied recently to linear elasticity and to other mixed formulations for incompressible fluid flows. Indeed, our approach is based on the introduction of the Galerkin least‐squares‐type terms arising from the corresponding constitutive and equilibrium equations, and from the Dirichlet boundary condition for the velocity, all of them multiplied by stabilization parameters. Then, we show that these parameters can be suitably chosen so that the resulting operator equation induces a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well posed for any choice of finite element subspaces. In particular, we can use continuous piecewise linear velocities, piecewise constant pressures, and rotated Raviart–Thomas elements for the stresses. Next, we derive reliable and efficient residual‐based a posteriori error estimators for the augmented mixed finite element schemes. In addition, several numerical experiments illustrating the performance of the augmented mixed finite element methods, confirming the properties of the a posteriori estimators, and showing the behavior of the associated adaptive algorithms are reported. The present work should be considered as a first step aiming finally to derive augmented mixed finite element methods for vorticity‐based formulations of the 3D Stokes problem. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of spectral and finite element methods applied to the study of the core‐annular flow in an undulating tube

A Galerkin/finite element and a pseudo‐spectral method, in conjunction with the primitive (velocity‐pressure) and streamfunction‐vorticity formulations, are tested for solving the two‐phase flow in a tube, which has a periodically varying, circular cross section. Two immiscible, incompressible, Newtonian fluids are arranged so that one of them is around the axis of the tube (core fluid) and the other one surrounds it (annular fluid). The physical and flow parameters are such that the interface between the two fluids remains continuous and single‐valued. This arrangement is usually referred to as Core‐Annular flow. A non‐orthogonal mapping is used to transform the uneven tube shape and the unknown, time dependent interface to fixed, cylindrical surfaces. With both methods and formulations, steady states are calculated first using the Newton–Raphson method. The most dangerous eigenvalues of the related linear stability problem are calculated using the Arnoldi method, and dynamic simulations are carried out using the implicit Euler method. It is shown that with a smooth tube shape the pseudo‐spectral method exhibits exponential convergence, whereas the finite element method exhibits algebraic convergence, albeit of higher order than expected from the relevant theory. Thus the former method, especially when coupled with the streamfunction‐vorticity formulation, is much more efficient. The finite element method becomes more advantageous when the tube shape contains a cusp, in which case the convergence rate of the pseudo‐spectral method deteriorates exhibiting algebraic convergence with the number of the axial spectral modes, whereas the convergence rate of the finite element method remains unaffected. Copyright © 2002 John Wiley & Sons, Ltd.     

Simulation of pressure‐tooling wire‐coating flow with Phan‐Thien/Tanner models

Annular pressure‐tooling extrusion is simulated for a low density polymer melt using a Taylor–Petrov–Galerkin finite element scheme. This represents industrial‐scale wire‐coating. Viscoelastic fluids are modeled via three forms of Phan‐Thien/Tanner (PTT) constitutive laws employed for short‐die and full specification pressure‐tooling. Effects of variation in Weissenberg number (We) and polymeric viscosity are investigated. Particular attention is paid to mesh refinement to predict accurate results. The impact of variation in shear‐thinning and strain‐softening properties is considered upon the modelling predictions. For the short‐die flow, the influence of the lack of strain softening is identified. For the full‐die flow and more severe deformation rates, the linear PTT model failed to converge. In contrast, the exponential PTT model is found to be more stable numerically and to adequately reflect the material response. Comparing short‐die and full‐die pressure‐tooling results, shear rates increase 10‐fold, while strain rates increase one hundred times. Copyright © 2002 John Wiley & Sons, Ltd.     

A second‐order time accurate finite element method for quasi‐incompressible viscous flows

A finite element method for quasi‐incompressible viscous flows is presented. An equation for pressure is derived from a second‐order time accurate Taylor–Galerkin procedure that combines the mass and the momentum conservation laws. At each time step, once the pressure has been determined, the velocity field is computed solving discretized equations obtained from another second‐order time accurate scheme and a least‐squares minimization of spatial momentum residuals. The terms that stabilize the finite element method (controlling wiggles and circumventing the Babuska–Brezzi condition) arise naturally from the process, rather than being introduced a priori in the variational formulation. A comparison between the present second‐order accurate method and our previous first‐order accurate formulation is shown. The method is also demonstrated in the computation of the leaky‐lid driven cavity flow and in the simulation of a crossflow past a circular cylinder. In both cases, good agreement with previously published experimental and computational results has been obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

A further work on multi‐phase two‐fluid approach for compressible multi‐phase flows

This paper is to continue our previous work Niu (Int. J. Numer. Meth. Fluids 2001; 36 :351–371) on solving a two‐fluid model for compressible liquid–gas flows using the AUSMDV scheme. We first propose a pressure–velocity‐based diffusion term originally derived from AUSMDV scheme Wada and Liou (SIAM J. Sci. Comput. 1997; 18 (3):633—657) to enhance its robustness. The scheme can be applied to gas and liquid fluids universally. We then employ the stratified flow model Chang and Liou (J. Comput. Physics 2007; 225 :240–873) for spatial discretization. By defining the fluids in different regions and introducing inter‐phasic force on cell boundary, the stratified flow model allows the conservation laws to be applied on each phase, and therefore, it is able to capture fluid discontinuities, such as the fluid interfaces and shock waves, accurately. Several benchmark tests are studied, including the Ransom's Faucet problem, 1D air–water shock tube problems, 2D shock‐water column and 2D shock‐bubble interaction problems. The results indicate that the incorporation of the new dissipation into AUSM⁺‐up scheme and the stratified flow model is simple, accurate and robust enough for the compressible multi‐phase flows. Copyright © 2008 John Wiley & Sons, Ltd.     

Evaluation of Smagorinsky‐based subgrid‐scale models in a finite‐volume computation

Smagorinsky‐based models are assessed in a turbulent channel flow simulation at Re_b=2800 and Re_b=12500. The Navier–Stokes equations are solved with three different grid resolutions by using a co‐located finite‐volume method. Computations are repeated with Smagorinsky‐based subgrid‐scale models. A traditional Smagorinsky model is implemented with a van Driest damping function. A dynamic model assumes a similarity of the subgrid and the subtest Reynolds stresses and an explicit filtering operation is required. A top‐hat test filter is implemented with a trapezoidal and a Simpson rule. At the low Reynolds number computation none of the tested models improves the results at any grid level compared to the calculations with no model. The effect of the subgrid‐scale model is reduced as the grid is refined. The numerical implementation of the test filter influences on the result. At the higher Reynolds number the subgrid‐scale models stabilize the computation. An analysis of an accurately resolved flow field reveals that the discretization error overwhelms the subgrid term at Re_b=2800 in the most part of the computational domain. Copyright © 2002 John Wiley & Sons, Ltd.     

Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress model

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three‐dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non‐staggered grid arrangement. The pressure–velocity coupling is handled by using the SIMPLEC‐algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central‐difference scheme. The hybrid scheme is used for solving the ε equation. The secondary flow generation using the RSM model is compared with a non‐linear k–ε model (non‐linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of incompressible two‐phase flows with a Boussinesq–Scriven interface stress tensor

We consider the numerical simulation of a three‐dimensional two‐phase incompressible flow with a viscous interface. The simulation is based on a sharp interface Navier–Stokes model and the Boussinesq–Scriven constitutive law for the interface viscous stress tensor. In the recent paper [Soft Matter 7, 7797–7804, 2011], a model problem with a spherical droplet in a Stokes Poiseuille flow with a Boussinesq–Scriven law for the surface viscosity has been analyzed. In that paper, relations for the droplet migration velocity are derived. We relate the results obtained with our numerical solver for the two‐phase Navier–Stokes model to these theoretical relations. Copyright © 2013 John Wiley & Sons, Ltd.     

A sensitivity analysis on the parameter of the GLS method for a second‐gradient theory of incompressible flow

Using a non‐conforming C⁰‐interior penalty method and the Galerkin least‐square approach, we develop a continuous–discontinuous Galerkin finite element method for discretizing fourth‐order incompressible flow problems. The formulation is weakly coercive for spaces that fail to satisfy the inf‐sup condition and consider discontinuous basis functions for the pressure field. We consider the results of a stability analysis through a lemma which indicates that there exists an optimal or quasi‐optimal least‐square stability parameter that depends on the polynomial degree used to interpolate the velocity and pressure fields, and on the geometry of the finite element in the mesh. We provide several numerical experiments illustrating such dependence, as well as the robustness of the method to deal with arbitrary basis functions for velocity and pressure, and the ability to stabilize large pressure gradients. We believe the results provided in this paper contribute for establishing a paradigm for future studies of the parameter of the Galerkin least square method for second‐gradient theory of incompressible flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. Copyright © 2004 John Wiley & Sons, Ltd.     

Spectral p‐multigrid discontinuous Galerkin solution of the Navier–Stokes equations

Discontinuous Galerkin (DG) methods are very well suited for the construction of very high‐order approximations of the Euler and Navier–Stokes equations on unstructured and possibly nonconforming grids, but are rather demanding in terms of computational resources. In order to improve the computational efficiency of this class of methods, a high‐order spectral element DG approximation of the Navier–Stokes equations coupled with a p‐multigrid solution strategy based on a semi‐implicit Runge–Kutta smoother is considered here. The effectiveness of the proposed approach in the solution of compressible shockless flow problems is demonstrated on 2D inviscid and viscous test cases by comparison with both a p‐multigrid scheme with non‐spectral elements and a spectral element DG approach with an implicit time integration scheme. Copyright © 2010 John Wiley & Sons, Ltd.