A coupled continuous and discontinuous finite element method for the incompressible flows

In this paper, we develop a coupled continuous Galerkin and discontinuous Galerkin finite element method based on a split scheme to solve the incompressible Navier–Stokes equations. In order to use the equal order interpolation functions for velocity and pressure, we decouple the original Navier–Stokes equations and obtain three distinct equations through the split method, which are nonlinear hyperbolic, elliptic, and Helmholtz equations, respectively. The hybrid method combines the merits of discontinuous Galerkin (DG) and finite element method (FEM). Therefore, DG is concerned to accomplish the spatial discretization of the nonlinear hyperbolic equation to avoid using the stabilization approaches that appeared in FEM. Moreover, FEM is utilized to deal with the Poisson and Helmholtz equations to reduce the computational cost compared with DG. As for the temporal discretization, a second‐order stiffly stable approach is employed. Several typical benchmarks, namely, the Poiseuille flow, the backward‐facing step flow, and the flow around the cylinder with a wide range of Reynolds numbers, are considered to demonstrate and validate the feasibility, accuracy, and efficiency of this coupled method. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first‐order upwind approximation for the viscoelastic stress. A non‐uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non‐linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss–Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd‐B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright © 2001 John Wiley & Sons, Ltd.     

Spectral element method for viscoelastic flows in a planar contraction channel

A new algorithm, which combines the spectral element method with elastic viscous splitting stress (EVSS) method, has been developed for viscoelastic fluid flows in a planar contraction channel. The system of spectral element approximations to the velocity, pressure, extra stress and the rate of deformation variables is solved by a preconditioned conjugate gradient method based on the Uzawa iteration procedure. The numerical approach is implemented on a planar four‐to‐one contraction channel for a fluid governed by an Oldroyd‐B constitutive equation. The behaviour of the Oldroyd‐B fluids in the contraction channel is investigated with various Weissenberg numbers. It is shown that numerical solutions obtained here agree well with experimental measurements and other numerical predictions. Copyright © 2003 John Wiley & Sons, Ltd.     

A finite volume approach for unsteady viscoelastic fluid flows

A finite volume, time‐marching for solving time‐dependent viscoelastic flow in two space dimensions for Oldroyd‐B and Phan Thien–Tanner fluids, is presented. A non‐uniform staggered grid system is used. The conservation and constitutive equations are solved using the finite volume method with an upwind scheme for the viscoelastic stresses and an hybrid scheme for the velocities. To calculate the pressure field, the semi‐implicit method for the pressure linked equation revised method is used. The discretized equations are solved sequentially, using the tridiagonal matrix algorithm solver with under‐relaxation. In both, the full approximation storage multigrid algorithm is used to speed up the convergence rate. Simulations of viscoelastic flows in four‐to‐one abrupt plane contraction are carried out. We will study the behaviour at the entrance corner of the four‐to‐one planar abrupt contraction. Using this solver, we show convergence up to a Weissenberg number We of 20 for the Oldroyd‐B model. No limiting Weissenberg number is observed even though a Phan Thien–Tanner model is used. Several numerical results are presented. Smooth and stable solutions are obtained for high Weissenberg number. Copyright © 2002 John Wiley & Sons, Ltd.     

Finite element method for viscoelastic flows based on the discrete adaptive viscoelastic stress splitting and the discontinuous Galerkin method: DAVSS-G/DG

Accurate and robust finite element methods for computing flows with differential constitutive equations require approximation methods that numerically preserve the ellipticity of the saddle point problem formed by the momentum and continuity equations and give numerically stable and accurate solutions to the hyperbolic constitutive equation. We present a new finite element formulation based on the synthesis of three ideas: the discrete adaptive splitting method for preserving the ellipticity of the momentum/continuity pair (the DAVSS formulation), independent interpolation of the components of the velocity gradient tensor (DAVSS-G), and application of the discontinuous Galerkin (DG) method for solving the constitutive equation. We call the method DAVSS-G/DG. The DAVSS-G/DG method is compared with several other methods for flow past a cylinder in a channel with the Oldroyd-B and Giesekus constitutive models. Results using the Streamline Upwind Petrov–Galerkin method (SUPG) show that introducing the adaptive splitting increases considerably the range of Deborah number (De) for convergence of the calculations over the well established EVSS-G formulation. When both formulations converge, the DAVSS-G and DEVSS-G methods give comparable results. Introducing the DG method for solution of the constitutive equation extends further the region of convergence without sacrificing accuracy. Calculations with the Oldroyd-B model are only limited by approximation of the almost singular gradients of the axial normal stress that develop near the rear stagnation point on the cylinder. These gradients are reduced in calculations with the Giesekus model. Calculations using the Giesekus model with the DAVSS-G/DG method can be continued to extremely large De and converge with mesh refinement.     

Time‐dependent finite element simulations of a shear‐thinning viscoelastic fluid with application to blood flow

A new finite element method is developed to simulate time‐dependent viscoelastic shear‐thinning flows characterized by the generalized Oldroyd‐B model. The focus of the algorithm is improved stability through a free‐energy dissipative scheme by using low‐order piecewise‐constant finite element approximations for stress. The algorithm is further modified by incorporating a pressure‐projection method, a DG‐upwinding scheme, a symmetric interior penalty DG method to solve the elliptic pressure‐update equation and a geometric multigrid preconditioner. The improved stability and cost to accuracy is compared when using higher order discontinuous bilinear approximation, where in addition, we consider the influence of a slope limiter for these elements. The algorithm is applied to the 2D start‐up‐driven cavity problem, and the stability of the free energy is illustrated and compared between element choices. An application of the model to modelling blood in small arterioles and channels is considered by simulating pulsatile blood flow through a stenotic arteriole. The individual influences of viscoelasticity and shear‐thinning within the generalized Oldroyd‐B model are investigated by comparing results to the Newtonian, generalized Newtonian and Oldroyd‐B models. Copyright © 2014 John Wiley & Sons, Ltd.     

On the treatment of transient area variation in 1D discontinuous Galerkin simulations of train‐induced pressure waves in tunnels

To simulate the pressure wave generated by a train travelling through a tunnel, we implement a discontinuous Galerkin (DG) method for the solution of the one‐dimensional equations of variable area flow. This formulation uses a spatial discretisation via Legendre polynomials of arbitrary degree, and the resulting semi‐discrete system is integrated using an explicit Runge–Kutta scheme. A simulation of subsonic steady flow in a nozzle shows that the scheme produces stable solutions, without the need for artificial dissipation, and that its performance is optimal for polynomial degrees between 5 and 7. However, when dealing with an unsteady area, we report the presence of numerical oscillations that are not due to the steep pressure fronts in the flow but rather to the projection of a moving area, with piecewise continuous derivatives onto a fixed grid. We propose a reformulation of the DG method to eliminate these oscillations that, put in simple terms, amount to splitting the integrals where the derivatives of the cross‐sectional area are discontinuous into subintegrals where they are continuous. The resulting method does not exhibit oscillations, and it is applied here to two practical cases involving train‐induced pressure waves in a tunnel. The first application is a validation of the DG method through comparison of its computational results with pressure data measured during transit at the Patchway tunnel near Bristol (UK). The second application is a study of the influence of the nose shape and length on the pressure wave gradients responsible for sonic boom at tunnel exit portals to show that the proposed modification is able to deal with realistic train shapes. Copyright © 2012 John Wiley & Sons, Ltd.     

High-order hybridisable discontinuous Galerkin method for the gas kinetic equation

ABSTRACT

The high-order hybridisable discontinuous Galerkin (HDG) method is used to find steady-state solution of gas kinetic equations on two-dimensional geometry. The velocity distribution function and its traces are approximated in piecewise polynomial space on triangular mesh and mesh skeleton, respectively. By employing a numerical flux derived from the upwind scheme and imposing its continuity on mesh skeleton, the global system for unknown traces is obtained with fewer coupled degrees of freedom, compared to the original DG method. The solutions of model equation for the Poiseuille flow through square channel show the higher order solver is faster than the lower order one. Moreover, the HDG scheme is more efficient than the original DG method when the degree of approximating polynomial is larger than 2. Finally, the developed scheme is extended to solve the Boltzmann equation with full collision operator, which can produce accurate results for shear-driven and thermally induced flows.     

A high‐order discontinuous Galerkin method for all‐speed flows

In this article, we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of steady solutions of the compressible fully coupled Reynolds‐averaged Navier–Stokes and k ? ω turbulence model equations for solving all‐speed flows. The system of equations is iterated to steady state by means of an implicit scheme. The DG solution is extended to the incompressible limit by implementing a low Mach number preconditioning technique. A full preconditioning approach is adopted, which modifies both the unsteady terms of the governing equations and the dissipative term of the numerical flux function by means of a new preconditioner, on the basis of a modified version of Turkel's preconditioning matrix. At sonic speed the preconditioner reduces to the identity matrix thus recovering the non‐preconditioned DG discretization. An artificial viscosity term is added to the DG discretized equations to stabilize the solution in the presence of shocks when piecewise approximations of order of accuracy higher than 1 are used. Moreover, several rescaling techniques are implemented in order to overcome ill‐conditioning problems that, in addition to the low Mach number stiffness, can limit the performance of the flow solver. These approaches, through a proper manipulation of the governing equations, reduce unbalances between residuals as a result of the dependence on the size of elements in the computational mesh and because of the inherent differences between turbulent and mean‐flow variables, influencing both the evolution of the Courant Friedrichs Lewy (CFL) number and the inexact solution of the linear systems. The performance of the method is demonstrated by solving three turbulent aerodynamic test cases: the flat plate, the L1T2 high‐lift configuration and the RAE2822 airfoil (Case 9). The computations are performed at different Mach numbers using various degrees of polynomial approximations to analyze the influence of the proposed numerical strategies on the accuracy, efficiency and robustness of a high‐order DG solver at different flow regimes. Copyright © 2014 John Wiley & Sons, Ltd.     

Stability/dispersion analysis of the discontinuous Galerkin linearized shallow‐water system

The frequency or dispersion relation for the discontinuous Galerkin mixed formulation of the 1‐D linearized shallow‐water equations is analysed, using several basic DG mixed schemes. The dispersion properties are compared analytically and graphically with those of the mixed continuous Galerkin formulation for piecewise‐linear bases on co‐located grids. Unlike the Galerkin case, the DG scheme does not exhibit spurious stationary pressure modes. However, spurious propagating modes have been identified in all the present discontinuous Galerkin formulations. Numerical solutions of a test problem to simulate fast gravity modes illustrate the theoretical results and confirm the presence of spurious propagating modes in the DG schemes. Copyright © 2005 John Wiley & Sons, Ltd.     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 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.     

Continuous adjoint-based error estimation and its application to adaptive discontinuous Galerkin method

An adaptive mesh refinement algorithm based on a continuous adjoint approach is developed. Both the primal equation and the adjoint equation are approximated with the discontinuous Galerkin (DG) method. The proposed adaptive algorithm is used in compressible Euler equations. Numerical tests are made to show the superiority of the proposed adaptive algorithm.     

Viscoelastic flow past confined objects using a micro–macro approach

This paper is concerned with the numerical prediction of viscoelastic flow past a cylinder in a channel and a sphere in a cylinder using molecular-based models. The basis of the numerical method employed is a micro–macro model in which the polymer dynamics is described by the evolution of an ensemble of Brownian configuration fields. The spectral element method is used to discretize the equations in space. Comparisons are made between the macroscopic simulations based on the Oldroyd B constitutive model and microscopic simulations based on Hookean dumbbells, and excellent agreement is found. The micro–macro approach can be used to simulate models, such as the finitely extensible nonlinear elastic (FENE) dumbbell model, which do not possess a closed-form constitutive equation. Numerical simulations are performed for the FENE model. The influence of the model parameters on the flow is described and, in particular, the dependence of the drag as a function of the Weissenberg number.     

A γ‐DGBGK scheme for compressible multi‐fluids

The paper presents a Discontinuous Galerkin γ‐BGK (γ‐DGBGK) method for compressible multicomponent flow simulations by coupling the discontinuous Galerkin method with a γ‐BGK scheme based on WENO limiters. In this γ‐DGBGK method, the construction of the flux in the DG method is based on the kinetic scheme which not only couples the convective and dissipative terms together, but also includes both discontinuous and continuous terms in the flux formulation at cell interfaces. WENO limiters are used to obtain uniform high‐order accuracy and sharp non‐oscillatory shock transition, and time accuracy obtained by integration for the flux function at the cell interface. Numerical examples in one and two space dimensions are presented to illustrate the robust and accuracy of the present scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

Discontinuous Galerkin finite element method applied to the coupled unsteady Stokes/Cahn-Hilliard equations

Two-phase flows driven by the interfacial dynamics are studied by tracking implicitly interfaces in the framework of the Cahn-Hilliard theory. The fluid dynamics is described by the Stokes equations with an additional source term in the momentum equation taking into account the capillary forces. A discontinuous Galerkin finite element method is used to solve the coupled Stokes/Cahn-Hilliard equations. The Cahn-Hilliard equation is treated as a system of two coupled equations corresponding to the advection-diffusion equation for the phase field and a nonlinear elliptic equation for the chemical potential. First, the variational formulation of the Cahn-Hilliard equation is presented. A numerical test is achieved showing the optimal order in error bounds. Second, the variational formulation in discontinuous Galerkin finite element approach of the Stokes equations is recalled, in which the same space of approximation is used for the velocity and the pressure with an adequate stabilization technique. The rates of convergence in space and time are evaluated leading to an optimal order in error bounds in space and a second order in time with a backward differentiation formula at the second order. Numerical tests devoted to two-phase flows are provided on ellipsoidal droplet retraction, on the capillary rising of a liquid in a tube, and on the wetting drop over a horizontal solid wall.     

A high‐order accurate discontinuous Galerkin finite element method for laminar low Mach number flows

In this paper we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of laminar flow simulations at low Mach numbers using an implicit scheme. The algorithm is based on the flux preconditioning approach, which modifies only the dissipative terms of the numerical flux. This formulation is quite simple to implement in existing implicit DG codes, it overcomes the time‐stepping restrictions of explicit multistage algorithms, is consistent in time and thus applicable to unsteady flows. The performance of the method is demonstrated by solving the flow around a NACA0012 airfoil and on a flat plate, at different low Mach numbers using various degrees of polynomial approximations. Computations with and without flux preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers. The time accurate solution of unsteady flow is also demonstrated by solving the vortex shedding behind a circular cylinder at the Reynolds number of 100. Copyright © 2012 John Wiley & Sons, Ltd.     

A kinetic energy preserving nodal discontinuous Galerkin spectral element method

In this work, we discuss the construction of a skew‐symmetric discontinuous Galerkin (DG) collocation spectral element approximation for the compressible Euler equations. Starting from the skew‐symmetric formulation of Morinishi, we mimic the continuous derivations on a discrete level to find a formulation for the conserved variables. In contrast to finite difference methods, DG formulations naturally have inter‐domain surface flux contributions due to the discontinuous nature of the approximation space. Thus, throughout the derivations we accurately track the influence of the surface fluxes to arrive at a consistent formulation also for the surface terms. The resulting novel skew‐symmetric method differs from the standard DG scheme by additional volume terms. Those volume terms have a special structure and basically represent the discretization error of the different product rules. We use the summation‐by‐parts (SBP) property of the Gauss–Lobatto‐based DG operator and show that the novel formulation is exactly conservative for the mass, momentum, and energy. Finally, an analysis of the kinetic energy balance of the standard DG discretization shows that because of aliasing errors, a nonzero transport source term in the evolution of the discrete kinetic energy mean value may lead to an inconsistent increase or decrease in contrast to the skew‐symmetric formulation. Furthermore, we derive a suitable interface flux that guarantees kinetic energy preservation in combination with the skew‐symmetric DG formulation. As all derivations require only the SBP property of the Gauss–Lobatto‐based DG collocation spectral element method operator and that the mass matrix is diagonal, all results for the surface terms can be directly applied in the context of multi‐domain diagonal norm SBP finite difference methods. Numerical experiments are conducted to demonstrate the theoretical findings. Copyright © 2014 John Wiley & Sons, Ltd.     

Implicit high-order discontinuous Galerkin method with HWENO type limiters for steady viscous flow simulations

Two types of implicit algorithms have been improved for high order discontinuous Galerkin (DG) method to solve compressible Navier-Stokes (NS) equations on triangular grids. A block lower-upper symmetric Gauss-Seidel (BLU-SGS) approach is implemented as a nonlinear iterative scheme. And a modified LU-SGS (LLU-SGS) approach is suggested to reduce the memory requirements while retain the good convergence performance of the original LU-SGS approach. Both implicit schemes have the significant advantage that only the diagonal block matrix is stored. The resulting implicit high-order DG methods are applied, in combination with Hermite weighted essentially non-oscillatory (HWENO) limiters, to solve viscous flow problems. Numerical results demonstrate that the present implicit methods are able to achieve significant efficiency improvements over explicit counterparts and for viscous flows with shocks, and the HWENO limiters can be used to achieve the desired essentially non-oscillatory shock transition and the designed high-order accuracy simultaneously.     

Well‐balancing issues related to the RKDG2 scheme for the shallow water equations

Discontinuous Galerkin (DG) finite element methods have salient features that are mainly highlighted by their locality, their easiness in balancing the flux and source term gradients and their component‐wise structure. In the light of this, this paper aims to provide insights into the well‐balancing property of a second‐order Runge–Kutta Discontinuous Galerkin (RKDG2) method. For this purpose, a Godunov‐type RKDG2 method is presented for solving the shallow water equations. The scheme is based on local DG linear approximations and does not entail any special treatment of the source terms in order to achieve well‐balanced numerical results. The performance of the present RKDG2 scheme in reproducing conserved solutions for both free surface and discharge over strongly irregular topography is demonstrated by applying to several hydraulic benchmarks. Meanwhile, the effects of different slope limiting procedures on the well‐balancing property are investigated and discussed. This work may provide useful guidelines for developing a well‐balanced RKDG2 numerical scheme for shallow water flow simulation. Copyright © 2009 John Wiley & Sons, Ltd.