Galerkin domain decomposition procedures for parabolic equations on rectangular domain

Galerkin domain decomposition procedures for parabolic equations with three cases of boundary conditions on rectangular domain are discussed. These procedures are non‐iterative and non‐overlapping ones. They rely on implicit Galerkin method in the sub‐domains and integral mean method on the inter‐domain boundaries to present explicit flux calculation. Thus, the parallelism can be achieved by the use of these procedures. Two kinds of approximating schemes are presented. Because of the explicit nature of the flux calculation, a less severe time‐step constraint is derived to preserve stability. To bound L²‐norm error estimates, new elliptic projections are established and analyzed. Numerical experiments are provided to confirm theoretical results. Copyright © 2009 John Wiley & Sons, Ltd.     

A spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for conjugate heat transfer applications

We present a spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for fluid–solid conjugate heat transfer applications. Using the discrete Boltzmann equation, we propose a numerical scheme for conjugate heat transfer applications on unstructured, non‐uniform grids. We employ a double‐distribution thermal lattice Boltzmann model to resolve flows with variable Prandtl (Pr) number. Based upon its finite element heritage, the spectral‐element discontinuous Galerkin discretization provides an effective means to model and investigate thermal transport in applications with complex geometries. Our solutions are represented by the tensor product basis of the one‐dimensional Legendre–Lagrange interpolation polynomials. A high‐order discretization is employed on body‐conforming hexahedral elements with Gauss–Lobatto–Legendre quadrature nodes. Thermal and hydrodynamic bounce‐back boundary conditions are imposed via the numerical flux formulation that arises because of the discontinuous Galerkin approach. As a result, our scheme does not require tedious extrapolation at the boundaries, which may cause loss of mass conservation. We compare solutions of the proposed scheme with an analytical solution for a solid–solid conjugate heat transfer problem in a 2D annulus and illustrate the capture of temperature continuities across interfaces for conductivity ratio γ > 1. We also investigate the effect of Reynolds (Re) and Grashof (Gr) number on the conjugate heat transfer between a heat‐generating solid and a surrounding fluid. Steady‐state results are presented for Re = 5?40 and Gr = 10⁵?10⁶. In each case, we discuss the effect of Re and Gr on the heat flux (i.e. Nusselt number Nu) at the fluid–solid interface. Our results are validated against previous studies that employ finite‐difference and continuous spectral‐element methods to solve the Navier–Stokes equations. Copyright © 2016 John Wiley & Sons, Ltd.     

A new general-purpose least-squares finite element model for steady incompressible low-viscosity laminar flow using isoparametric C1-continuous Hermite elements

This paper deals with a critical evaluation of various finite element models for low-viscosity laminar incompressible flow in geometrically complex domains. These models use Galerkin weighted residuals UVP, continuous penalty, discrete penalty and least-squares procedures. The model evaluations are based on the use of appropriate tensor product Lagrange and simplex quadratic triangular elements and a newly developed isoparametric Hermite element. All of the described models produce very accurate results for horizontal flows. In vertical flow domains, however, two different cases can be recognized. Downward flows, i.e. when the gravitational force is in the direction of the flow, usually do not present any special problem. In contrast, laminar flow of low-viscosity Newtonian fluids where the gravitational force is acting in the direction opposite to the flow presents a difficult case. We show that only by using the least-squares method in conjunction with C¹-continuous Hermite elements can this type of laminar flow be modelled accurately. The problem of smooth isoparametric mapping of C¹ Hermite elements, which is necessary in dealing with geometrically complicated domains, is tackled by means of an auxiliary optimization procedure. We conclude that the least-squares method in combination with isoparmetric Hermite elements offers a new general-purpose modelling technique which can accurately simulate all types of low-viscosity incompressible laminar flow in complex domains.     

The locally conservative Galerkin (LCG) method for solving the incompressible Navier–Stokes equations

In this paper, the locally conservative Galerkin (LCG) method (Numer. Heat Transfer B Fundam. 2004; 46 :357–370; Int. J. Numer. Methods Eng. 2007) has been extended to solve the incompressible Navier–Stokes equations. A new correction term is also incorporated to make the formulation to give identical results to that of the continuous Galerkin (CG) method. In addition to ensuring element‐by‐element conservation, the method also allows solution of the governing equations over individual elements, independent of the neighbouring elements. This is achieved within the CG framework by breaking the domain into elemental sub‐domains. Although this allows discontinuous trial function field, we have carried out the formulation using the continuous trial function space as the basis. Thus, the changes in the existing CFD codes are kept to a minimum. The edge fluxes, establishing the continuity between neighbouring elements, are calculated via a post‐processing step during the time‐stepping operation. Therefore, the employed formulation needs to be carried out using either a time‐stepping or an equivalent iterative scheme that allows post‐processing of fluxes. The time‐stepping algorithm employed in this paper is based on the characteristic‐based split (CBS) scheme. Both steady‐ and unsteady‐state examples presented show that the element‐by‐element formulation employed is accurate and robust. Copyright © 2007 John Wiley & Sons, Ltd.     

An efficient high order direct ALE ADER finite volume scheme with a posteriori limiting for hydrodynamics and magnetohydrodynamics

In this paper, we present a new family of direct arbitrary–Lagrangian–Eulerian (ALE) finite volume schemes for the solution of hyperbolic balance laws on unstructured meshes in multiple space dimensions. The scheme is designed to be high‐order accurate both in space and time, and the mesh motion, which provides the new mesh configuration at the next time step, is taken into account in the final finite volume scheme that is based directly on a space‐time conservation formulation of the governing PDE system. To improve the computational efficiency of the algorithm, high order of accuracy in space is achieved using the a posteriori MOOD limiting strategy that allows the reconstruction procedure to be carried out with only one reconstruction stencil for any order of accuracy. We rely on an element‐local space‐time Galerkin finite element predictor on moving curved meshes to obtain a high‐order accurate one‐step time discretization, while the mesh velocity is computed by means of a suitable nodal solver algorithm that might also be supplemented with a local rezoning procedure to improve the mesh quality. Next, the old mesh configuration at time level tⁿ is connected to the new one at t^n + 1 by straight edges, hence providing unstructured space‐time control volumes, on the boundary of which the numerical flux has to be integrated. Here, we adopt a quadrature‐free integration, in which the space‐time boundaries of the control volumes are split into simplex sub‐elements that yield constant space‐time normal vectors and Jacobian matrices. In this way, the integrals over the simplex sub‐elements can be evaluated once and for all analytically during a preprocessing step. We apply the new high‐order direct ALE algorithm to the Euler equations of compressible gas dynamics (also referred to as hydrodynamics equations) as well as to the magnetohydrodynamics equations and we solve a set of classical test problems in two and three space dimensions. Numerical convergence rates are provided up to fifth order of accuracy in 2D and 3D for both hyperbolic systems considered in this paper. Finally, the efficiency of the new method is measured and carefully compared against the original formulation of the algorithm that makes use of a WENO reconstruction technique and Gaussian quadrature formulae for the flux integration: depending on the test problem, the new class of very efficient direct ALE schemes proposed in this paper can run up to ≈12 times faster in the 3D case. Copyright © 2016 John Wiley & Sons, Ltd.     

A numerical study of heat island flows: Stationary solutions

We present two‐dimensional numerical simulations of a natural convection problem in an unbounded domain. The flow circulation is induced by a heat island located on the ground and thermal stratification is applied in the vertical direction. The main effect of this stably stratified environment is to induce the propagation of thermal perturbations in the horizontal direction far from the local thermal source. Numerical stationary solutions at Ra?10⁵ are computed in large elongated computational domains: convergence with respect to the domain sizes is investigated at different resolutions. On fine grids, with mesh size , a thermal sponge layer is added at the vertical boundaries: this local damping technique improves the convergence with respect to the domain length. Boussinesq equations are discretized with a second‐order finite volume scheme on a staggered grid combined with a second‐order projection method for the time integration. Copyright © 2008 John Wiley & Sons, Ltd.     

流体饱和多孔介质动力问题的显式时域解法

为了简化分析,Zienkiewicz等基于Biot理论,在忽略流体相对于土骨架运动的加速度条件下,建立了以土骨架位移u和孔隙流体压力p为基本变量的u-p格式饱和两相介质动力方程。针对该u-p方程,在空间上,采用伽辽金法有限元离散,并结合对角化形式的质量矩阵和流体压缩矩阵,忽略相邻结点间的惯性和流体压缩量间的耦合作用。在时域内,基于杜修力等提出的显式算法和Euler预估-校正法,建立了一种具有二阶精度的全显式时域积分法。采用一维饱和土模型,对比提出算法的数值解与Simon方法的解析解,发现两者吻合良好,验证了本文方法的正确性。并分析了饱和土二维动力问题,以及渗透系数和排水条件对饱和土动力响应的影响。     

Characteristic Galerkin method for convection-diffusion equations and implicit algorithm using precise integration

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability. The project sponsored by the State Scientific and Technological Commission of China through “China State Key Project: the Theory and Methodology for Scientific and Engineering Computations with Large Scale”, the National Natural Science Foundation of China and the European Commission Research Project CI1^*CT94-0014.     

自然单元法计算裂纹与材料边界问题

提出一种新的非凸边界上自然单元法形函数计算方法,通过边界结点限制点对间的邻点关系,对包括裂纹和材料边界在内的各种类型的非凸边界具有统一的处理原则,所得到的近似函数在边界结点间具有线性插值性.     

Connections between the discontinuous Galerkin method and high‐order flux reconstruction schemes

With high‐order methods becoming more widely adopted throughout the field of computational fluid dynamics, the development of new computationally efficient algorithms has increased tremendously in recent years. One of the most recent methods to be developed is the flux reconstruction approach, which allows various well‐known high‐order schemes to be cast within a single unifying framework. Whilst a connection between flux reconstruction and the more widely adopted discontinuous Galerkin method has been established elsewhere, it still remains to fully investigate the explicit connections between the many popular variants of the discontinuous Galerkin method and the flux reconstruction approach. In this work, we closely examine the connections between three nodal versions of tensor‐product discontinuous Galerkin spectral element approximations and two types of flux reconstruction schemes for solving systems of conservation laws on quadrilateral meshes. The different types of discontinuous Galerkin approximations arise from the choice of the solution nodes of the Lagrange basis representing the solution and from the quadrature approximation used to integrate the mass matrix and the other terms of the discretization. By considering both linear and nonlinear advection equations on a regular grid, we examine the mathematical properties that connect these discretizations. These arguments are further confirmed by the results of an empirical numerical study. Copyright © 2014 John Wiley & Sons, Ltd.     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

Hybridizable discontinuous Galerkin p‐adaptivity for wave propagation problems

A p‐adaptive hybridizable discontinuous Galerkin method for the solution of wave problems is presented in a challenging engineering problem. Moreover, its performance is compared with a high‐order continuous Galerkin. The hybridization technique allows to reduce the coupled degrees of freedom to only those on the mesh element boundaries, whereas the particular choice of the numerical fluxes opens the path to a superconvergent postprocessed solution. This superconvergent postprocessed solution is used to construct a simple and inexpensive error estimator. The error estimator is employed to obtain solutions with the prescribed accuracy in the area (or areas) of interest and also drives a proposed iterative mesh adaptation procedure. The proposed method is applied to a nonhomogeneous scattering problem in an unbounded domain. This is a challenging problem because, on the one hand, for high frequencies, numerical difficulties are an important issue because of the loss of the ellipticity and the oscillatory behavior of the solution. And on the other hand, it is applied to real harbor agitation problems. That is, the mild slope equation in frequency domain (Helmholtz equation with nonconstant coefficients) is solved on real geometries with the corresponding perfectly matched layer to damp the diffracted waves. The performance of the method is studied on two practical examples. The adaptive hybridizable discontinuous Galerkin method exhibits better efficiency compared with a high‐order continuous Galerkin method using static condensation of the interior nodes. Copyright © 2013 John Wiley & Sons, Ltd.     

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

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

Time‐accurate intermediate boundary conditions for large eddy simulations based on projection methods

Projection methods are among the most adopted procedures for solving the Navier–Stokes equations system for incompressible flows. In order to simplify the numerical procedures, the pressure–velocity de‐coupling is often obtained by adopting a fractional time‐step method. In a specific formulation, suitable for the incompressible flows equations, it is based on a formal decomposition of the momentum equation, which is related to the Helmholtz–Hodge Decomposition theorem of a vector field in a finite domain. Owing to the continuity constraint also in large eddy simulation of turbulence, as happens for laminar solutions, the filtered pressure characterizes itself only as a Lagrange multiplier, not a thermodynamic state variable. The paper illustrates the implications of adopting such procedures when the decoupling is performed onto the filtered equations system. This task is particularly complicated by the discretization of the time integral of the sub‐grid scale tensor. A new proposal for developing time‐accurate and congruent intermediate boundary conditions is addressed. Several tests for periodic and non‐periodic channel flows are presented. This study follows and completes the previous ones reported in (Int. J. Numer. Methods Fluids 2003; 42, 43 ). Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Local discontinuous Galerkin method for a nonlocal viscous conservation laws

The purpose of this article is to investigate high‐order numerical approximations of scalar conservation laws with nonlocal viscous term. The viscous term is given in the form of convolution in space variable. With the help of the characteristic of viscous term, we design a semidiscrete local discontinuous Galerkin (LDG) method to solve the nonlocal model. We prove stability and convergence of semidiscrete LDG method in L² norm. The theoretical analysis reveals that the present numerical scheme is stable with optimal convergence order for the linear case, and it is stable with sub‐optimal convergence order for nonlinear case. To demonstrate the validity and accuracy of our scheme, we test the Burgers equation with two typical nonlocal fractional viscous terms. The numerical results show the convergence order accuracy in space for both linear and nonlinear cases. Some numerical simulations are provided to show the robustness and effectiveness of the present numerical scheme.     

An efficient and stable finite element solver of higher order in space and time for nonstationary incompressible flow

In this paper, we present fully implicit continuous Galerkin–Petrov (cGP) and discontinuous Galerkin (dG) time‐stepping schemes for incompressible flow problems which are, in contrast to standard approaches like for instance the Crank–Nicolson scheme, of higher order in time. In particular, we analyze numerically the higher order dG(1) and cGP(2) methods, which are super convergent of third, resp., fourth order in time, whereas for the space discretization, the well‐known LBB‐stable finite element pair of third‐order accuracy is used. The discretized systems of nonlinear equations are treated by using the Newton method, and the associated linear subproblems are solved by means of a monolithic (geometrical) multigrid method with a blockwise Vanka‐like smoother treating all components simultaneously. We perform nonstationary simulations (in 2D) for two benchmarking configurations to analyze the temporal accuracy and efficiency of the presented time discretization schemes w.r.t. CPU and numerical costs. As a first test problem, we consider a classical ‘flow around cylinder’ benchmark. Here, we concentrate on the nonstationary behavior of the flow patterns with periodic oscillations and examine the ability of the different time discretization schemes to capture the dynamics of the flow. As a second test case, we consider the nonstationary ‘flow through a Venturi pipe’. The objective of this simulation is to control the instantaneous and mean flux through this device. Copyright © 2013 John Wiley & Sons, Ltd.     

Resolution to a three‐dimensional physical oceanographic problem using the non‐linear Galerkin method

This paper presents the numerical results concerning the adaptation of the non‐linear Galerkin method to three‐dimensional geophysical fluid equations. This method was developed by Marion and Temam to solve the Navier–Stokes two‐dimensional equations. It allows a substantial decrease in calculation costs due to the application of an appropriate treatment to each mode based on its position in the spectrum. The large scales involved in the study of geophysical flow require that the earth's rotational effects and the existence of a high degree of stratification be taken into account. These phenomena play an important role in the distribution of the energy spectrum. It is shown here that the non‐linear Galerkin method is very well‐suited to the treatment of these phenomena. First, the method for the particular situation of a rigid‐lid with a flat bottom is validated, for which the functional basis used is particularly well‐adapted. Then the more general case of a domain exhibiting variable bathymetry is presented, which necessitates the use of the transformation σ, thus providing a study domain with a cylindrical configuration. Copyright © 1999 John Wiley & Sons, Ltd.     

Analysis of a discontinuous Galerkin approximation of the Stokes problem based on an artificial compressibility flux

In this work, we propose and analyse a discontinuous Galerkin (DG) method for the Stokes problem based on an artificial compressibility numerical flux. A crucial step in the definition of a DG method is the choice of the numerical fluxes, which affect both the accuracy and the order of convergence of the method. We propose here to treat the viscous and the inviscid terms separately. The former is discretized using the well‐known BRMPS method. For the latter, the problem is locally modified by adding an artificial compressibility term of the form (1/c²)(?p/?t) for the sole purpose of interface flux computation. The flux is obtained as the exact solution of a local Riemann problem. The analysis of the method extends the well‐established strategies for the DG discretization of the Laplacian to the resulting partially coercive problem. Copyright © 2007 John Wiley & Sons, Ltd.     

A numerical analysis of a class of generalized Boussinesq‐type equations using continuous/discontinuous FEM

An improved class of Boussinesq systems of an arbitrary order using a wave surface elevation and velocity potential formulation is derived. Dissipative effects and wave generation due to a time‐dependent varying seabed are included. Thus, high‐order source functions are considered. For the reduction of the system order and maintenance of some dispersive characteristics of the higher‐order models, an extra O(μ²ⁿ⁺²) term (n ∈ ?) is included in the velocity potential expansion. We introduce a nonlocal continuous/discontinuous Galerkin FEM with inner penalty terms to calculate the numerical solutions of the improved fourth‐order models. The discretization of the spatial variables is made using continuous P₂ Lagrange elements. A predictor‐corrector scheme with an initialization given by an explicit Runge–Kutta method is also used for the time‐variable integration. Moreover, a CFL‐type condition is deduced for the linear problem with a constant bathymetry. To demonstrate the applicability of the model, we considered several test cases. Improved stability is achieved. Copyright © 2011 John Wiley & Sons, Ltd.