Numerical modelling of shear-dependent mass transfer in large arteries

A numerical scheme for the simulation of blood flow and transport processes in large arteries is presented. Blood flow is described by the unsteady 3D incompressible Navier–Stokes equations for Newtonian fluids; solute transport is modelled by the advection–diffusion equation. The resistance of the arterial wall to transmural transport is described by a shear-dependent wall permeability model. The finite element formulation of the Navier–Stokes equations is based on an operator-splitting method and implicit time discretization. The streamline upwind/Petrov–Galerkin (SUPG) method is applied for stabilization of the advective terms in the transport equation and in the flow equations. A numerical simulation is carried out for pulsatile mass transport in a 3D arterial bend to demonstrate the influence of arterial flow patterns on wall permeability characteristics and transmural mass transfer. The main result is a substantial wall flux reduction at the inner side of the curved region. © 1997 John Wiley & Sons, Ltd.     

SUPG finite element method based on penalty function for lid-driven cavity flow up to $$Re = 27500$$

A streamline upwind/Petrov–Galerkin(SUPG)finite element method based on a penalty function is proposed for steady incompressible Navier–Stokes equations.The SUPG stabilization technique is employed for the formulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pressure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.     

A two‐phase flow interface capturing finite element method

A new interface capturing algorithm is proposed for the finite element simulation of two‐phase flows. It relies on the solution of an advection equation for the interface between the two phases by a streamline upwind Petrov–Galerkin (SUPG) scheme combined with an adaptive mesh refinement procedure and a filtering technique. This method is illustrated in the case of a Rayleigh–Taylor two‐phase flow problem governed by the Stokes equations. Copyright © 2006 John Wiley & Sons, Ltd.     

Time‐accurate stabilized finite‐element model for weakly nonlinear and weakly dispersive water waves

Introduction of a time‐accurate stabilized finite‐element approximation for the numerical investigation of weakly nonlinear and weakly dispersive water waves is presented in this paper. To make the time approximation match the order of accuracy of the spatial representation of the linear triangular elements by the Galerkin finite‐element method, the fourth‐order time integration of implicit multistage Padé method is used for the development of the numerical scheme. The streamline‐upwind Petrov–Galerkin (SUPG) method with crosswind diffusion is employed to stabilize the scheme and suppress the spurious oscillations, usually common in the numerical computation of convection‐dominated flow problems. The performance of numerical stabilization and accuracy is addressed. Treatments of various boundary conditions, including the open boundary conditions, the perfect reflecting boundary conditions along boundaries with irregular geometry, are also described. Numerical results showing the comparisons with analytical solutions, experimental measurements, and other published numerical results are presented and discussed. Copyright © 2007 John Wiley & Sons, Ltd.     

A stabilized formulation for the advection–diffusion equation using the Generalized Finite Element Method

This paper presents a stable formulation for the advection–diffusion equation based on the Generalized (or eXtended) Finite Element Method, GFEM (or X‐FEM). Using enrichment functions that represent the exponential character of the exact solution, smooth numerical solutions are obtained for problems with steep gradients and high Péclet numbers in one‐ and two‐dimensions. In contrast with traditional stabilized methods that require the construction of stability parameters and stabilization terms, the present work avoids numerical instabilities by improving the classical Galerkin solution with enrichment functions (that need not be polynomials) using GFEM, which is an instance of the partition of unity framework. This work also presents a strategy for constructing enrichment functions for problems involving complex geometries by employing a global–local‐type approach. Representative numerical results are presented to illustrate the performance of the proposed method. Copyright © 2010 John Wiley & Sons, Ltd.     

Large‐eddy simulation with complex 2‐D geometries using a parallel finite‐element/spectral algorithm

A parallel stabilized finite‐element/spectral formulation is presented for incompressible large‐eddy simulation with complex 2‐D geometries. A unique discretization scheme is developed consisting of a streamline‐upwind Petrov–Galerkin/Pressure‐Stabilized Petrov–Galerkin (SUPG/PSPG) finite‐element discretization in the 2‐D plane with a collocated spectral/pseudospectral discretization in the out‐of‐plane direction. This formulation provides an efficient approach for solving 3‐D flows over arbitrary 2‐D geometries. Utilizing this discretization and through explicit temporal treatment of the non‐linear terms, the system of equations for each Fourier mode is decoupled within each time step. A novel parallelization approach is then taken, where the computational work is partitioned in Fourier space. A validation of the algorithm is presented via comparison of results for flow past a circular cylinder with published values for Re=195, 300, and 3900. Copyright © 2003 John Wiley & Sons, Ltd.     

Analysis of flow in the plate‐spiral of a reaction turbine using a streamline upwind Petrov–Galerkin method

The prediction of the flow field in a novel spiral casing has been accomplished. Hydraulic turbine manufacturers are considering the potential of using a special type of spiral casing because of the easier manufacturing process involved in its fabrication. These special spiral casings are known as plate‐spirals. Numerical simulation of complex three‐dimensional flow through such spiral casings has been accomplished using a finite element method (FEM). An explicit Eulerian velocity correction scheme has been deployed to solve the Reynolds‐average Navier–Stokes equations. The simulation has been performed to describe the flow in high Reynolds number (10⁶) regimes. For spatial discretization, a streamline upwind Petrov–Galerkin (SUPG) technique has been used. The velocity field and the pressure distribution inside the spiral casing reveal meaningful results. Copyright © 2000 John Wiley & Sons, Ltd.     

Generalized Fourier analyses of the advection–diffusion equation—Part I: one‐dimensional domains

This paper presents a detailed multi‐methods comparison of the spatial errors associated with finite difference, finite element and finite volume semi‐discretizations of the scalar advection–diffusion equation. The errors are reported in terms of non‐dimensional phase and group speed, discrete diffusivity, artificial diffusivity, and grid‐induced anisotropy. It is demonstrated that Fourier analysis provides an automatic process for separating the discrete advective operator into its symmetric and skew‐symmetric components and characterizing the spectral behaviour of each operator. For each of the numerical methods considered, asymptotic truncation error and resolution estimates are presented for the limiting cases of pure advection and pure diffusion. It is demonstrated that streamline upwind Petrov–Galerkin and its control‐volume finite element analogue, the streamline upwind control‐volume method, produce both an artificial diffusivity and a concomitant phase speed adjustment in addition to the usual semi‐discrete artifacts observed in the phase speed, group speed and diffusivity. The Galerkin finite element method and its streamline upwind derivatives are shown to exhibit super‐convergent behaviour in terms of phase and group speed when a consistent mass matrix is used in the formulation. In contrast, the CVFEM method and its streamline upwind derivatives yield strictly second‐order behaviour. In Part II of this paper, we consider two‐dimensional semi‐discretizations of the advection–diffusion equation and also assess the affects of grid‐induced anisotropy observed in the non‐dimensional phase speed, and the discrete and artificial diffusivities. Although this work can only be considered a first step in a comprehensive multi‐methods analysis and comparison, it serves to identify some of the relative strengths and weaknesses of multiple numerical methods in a common analysis framework. Published in 2004 by John Wiley & Sons, Ltd.     

Development and validation of a SUPG finite element scheme for the compressible Navier–Stokes equations using a modified inviscid flux discretization

This paper considers the streamline‐upwind Petrov–Galerkin (SUPG) method applied to the unsteady compressible Navier–Stokes equations in conservation‐variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, and the second‐order accurate time discretization are presented. The numerical method is discussed in detail. The performance of the algorithm is then investigated by considering inviscid flow past a circular cylinder. Validation of the finite element formulation via comparisons with experimental data for high‐Mach number perfect gas laminar flows is presented, with a specific focus on comparisons with experimentally measured skin friction and convective heat transfer on a 15° compression ramp. Copyright © 2007 John Wiley & Sons, Ltd.     

Comparison of Galerkin and control volume finite element for advection–diffusion problems

The control volume finite element method (CVFEM) was developed to combine the local numerical conservation property of control volume methods with the unstructured grid and generality of finite element methods (FEMs). Most implementations of CVFEM include mass‐lumping and upwinding techniques typical of control volume schemes. In this work we compare, via numerical error analysis, CVFEM and FEM utilizing consistent and lumped mass implementations, and stabilized Petrov–Galerkin streamline upwind schemes in the context of advection–diffusion processes. For this type of problem, we find no apparent advantage to the local numerical conservation aspect of CVFEM as compared to FEM. The stabilized schemes improve accuracy and degree of positivity on coarse grids, and also reduce iteration counts for advection‐dominated problems. Published in 2005 by John Wiley & Sons, Ltd.     

3D macrosegregation simulation with anisotropic remeshing

The article presents a three-dimensional coupled numerical solution of momentum, mass, energy and solute conservation equations, for binary alloy solidification. The interdendritic flow in the mushy zone is assumed to obey the Darcy's law. Microsegregation is governed by the lever rule, assuming local equilibrium at phase interfaces. The resulting energy and solute advection–diffusion equations are solved using the Streamline-Upwind/Petrov–Galerkin (SUPG) finite element method. A SUPG-PSPG velocity-pressure formulation is applied for the momentum equation. The full algorithm was implemented in the 3D code THERCAST, together with an anisotropic remeshing method. Two applications have been considered: a small ingot of Pb-48wt%Sn alloy and a large steel ingot. The numerical results of these two cases are presented with the evolution of temperature, liquid velocity, and solute concentration fields during solidification. To cite this article: S. Gouttebroze et al., C. R. Mecanique 335 (2007).     

A least‐square weighted residual method for level set formulation

In this paper, a least‐square weighted residual method (LSWRM) for level set (LS) formulation is introduced to achieve interface capturing in two‐dimensional (2D) and three‐dimensional (3D) problems. An LSWRM was adopted for two semi‐discretized advection and reinitialization equations of the LS formulation. The present LSWRM provided good mathematical properties such as natural numerical diffusion and the symmetry of the resulting algebraic systems for the advection and reinitialization equations. The proposed method was validated by solving some 2D and 3D benchmark problems such as those involving a rotating slotted disk, the rotation of a slotted sphere, and a time‐reversed single‐vortex flow and a deformation problem of a spherical fluid. The numerical results were compared with those obtained from essentially non‐oscillatory type formulations and particle LS methods. Further, the proposed LSWRM for the LS formulation was coupled with a splitting finite element method code to solve the incompressible Navier–Stokes equations, and then, the collapse of a 3D broken dam flow was well simulated; in the simulation, the entrapping of air and the splashing of the surge front of water were reproduced. The mass conservation of the present method was found to be satisfactory during the entire simulation. Copyright © 2011 John Wiley & Sons, Ltd.     

AN ENTROPY VARIABLE FORMULATION AND APPLICATIONS FOR THE TWO-DIMENSIONAL SHALLOW WATER EQUATIONS

A new symmetric formulation of the two-dimensional shallow water equations and a streamline upwind Petrov–Galerkin (SUPG) scheme are developed and tested. The symmetric formulation is constructed by means of a transformation of dependent variables derived from the relation for the total energy of the water column. This symmetric form is well suited to the SUPG approach as seen in analogous treatments of gas dynamics problems based on entropy variables. Particulars related to the construction of the upwind test functions and an appropriate discontinuity-capturing operator are included. A formal extension to the viscous, dissipative problem and a stability analysis are also presented. Numerical results for shallow water flow in a channel with (a) a step transition, (b) a curved wall transition and (c) a straight wall transition are compared with experimental and other computational results from the literature.     

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.     

Finite element solution of three‐dimensional turbulent flows applied to mold‐filling problems

This paper presents a finite element solution algorithm for three‐dimensional isothermal turbulent flows for mold‐filling applications. The problems of interest present unusual challenges for both the physical modelling and the solution algorithm. High‐Reynolds number transient turbulent flows with free surfaces have to be computed on complex three‐dimensional geometries. In this work, a segregated algorithm is used to solve the Navier–Stokes, turbulence and front‐tracking equations. The streamline–upwind/Petrov–Galerkin method is used to obtain stable solutions to convection‐dominated problems. Turbulence is modelled using either a one‐equation turbulence model or the κ–ε two‐equation model with wall functions. Turbulence equations are solved for the natural logarithm of the turbulence variables. The change of dependent variables allows for a robust solution algorithm and good predictions even on coarse meshes. This is very important in the case of large three‐dimensional applications for which highly refined meshes result in untreatable large numbers of elements. The position of the flow front in the mold cavity is computed using a level set approach. Finally, equations are integrated in time using an implicit Euler scheme. The methodology presents the robustness and cost effectiveness needed to tackle complex industrial applications. Copyright © 2000 John Wiley & Sons, Ltd.     

Modified Galerkin Method for Unsteady Two-Dimensional Viscous Fluid Flow Analysis

The Modified Galerkin Method (MGM) has been proposed as one of the most efficient methods for two-dimensional convection-diffusion equations. In the MGM, the non-symmetric matrices, which are derived from the convection term in the Galerkin formulation, are not used, and an artificial diffusion is introduced through an error analysis approach to improve its discretization accuracy in both time and space directions. In this study, the MGM is applied for two-dimensional viscous fluid flow analysis, and the driven cavity flow problems are solved up to Reynolds number of 10,000 using the vorticity-stream function formulation and non-uniform meshes. The results show the effectiveness of MGM.     

Adiabatic shock capturing in perfect gas hypersonic flows

This paper considers the streamline‐upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier–Stokes equations in conservation‐variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock‐capturing operator, which is required to regularize the scheme in the presence of strong, shock‐induced gradients. A standard shock‐capturing operator that has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator that is designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite‐dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock‐capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator. Published in 2009 by John Wiley & Sons, Ltd.     

A multi‐dimensional monotonic finite element model for solving the convection–diffusion‐reaction equation

In this paper, we develop a finite element model for solving the convection–diffusion‐reaction equation in two dimensions with an aim to enhance the scheme stability without compromising consistency. Reducing errors of false diffusion type is achieved by adding an artificial term to get rid of three leading mixed derivative terms in the Petrov–Galerkin formulation. The finite element model of the Petrov–Galerkin type, while maintaining convective stability, is modified to suppress oscillations about the sharp layer by employing the M‐matrix theory. To validate this monotonic model, we consider test problems which are amenable to analytic solutions. Good agreement is obtained with both one‐ and two‐dimensional problems, thus validating the method. Other problems suitable for benchmarking the proposed model are also investigated. Copyright © 2002 John Wiley & Sons, Ltd.     

Higher‐order bounded differencing schemes for compressible and incompressible flows

In recent years, three higher‐order (HO) bounded differencing schemes, namely AVLSMART, CUBISTA and HOAB that were derived by adopting the normalized variable formulation (NVF), have been proposed. In this paper, a comparative study is performed on these schemes to assess their numerical accuracy, computational cost as well as iterative convergence property. All the schemes are formulated on the basis of a new dual‐formulation in order to facilitate their implementations on unstructured meshes. Based on the proposed dual‐formulation, the net effective blending factor (NEBF) of a high‐resolution (HR) scheme can now be measured and its relevance on the accuracy and computational cost of a HR scheme is revealed on three test problems: (1) advection of a scalar step‐profile; (2) 2D transonic flow past a circular arc bump; and (3) 3D lid‐driven incompressible cavity flow. Both density‐based and pressure‐based methods are used for the computations of compressible and incompressible flow, respectively. Computed results show that all the schemes produce solutions which are nearly as accurate as the third‐order QUICK scheme; however, without the unphysical oscillations which are commonly inherited from the HO linear differencing scheme. Generally, it is shown that at higher value of NEBF, a HR scheme can attain better accuracy at the expense of computational cost. Copyright © 2006 John Wiley & Sons, Ltd.