A finite element thermally coupled flow formulation for phase‐change problems

A finite element, thermally coupled incompressible flow formulation considering phase‐change effects is presented. This formulation accounts for natural convection, temperature‐dependent material properties and isothermal and non‐isothermal phase‐change models. In this context, the full Navier–Stokes equations are solved using a generalized streamline operator (GSO) technique. The highly non‐linear phase‐change effects are treated with a temperature‐based algorithm, which provides stability and convergence of the numerical solution. The Boussinesq approximation is used in order to consider the temperature‐dependent density variation. Furthermore, the numerical solution of the coupled problem is approached with a staggered incremental‐iterative solution scheme, such that the convergence criteria are written in terms of the residual vectors. Finally, this formulation is used for the solutions of solidification and melting problems validating some numerical results with other existing solutions obtained with different methodologies. Copyright © 2000 John Wiley & Sons, Ltd.     

A class of higher order compact schemes for the unsteady two‐dimensional convection–diffusion equation with variable convection coefficients

A class of higher order compact (HOC) schemes has been developed with weighted time discretization for the two‐dimensional unsteady convection–diffusion equation with variable convection coefficients. The schemes are second or lower order accurate in time depending on the choice of the weighted average parameter μ and fourth order accurate in space. For 0.5?μ?1, the schemes are unconditionally stable. Unlike usual HOC schemes, these schemes are capable of using a grid aspect ratio other than unity. They efficiently capture both transient and steady solutions of linear and nonlinear convection–diffusion equations with Dirichlet as well as Neumann boundary condition. They are applied to one linear convection–diffusion problem and three flows of varying complexities governed by the two‐dimensional incompressible Navier–Stokes equations. Results obtained are in excellent agreement with analytical and established numerical results. Overall the schemes are found to be robust, efficient and accurate. Copyright © 2002 John Wiley & Sons, Ltd.     

A short note on the general boundary element method for viscous flows with high Reynolds number

In this paper, the general boundary element method and the parallel computation are employed to solve laminar viscous flows in a driven square cavity, governed by the exact Navier–Stokes equations. Using the solution at Re=0 as the initial approximation, the convergent numerical results for high Reynolds number at Re=7500 are obtained, for the first time, by the boundary element method. This verifies the validity and great potential of the general boundary element method for highly non‐linear problems, which may greatly enlarge application regions of the boundary element method in science and engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

A collocated grid,projection method for time‐accurate calculation of low‐Mach number variable density flows in general curvilinear coordinates

A time‐accurate algorithm is proposed for low‐Mach number, variable density flows on curvilinear grids. Spatial discretization is performed on collocated grid that offers computational simplicity in curvilinear coordinates. The flux interpolation technique is used to avoid the pressure odd–even decoupling of the collocated grid arrangement. To increase the stability of the method, a two‐step predictor–corrector time integration scheme is employed. At each step, the projection method is used to calculate the hydrodynamic pressure and to satisfy the continuity equation. The robustness and accuracy of the method is illustrated with a series of numerical experiments including thermally driven cavity, polar cavity, three‐dimensional cavity, and direct numerical simulation of non‐isothermal turbulent channel flow. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical study of projection algorithms in the finite element simulation of three‐dimensional viscous incompressible flow

In this work a comparative study of two versions of the projection algorithm used either for time integration or as an iterative method to solve the three‐dimensional incompressible Navier–Stokes equations is presented. It is also shown that these projection algorithms combined with the finite element method are particularly suited for the treatment of outflow boundary conditions in the context of external flows. This assertion is illustrated by means of some numerical examples in which five types of boundary conditions are compared. The scheme is applied to simulate the flow past a cylinder clamped on two fixed parallel solid walls. Comparison with experimental data available for this problem shows good agreement of the velocity and pressure fields, both computed with continuous piecewise linear elements. Copyright © 1999 John Wiley & Sons, Ltd.     

Pseudostress–velocity formulation for incompressible Navier–Stokes equations

This paper presents a numerical algorithm using the pseudostress–velocity formulation to solve incompressible Newtonian flows. The pseudostress–velocity formulation is a variation of the stress–velocity formulation, which does not require symmetric tensor spaces in the finite element discretization. Hence its discretization is greatly simplified. The discrete system is further decoupled into an H ( div ) problem for the pseudostress and a post‐process resolving the velocity. This can be done conveniently by using the penalty method for steady‐state flows or by using the time discretization for nonsteady‐state flows. We apply this formulation to the 2D lid‐driven cavity problem and study its grid convergence rate. Also, computational results of the time‐dependent‐driven cavity problem and the flow past rectangular problem are reported. Copyright © 2009 John Wiley & Sons, Ltd.     

Preconditioned Krylov subspace methods used in solving two‐dimensional transient two‐phase flows

This paper investigates the performance of preconditioned Krylov subspace methods used in a previously presented two‐fluid model developed for the simulation of separated and intermittent gas–liquid flows. The two‐fluid model has momentum and mass balances for each phase. The equations comprising this model are solved numerically by applying a two‐step semi‐implicit time integration procedure. A finite difference numerical scheme with a staggered mesh is used. Previously, the resulting linear algebraic equations were solved by a Gaussian band solver. In this study, these algebraic equations are also solved using the generalized minimum residual (GMRES) and the biconjugate gradient stabilized (Bi‐CGSTAB) Krylov subspace iterative methods preconditioned with incomplete LU factorization using the ILUT(p, τ) algorithm. The decrease in the computational time using the iterative solvers instead of the Gaussian band solver is shown to be considerable. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical heat transfer in a cavity with a solar control coating deposited to a vertical semitransparent wall

A transient two‐dimensional computational model of combined natural convection, conduction, and radiation in a cavity with an aspect ratio of one, containing air as a laminar and non‐participating fluid, is presented. The cavity has two opaque adiabatic horizontal walls, one opaque isothermal vertical wall, and an opposite semitransparent wall, which consists of a 6‐mm glass sheet with a solar control coating of SnS–Cu_xS facing the cavity. The semitransparent wall also exchanges heat by convection and radiation from its external surface to the surroundings and allows solar radiation pass through into the interior of the cavity. The momentum and energy equations in the transient state were solved by finite differences using the alternating direction implicit (ADI) technique. The transient conduction equation and the radiative energy flux boundary conditions are coupled to these equations. The results in this paper are limited to the following conditions: 10⁴≤Gr≤10⁶, an isothermal vertical cold wall of 21°C, outside air temperatures in the range 30°C≤T₀≤40°C and incident solar radiation of AM2 (750 W m⁻²) normal to the semitransparent wall. The model allows calculation of the redistribution of the absorbed component of solar radiation to the inside and outside of the cavity. The influences of the time step and mesh size were considered. Using arguments of energy balance in the cavity, it was found that the percentage difference was less than 4 per cent, showing a possible total numerical error less than this number. For Gr=10⁶ a wave appeared in the upper side of the cavity, suggesting the influence of the boundary walls over the air flow inside the cavity. A Nusselt number correlation as a function of the Rayleigh number is presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Spectral analysis of flux vector splitting finite volume methods

New results are presented here for finite volume (FV) methods that use flux vector splitting (FVS) along with higher‐order reconstruction schemes. Apart from spectral accuracy of the resultant methods, the numerical stability is investigated which restricts the allowable time step or the Courant–Friedrich–Lewy (CFL) number. Also the dispersion relation preservation (DRP) property of various spatial and temporal discretization schemes is investigated. The DRP property simultaneously fixes space and time steps. This aspect of numerical schemes is important for simulation of high‐Reynolds number flows, compressible flows with shock(s) and computational aero‐acoustics. It is shown here that for direct numerical simulation applications, the DRP property is more restrictive than stability criteria. Copyright © 2001 John Wiley & Sons, Ltd.     

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

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

Preconditioned finite element algorithms for 3D Stokes flows

Preconditioned conjugate gradient algorithms for solving 3D Stokes problems by stable piecewise discontinuous pressure finite elements are presented. The emphasis is on the preconditioning schemes and their numerical implementation for use with Hermitian based discontinuous pressure elements. For the piecewise constant discontinuous pressure elements, a variant implementation of the preconditioner proposed by Cahouet and Chabard for the continuous pressure elements is employed. For the piecewise linear discontinuous pressure elements, a new preconditioner is presented. Numerical examples are presented for the cubic lid-driven cavity problem with two representative elements, i.e. the Q2-PO and the Q2-P1 brick elements. Numerical results show that the preconditioning schemes are very effective in reducing the number of pressure iterations at very reasonable costs. It is also shown that they are insensitive to the mesh Reynolds number except for nearly steady flows (Re_m → 0) and are almost independent of mesh sizes. It is demonstrated that the schemes perform reasonably well on non-uniform meshes.     

Optimal shape design for the time‐dependent Navier–Stokes flow

This paper is concerned with the problem of shape optimization of two‐dimensional flows governed by the time‐dependent Navier–Stokes equations. We derive the structures of shape gradients for time‐dependent cost functionals by using the state derivative and its associated adjoint state. Finally, we apply a gradient‐type algorithm to our problem, and numerical examples show that our theory is useful for practical purposes and the proposed algorithm is feasible in low Reynolds number flows. Copyright © 2007 John Wiley & Sons, Ltd.     

Iterative solvers for quadratic discretizations of the generalized Stokes problem

We present in this paper various iterative methods for the solution of large linear and non‐linear systems resulting from the discretization of the generalized Stokes problem. A second‐order (O(h²)) P₂‐P₁ mixed finite element is used for the approximation of the velocity and the pressure. Solution strategies based on conjugate gradient‐like methods, the Uzawa's and Arrow–Hurwicz's methods are presented. Schur complement methods are also explored in the context of a hierarchical decomposition of the velocity field. The ever present preconditioning problem is also addressed. The performance of these iterative methods will be discussed on complex flows of industrial interest. Copyright © 2004 John Wiley & Sons, Ltd.     

Parallel three‐dimensional simulation of the injection molding process

This paper describes the development of a parallel three‐dimensional unstructured non‐isothermal flow solver for the simulation of the injection molding process. The numerical model accounts for multiphase flow in which the melt and air regions are considered to be a continuous incompressible fluid with distinct physical properties. This aspect avoids the complex reconstruction of the interface. A collocated finite volume method is employed, which can switch between first‐ and second‐order accuracy in both space and time. The pressure implicit with splitting of operators algorithm is used to compute the transient flow variables and couple velocity and pressure. The temperature equation is solved using a transport equation with convection and diffusion terms. An upwind differencing scheme is used for the discretization of the convection term to enforce a bounded solution. In order to capture the sharp interface, a bounded compressive high‐resolution scheme is employed. Parallelization of the code is achieved using the PETSc framework and a single program multiple data message passing model. Predicted numerical solutions for several example problems are considered. The first case validates the solution algorithm for moderate Reynolds number flows using a structured mesh. The second case employs an unstructured hybrid mesh showing the capability of the solver to describe highly viscous flows closer to realistic injection molding conditions. The final case presents the non‐isothermal filling of a thick cavity using three mesh sizes and up to 80 processors to assess parallel performance. The proposed algorithm is shown to have good accuracy and scalability. Copyright © 2008 John Wiley & Sons, Ltd.     

Boundary element analysis of driven cavity flow for low and moderate Reynolds numbers

A boundary element method for steady two‐dimensional low‐to‐moderate‐Reynolds number flows of incompressible fluids, using primitive variables, is presented. The velocity gradients in the Navier–Stokes equations are evaluated using the alternatives of upwind and central finite difference approximations, and derivatives of finite element shape functions. A direct iterative scheme is used to cope with the non‐linear character of the integral equations. In order to achieve convergence, an underrelaxation technique is employed at relatively high Reynolds numbers. Driven cavity flow in a square domain is considered to validate the proposed method by comparison with other published data. Copyright © 2001 John Wiley & Sons, Ltd.     

A spectral quadrilateral multidomain penalty method model for high Reynolds number incompressible stratified flows

A two‐dimensional quadrilateral spectral multidomain penalty method (SMPM) model has been developed for the simulation of high Reynolds number incompressible stratified flows. The implementation of higher‐order quadrilateral subdomains renders this model a nontrivial extension of a one‐dimensional subdomain SMPM model built for the simulation of the same type of flows in vertically nonperiodic domains (Diamessis et al., J. Comp. Phys, 202 :298‐322, 2005). The nontrivial aspects of this extension consist of the implementation of subdomain corners, the penalty formulation of the pressure Poisson equation (PPE), and, most importantly, the treatment of specific challenges that arise in the iterative solution of the SMPM‐discretized PPE. The two primary challenges within the framework of the iterative solution of the PPE are its regularization to ensure the consistency of the associated linear system of equations and the design of an appropriate two‐level preconditioner. A qualitative and quantitative assessment of the accuracy, efficiency, and stability of the quadrilateral SMPM solver is provided through its application to the standard benchmarks of the Taylor vortex, lid‐driven cavity, and double shear layer. The capacity of the flow solver for the study of environmental stratified flow processes is shown through the simulation of long‐distance propagation of an internal solitary wave of depression in a manner that is free of numerical dispersion and dissipation. The methods and results presented in this paper make it a point of reference for future studies oriented toward the reliable application of the quadrilateral SMPM model to more complex environmental stratified flow process studies. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical approximation of optimal control of unsteady flows using SQP and time decomposition

In this paper, we present numerical approximations of optimal control of unsteady flow problems using sequential quadratic programming method (SQP) and time domain decomposition. The SQP method is considered superior due to its fast convergence and its ability to take advantage of existing numerical techniques for fluid flow problems. It iteratively solves a sequence of linear quadratic optimal control problems converging to the solution of the non‐linear optimal control problem. The solution to the linear quadratic problem is characterized by the Karush–Kuhn–Tucker (KKT) optimality system which in the present context is a formidable system to solve. As a remedy various time domain decompositions, inexact SQP implementations and block iterative methods to solve the KKT systems are examined. Numerical results are presented showing the efficiency and feasibility of the algorithms. Copyright © 2004 John Wiley & Sons, Ltd.     

Computer‐algebra multiple‐timescale method for spatially periodic thin‐film viscous‐flow problems

A flexible, fully automated, computer‐algebra algorithm is developed for solving a class of non‐linear partial‐differential evolution equations arising frequently in the modeling of two‐dimensional transient free‐surface viscous thin‐film flows. The method, which is formulated for solving spatially periodic problems, is based upon an explicit multiple‐timescale asymptotic approximation of the thin‐film thickness. It admits the resolution of diverse physical phenomena by employing a finite geometric progression of increasingly slow timescales. The method is implemented on a challenging test problem comprising the evolution of an annular film of viscous liquid, with a free surface, adhering to the exterior of a horizontal rotating circular cylinder; as a model for numerous industrially motivated coating flows, this benchmark problem has been analyzed in diverse numerical and theoretical studies, against whose results those of the present method are compared. The explicit algebraic form of the solution admits a study of large‐time evolutionary dynamics that lies beyond the reach of considerably more expensive conventional numerical solvers, thereby shedding new light on the hitherto‐undiscovered explicit dependence of large‐time evolutionary fluid dynamics in terms of independent parameters describing gravitational and capillary effects. The results obtained from the new computer‐algebra procedure are demonstrated to be in good agreement with those obtained from a bespoke efficient numerical integration method that is spectrally accurate in space and 8th‐order (Runge–Kutta) in time. Newly discovered mechanisms describing the decay of free‐surface wave modes, from arbitrary initial conditions to the steady state, are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of preconditioned iterative solvers for incompressible flow problems

Solving efficiently the incompressible Navier–Stokes equations is a major challenge, especially in the three‐dimensional case. The approach investigated by Elman et al. (Finite Elements and Fast Iterative Solvers. Oxford University Press: Oxford, 2005) consists in applying a preconditioned GMRES method to the linearized problem at each iteration of a nonlinear scheme. The preconditioner is built as an approximation of an ideal block‐preconditioner that guarantees convergence in 2 or 3 iterations. In this paper, we investigate the numerical behavior for the three‐dimensional lid‐driven cavity problem with wedge elements; the ultimate motivation of this analysis is indeed the development of a preconditioned Krylov solver for stratified oceanic flows which can be efficiently tackled using such meshes. Numerical results for steady‐state solutions of both the Stokes and the Navier–Stokes problems are presented. Theoretical bounds on the spectrum and the rate of convergence appear to be in agreement with the numerical experiments. Sensitivity analysis on different aspects of the structure of the preconditioner and the block decomposition strategies are also discussed. Copyright © 2011 John Wiley & Sons, Ltd.