A Petrov–Galerkin finite element model for one‐dimensional fully non‐linear and weakly dispersive wave propagation

A new finite element method is presented to solve one‐dimensional depth‐integrated equations for fully non‐linear and weakly dispersive waves. For spatial integration, the Petrov–Galerkin weighted residual method is used. The weak forms of the governing equations are arranged in such a way that the shape functions can be piecewise linear, while the weighting functions are piecewise cubic with C²‐continuity. For the time integration an implicit predictor–corrector iterative scheme is employed. Within the framework of linear theory, the accuracy of the scheme is discussed by considering the truncation error at a node. The leading truncation error is fourth‐order in terms of element size. Numerical stability of the scheme is also investigated. If the Courant number is less than 0.5, the scheme is unconditionally stable. By increasing the number of iterations and/or decreasing the element size, the stability characteristics are improved significantly. Both Dirichlet boundary condition (for incident waves) and Neumann boundary condition (for a reflecting wall) are implemented. Several examples are presented to demonstrate the range of applicabilities and the accuracy of the model. Copyright © 2001 John Wiley & Sons, Ltd.     

A high‐order Petrov–Galerkin finite element method for the classical Boussinesq wave model

A high‐order Petrov–Galerkin finite element scheme is presented to solve the one‐dimensional depth‐integrated classical Boussinesq equations for weakly non‐linear and weakly dispersive waves. Finite elements are used both in the space and the time domains. The shape functions are bilinear in space–time, whereas the weighting functions are linear in space and quadratic in time, with C⁰‐continuity. Dispersion correction and a highly selective dissipation mechanism are introduced through additional streamline upwind terms in the weighting functions. An implicit, conditionally stable, one‐step predictor–corrector time integration scheme results. The accuracy and stability of the non‐linear discrete equations are investigated by means of a local Taylor series expansion. A linear spectral analysis is used for the full characterization of the predictor–corrector inner iterations. Based on the order of the analytical terms of the Boussinesq model and on the order of the numerical discretization, it is concluded that the scheme is fourth‐order accurate in terms of phase velocity. The dissipation term is third order only affecting the shortest wavelengths. A numerical convergence analysis showed a second‐order convergence rate in terms of both element size and time step. Four numerical experiments are addressed and their results are compared with analytical solutions or experimental data available in the literature: the propagation of a solitary wave, the oscillation of a flat bottom closed basin, the oscillation of a non‐flat bottom closed basin, and the propagation of a periodic wave over a submerged bar. Copyright © 2008 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.     

Positivity‐preserving,flux‐limited finite‐difference and finite‐element methods for reactive transport

A new class of positivity‐preserving, flux‐limited finite‐difference and Petrov–Galerkin (PG) finite‐element methods are devised for reactive transport problems.The methods are similar to classical TVD flux‐limited schemes with the main difference being that the flux‐limiter constraint is designed to preserve positivity for problems involving diffusion and reaction. In the finite‐element formulation, we also consider the effect of numerical quadrature in the lumped and consistent mass matrix forms on the positivity‐preserving property. Analysis of the latter scheme shows that positivity‐preserving solutions of the resulting difference equations can only be guaranteed if the flux‐limited scheme is both implicit and satisfies an additional lower‐bound condition on time‐step size. We show that this condition also applies to standard Galerkin linear finite‐element approximations to the linear diffusion equation. Numerical experiments are provided to demonstrate the behavior of the methods and confirm the theoretical conditions on time‐step size, mesh spacing, and flux limiting for transport problems with and without nonlinear reaction. Copyright © 2003 John Wiley & Sons, Ltd.     

Interface‐preserving level‐set reinitialization for DG‐FEM

This paper presents a contribution to level‐set reinitialization in the context of discontinuous Galerkin finite element methods. We focus on high‐order polynomials for the discretization and level set geometries, which are comparable to the element size. In contrast to hyperbolic and geometric reinitialization techniques, our method relies on solving a nonlinear elliptic PDE iteratively. We critically compare two different variants of the algorithm experimentally in numerical studies. The results demonstrate that the method is stable for nontrivial test cases and shows high‐order accuracy. Copyright © 2016 John Wiley & Sons, Ltd.     

Further experiences with computing non‐hydrostatic free‐surface flows involving water waves

A semi‐implicit, staggered finite volume technique for non‐hydrostatic, free‐surface flow governed by the incompressible Euler equations is presented that has a proper balance between accuracy, robustness and computing time. The procedure is intended to be used for predicting wave propagation in coastal areas. The splitting of the pressure into hydrostatic and non‐hydrostatic components is utilized. To ease the task of discretization and to enhance the accuracy of the scheme, a vertical boundary‐fitted co‐ordinate system is employed, permitting more resolution near the bottom as well as near the free surface. The issue of the implementation of boundary conditions is addressed. As recently proposed by the present authors, the Keller‐box scheme for accurate approximation of frequency wave dispersion requiring a limited vertical resolution is incorporated. The both locally and globally mass conserved solution is achieved with the aid of a projection method in the discrete sense. An efficient preconditioned Krylov subspace technique to solve the discretized Poisson equation for pressure correction with an unsymmetric matrix is treated. Some numerical experiments to show the accuracy, robustness and efficiency of the proposed method are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

Finite analytic boundary condition for a higher‐order finite analytic scheme

A higher‐order finite analytic scheme based on one‐dimensional finite analytic solutions is used to discretize three‐dimensional equations governing turbulent incompressible free surface flow. In order to preserve the accuracy of the numerical scheme, a new, finite analytic boundary condition is proposed for an accurate numerical solution of the partial differential equation. This condition has higher‐order accuracy. Thus, the same order of accuracy is used for the boundary. Boundary conditions were formulated and derived for fluid inflow, outflow, impermeable surfaces and symmetry planes. The derived boundary conditions are treated implicitly and updated with the solution of the problem. The basic idea for the derivation of boundary conditions was to use the discretized form of the governing equations for the fluid flow simplified on the boundaries and flow information. To illustrate the influence of the higher‐order effects at the boundaries, another, lower‐order finite analytic boundary condition, is suggested. The simulations are performed to demonstrate the validity of the present scheme and boundary conditions for a Wigley hull advancing in calm water. Copyright © 2005 John Wiley & Sons, Ltd.     

A modified Galerkin/finite element method for the numerical solution of the Serre‐Green‐Naghdi system

A new modified Galerkin/finite element method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low‐order Lagrange finite element spaces, despite the fact that the system contains third order spatial partial derivatives for the depth averaged velocity of the fluid. After studying the efficacy and the conservation properties of the new numerical method, we proceed with the validation of the new numerical model and boundary conditions by comparing the numerical solutions with laboratory experiments and with available theoretical asymptotic results. Copyright © 2016 John Wiley & Sons, Ltd.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

A FE‐based algorithm for the inverse natural convection problem

Several numerical algorithms for solving inverse natural convection problems are revisited and studied. Our aim is to identify the unknown strength of a time‐varying heat source via a set of coupled nonlinear partial differential equations obtained by the so‐called finite element consistent splitting scheme (CSS) in order to get a good approximation of the unknown heat source from both the measured data and model results, by minimizing a functional that measures discrepancies between model and measured data. Viewed as an optimization problem, the solutions are obtained by means of the conjugate gradient method. A second‐order CSS in time involving the direct problem, the adjoint problem, the sensitivity problem and a system of sensitivity functions is used in order to enhance the numerical accuracy obtained for the unknown heat source function. A spatial discretization of all field equations is implemented using equal‐order and mixed finite element methods. Numerical experiments validate the proposed optimization algorithms that are in good agreement with the existing results. Copyright © 2010 John Wiley & Sons, Ltd.     

A high‐order mass‐lumping procedure for B‐spline collocation method with application to incompressible flow simulations

This paper presents new developments of the staggered spline collocation method for cost‐effective solution to the incompressible Navier–Stokes equations. Maximal decoupling of the velocity and the pressure is obtained by using the fractional step method of Gresho and Chan, allowing the solution to sparse elliptic problems only. In order to preserve the high‐accuracy of the B‐spline method, this fractional step scheme is used in association with a sparse approximation to the inverse of the consistent mass matrix. Such an approximation is constructed from local spline interpolation method, and represents a high‐order generalization of the mass‐lumping technique of the finite‐element method. A numerical investigation of the accuracy and the computational efficiency of the resulting semi‐consistent spline collocation schemes is presented. These schemes generate a stable and accurate unsteady Navier–Stokes solver, as assessed by benchmark computations. Copyright © 2003 John Wiley & Sons, Ltd.     

Extension of domain‐free discretization method to simulate compressible flows over fixed and moving bodies

This paper is the first endeavour to present the local domain‐free discretization (DFD) method for the solution of compressible Navier–Stokes/Euler equations in conservative form. The discretization strategy of DFD is that for any complex geometry, there is no need to introduce coordinate transformation and the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior dependent points are updated at each time step to impose the wall boundary condition by the approximate form of solution near the boundary. Some points inside the solution domain are constructed for the approximate form of solution, and the flow variables at constructed points are evaluated by the linear interpolation on triangles. The numerical schemes used in DFD are the finite element Galerkin method for spatial discretization and the dual‐time scheme for temporal discretization. Some numerical results of compressible flows over fixed and moving bodies are presented to validate the local DFD method. Copyright © 2006 John Wiley & Sons, Ltd.     

A local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.     

A local domain‐free discretization method to simulate three‐dimensional compressible inviscid flows

In this paper, the domain‐free discretization method (DFD) is extended to simulate the three‐dimensional compressible inviscid flows governed by Euler equations. The discretization strategy of DFD is that the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior‐dependent points are updated at each time step by extrapolation along the wall normal direction in conjunction with the wall boundary conditions and the simplified momentum equation in the vicinity of the wall. Spatial discretization is achieved with the help of the finite element Galerkin approximation. The concept of ‘osculating plane’ is adopted, with which the local DFD can be easily implemented for the three‐dimensional case. Geometry‐adaptive tetrahedral mesh is employed for three‐dimensional calculations. Finally, we validate the DFD method for three‐dimensional compressible inviscid flow simulations by computing transonic flows over the ONERA M6 wing. Comparison with the reference experimental data and numerical results on boundary‐conforming grid was displayed and the results show that the present DFD results compare very well with the reference data. Copyright © 2009 John Wiley & Sons, Ltd.     

An improved KGF‐SPH with a novel discrete scheme of Laplacian operator for viscous incompressible fluid flows

The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

FEM‐BEM coupling for the exterior Stokes problem with non‐conforming finite elements and an application to small droplet deformation dynamics

A non‐conforming, discontinuous Galerkin finite element–boundary element coupling procedure is presented for the exterior planar Stokes problem. The novel coupled formulation is developed using that for the conforming case as a guide to the introduction of extra mortar variables used to couple a discontinuous interior finite element solution with a continuous exterior boundary element solution. Convergence results for the new scheme are presented, for a range of different interior penalties, on computational domains discretized with regular structured meshes. To illustrate an application, the excitations required to model two‐phase droplet deformations in an extensional flow, under simple surface tension, with the new scheme are also presented. For a selection of different drop viscocities and exterior flows, with and without a rotational component, the progression to a steady‐state deformation of initially undeformed circular drops is calculated and the results compared with those from both a conforming FEM‐BEM equivalent scheme and from a small perturbation analysis where available. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Developing a unified FVE‐ALE approach to solve unsteady fluid flow with moving boundaries

In this study, an arbitrary Lagrangian–Eulerian (ALE) approach is incorporated with a mixed finite‐volume–element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non‐stationary meshes. The method collects the advantages of both finite‐volume and finite‐element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical‐influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving triangular and quadrilateral elements, either partially or entirely. The use of hybrid FE grids facilitates this achievement. To show the robustness of the unified algorithm, we examine both the first‐ and the second‐order temporal stencils. The accuracy and performance of the extended method are evaluated via simulating the unsteady flow fields around a fixed cylinder, a transversely oscillating cylinder, and in a channel with an indented wall. The numerical results presented demonstrate significant accuracy benefits for the new hybrid method on coarse meshes and where large time steps are taken. Of importance, the current method yields the second‐order temporal accuracy when the second‐order stencil is used to discretize the unsteady terms. Copyright © 2009 John Wiley & Sons, Ltd.     

A deconvolution‐based fourth‐order finite volume method for incompressible flows on non‐uniform grids

This paper is concerned with the development of a new high‐order finite volume method for the numerical simulation of highly convective unsteady incompressible flows on non‐uniform grids. Specifically, both a high‐order fluxes integration and the implicit deconvolution of the volume‐averaged field are considered. This way, the numerical solution effectively stands for a fourth‐order approximation of the point‐wise one. Moreover, the procedure is developed in the framework of a projection method for the pressure–velocity decoupling, while originally deriving proper high‐order intermediate boundary conditions. The entire numerical procedure is discussed in detail, giving particular attention to the consistent discretization of the deconvolution operation. The present method is also cast in the framework of approximate deconvolution modelling for large‐eddy simulation. The overall high accuracy of the method, both in time and space, is demonstrated. Finally, as a model of real flow computation, a two‐dimensional time‐evolving mixing layer is simulated, with and without sub‐grid scales modelling. Copyright © 2003 John Wiley & Sons, Ltd.