Search

This work presents a hybrid pressure face-centred finite volume (FCFV) solver to simulate steady-state incompressible Navier-Stokes flows. The method leverages the robustness, in the incompressible limit, of the hybridisable discontinuous Galerkin paradigm for compressible and weakly compressible flows to derive the formulation of a novel, low-order face-based discretisation. The incompressibility constraint is enforced in a weak sense, by introducing an inter-cell mass flux defined in terms of a new, hybrid variable, representing the pressure at the cell faces. This results in a new hybridisation strategy where cell variables (velocity, pressure and deviatoric strain rate tensor) are expressed as a function of velocity and pressure at the barycentre of the cell faces. The hybrid pressure formulation provides first-order convergence of all variables, including the stress, independently of cell type, stretching and distortion. Numerical benchmarks of Navier-Stokes flows at low and moderate Reynolds numbers, in two and three dimensions, are presented to evaluate accuracy and robustness of the method. In particular, the hybrid pressure formulation outperforms the FCFV method when convective effects are relevant, achieving accurate predictions on significantly coarser meshes., Comment: 44 pages, 18 figures, 6 tables

Purpose: This study aims to assess the accuracy of degree adaptive strategies in the context of incompressible Navier-Stokes flows using the high order hybridisable discontinuous Galerkin (HDG) method. Design/methodology/approach: The work presents a series of numerical examples to show the inability of standard degree adaptive processes to accurate capture aerodynamic quantities of interest, in particular the drag. A new conservative projection is proposed and the results between a standard degree adaptive procedure and the adaptive process enhanced with this correction are compared. The examples involve two transient problems where flow vortices or a gust needs to be accurately propagated over long distances. \noindent \textbf{}Findings:polynomials with a lower degree. Due to the coupling of velocity-pressure in incompressible flows, the violation of the incompressibility constraint leads to inaccurate pressure fields in the wake that have a sizeable effect on the drag. The new conservative projection proposed is found to remove all the numerical artefacts shown by the standard adaptive process. Originality/value: This work proposes a new conservative projection for the degree adaptive process. The projection does not introduce a significant overhead because it requires to solve an element-by-element problem and only for those elements where the adaptive process lowers the degree of approximation. Numerical results show that with the proposed projection non-physical oscillations in the drag disappear and the results are in good agreement with reference solutions., Comment: 35 pages, 24 figures, 1 table

This work presents a methodology to predict a near-optimal spacing function, which defines the element sizes, suitable to perform steady RANS turbulent viscous flow simulations. The strategy aims at utilising existing high fidelity simulations to compute a target spacing function and train an artificial neural network (ANN) to predict the spacing function for new simulations, either unseen operating conditions or unseen geometric configurations. Several challenges induced by the use of highly stretched elements are addressed. The final goal is to substantially reduce the time and human expertise that is nowadays required to produce suitable meshes for simulations. Numerical examples involving turbulent compressible flows in two dimensions are used to demonstrate the ability of the trained ANN to predict a suitable spacing function. The influence of the NN architecture and the size of the training dataset are discussed. Finally, the suitability of the predicted meshes to perform simulations is investigated.

This work develops, for the first time, a face-centred finite volume (FCFV) solver for the simulation of laminar and turbulent viscous incompressible flows. The formulation relies on the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the negative Spalart-Allmaras (SA) model and three novel convective stabilisations, inspired by Riemann solvers, are derived and compared numerically. The resulting method achieves first-order convergence of the velocity, the velocity-gradient tensor and the pressure. FCFV accurately predicts engineering quantities of interest, such as drag and lift, on unstructured meshes and, by avoiding gradient reconstruction, the method is less sensitive to mesh quality than other FV methods, even in the presence of highly distorted and stretched cells. A monolithic and a staggered solution strategies for the RANS-SA system are derived and compared numerically. Numerical benchmarks, involving laminar and turbulent, steady and transient cases are used to assess the performance, accuracy and robustness of the proposed FCFV method., Comment: 42 pages, 25 figures, 4 tables

This work presents the face-centred finite volume (FCFV) paradigm for the simulation of compressible flows. The FCFV method defines the unknowns at the face barycentre and uses a hybridisation procedure to eliminate all the degrees of freedom inside the cells. In addition, Riemann solvers are defined implicitly within the expressions of the numerical fluxes. The resulting methodology provides first-order accurate approximations of the conservative quantities, i.e. density, momentum and energy, as well as of the viscous stress tensor and of the heat flux, without the need of any gradient reconstruction procedure. Hence, the FCFV solver preserves the accuracy of the approximation in presence of distorted and highly stretched cells, providing a solver insensitive to mesh quality. In addition, FCFV is capable of constructing non-oscillatory approximations of sharp discontinuities without resorting to shock capturing or limiting techniques. For flows at low Mach number, the method is robust and is capable of computing accurate solutions in the incompressible limit without the need of introducing specific pressure correction strategies. A set of 2D and 3D benchmarks of external flows is presented to validate the methodology in different flow regimes, from inviscid to viscous laminar flows, from transonic to subsonic incompressible flows, demonstrating its potential to handle compressible flows in realistic scenarios., Comment: 38 pages, 22 figures

This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems., Comment: 90 pages, 51 figures

This work presents a review of high-order hybridisable discontinuous Galerkin (HDG) methods in the context of compressible flows. Moreover, an original unified framework for the derivation of Riemann solvers in hybridised formulations is proposed. This framework includes, for the first time in an HDG context, the HLL and HLLEM Riemann solvers as well as the traditional Lax-Friedrichs and Roe solvers. HLL-type Riemann solvers demonstrate their superiority with respect to Roe in supersonic cases due to their positivity preserving properties. In addition, HLLEM specifically outstands in the approximation of boundary layers because of its shear preservation, which confers it an increased accuracy with respect to HLL and Lax-Friedrichs. A comprehensive set of relevant numerical benchmarks of viscous and inviscid compressible flows is presented. The test cases are used to evaluate the competitiveness of the resulting high-order HDG scheme with the aforementioned Riemann solvers and equipped with a shock treatment technique based on artificial viscosity., Comment: 60 pages, 31 figures. arXiv admin note: substantial text overlap with arXiv:1912.00044

Reduced order models (ROM) are commonly employed to solve parametric problems and to devise inexpensive response surfaces to evaluate quantities of interest in real-time. There are many families of ROMs in the literature and choosing among them is not always a trivial task. This work presents a comparison of the performance of a priori and a posteriori proper generalised decomposition (PGD) algorithms for an incompressible Stokes flow problem in a geometrically parametrised domain. This problem is particularly challenging as the geometric parameters affect both the solution manifold and the computational spatial domain. The difficulty is further increased because multiple geometric parameters are considered and extended ranges of values are analysed for the parameters and this leads to significant variations in the flow features. Using a set of numerical experiments involving geometrically parametrised microswimmers, the two PGD algorithms are extensively compared in terms of their accuracy and their computational cost, expressed as a function of the number of full-order solves required., Comment: 58 pages, 31 figures

This paper proposes a novel computational framework for the solution of geometrically parametrised flow problems governed by the Stokes equation. The proposed method uses a high-order hybridisable discontinuous Galerkin formulation and the proper generalised decomposition rationale to construct an off-line solution for a given set of geometric parameters. The generalised solution contains the information for all the geometric parameters in a user-defined range and it can be used to compute sensitivities. The proposed approach circumvents many of the weaknesses of other approaches based on the proper generalised decomposition for computing generalised solutions of geometrically parametrised problems. Four numerical examples show the optimal approximation properties of the proposed method and demonstrate its applicability in two and three dimensions., Comment: 51 pages, 30 figures

An a priori reduced order method based on the proper generalised decomposition (PGD) is proposed to compute parametric solutions involving turbulent incompressible flows of interest in an industrial context, using OpenFOAM. The PGD framework is applied for the first time to the incompressible Navier-Stokes equations in the turbulent regime, to compute a generalised solution for velocity, pressure and turbulent viscosity, explicitly depending on the design parameters of the problem. In order to simulate flows of industrial interest, a minimally intrusive implementation based on OpenFOAM SIMPLE algorithm applied to the Reynolds-averaged Navier-Stokes equations with the Spalart-Allmaras turbulence model is devised. The resulting PGD strategy is applied to parametric flow control problems and achieves both qualitative and quantitative agreement with the full order OpenFOAM solution for convection-dominated fully-developed turbulent incompressible flows, with Reynolds number up to one million., Comment: 44 pages, 13 figures, 2 tables

A second-order face-centred finite volume strategy on general meshes is proposed. The method uses a mixed formulation in which a constant approximation of the unknown is computed on the faces of the mesh. Such information is then used to solve a set of problems, independent cell-by-cell, to retrieve the local values of the solution and its gradient. The main novelty of this approach is the introduction of a new basis function, utilised for the linear approximation of the primal variable in each cell. Contrary to the commonly used nodal basis, the proposed basis is suitable for computations on general meshes, including meshes with different cell types. The resulting approach provides second-order accuracy for the solution and first-order for its gradient, without the need of reconstruction procedures, is robust in the incompressible limit and insensitive to cell distortion and stretching. The second-order accuracy of the solution is exploited to devise an automatic mesh adaptivity strategy. An efficient error indicator is obtained from the computation of one extra local problem, independent cell-by-cell, and is used to drive mesh adaptivity. Numerical examples illustrating the approximation properties of the method and of the mesh adaptivity procedure are presented. The potential of the proposed method with automatic mesh adaptation is demonstrated in the context of microfluidics., Comment: 32 pages, 23 figures, 6 tables

This work proposes a high-order hybridised discontinuous Galerkin (HDG) formulation of the Harten-Lax-Van Leer (HLL) Riemann solver for compressible flows. A unified framework is introduced to present Lax-Friedrichs, Roe and HLL Riemann solvers via appropriate definitions of the HDG numerical fluxes. The resulting high-order HDG method with HLL Riemann solver is evaluated through a set of numerical simulations of inviscid compressible flows in different regimes, from subsonic isentropic flows to transonic and supersonic problems with shocks. The accuracy of the proposed method is comparable with the one of Lax-Friedrichs and Roe numerical fluxes in subsonic and transonic flows. The superior performance of HLL is highlighted in supersonic cases, where the method provides extra robustness, being able to produce positivity preserving approximations without the need of any user-defined entropy fix., Comment: An extended and updated version of the paper covering both inviscid and viscous compressible flows is available at arXiv:2009.06396

A hybridizable discontinuous Galerkin (HDG) formulation of the linearized incompressible Navier-Stokes equations, known as Oseen equations, is presented. The Cauchy stress formulation is considered and the symmetry of the stress tensor and the mixed variable, namely the scaled strain-rate tensor, is enforced pointwise via Voigt notation. Using equal-order polynomial approximations of degree k for all variables, HDG provides a stable discretization. Moreover, owing to Voigt notation, optimal convergence of order k+1 is obtained for velocity, pressure and strain-rate tensor and a local postprocessing strategy is devised to construct an approximation of the velocity superconverging with order k+2, even for low-order polynomial approximations. A tutorial for the numerical solution of incompressible flow problems using HDG is presented, with special emphasis on the technical details required for its implementation., Comment: 35 pages, 13 figures

Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV)., Comment: 29 pages, 7 figures

A second-order face-centred finite volume method (FCFV) is proposed. Contrary to the more popular cell-centred and vertex-centred finite volume (FV) techniques, the proposed method defines the solution on the faces of the mesh (edges in two dimensions). The method is based on a mixed formulation and therefore considers the solution and its gradient as independent unknowns. They are computed solving an element-by-element problem after the solution at the faces is determined. The proposed approach avoids the need of reconstructing the solution gradient, as required by cell-centred and vertex-centred FV methods. This strategy leads to a method that is insensitive to mesh distortion and stretching. The current method is second-order and requires the solution of a global system of equations of identical size and identical number of non-zero elements when compared to the recently proposed first-order FCFV. The formulation is presented for Poisson and Stokes problems. Numerical examples are used to illustrate the approximation properties of the method as well as to demonstrate its potential in three dimensional problems with complex geometries. The integration of a mesh adaptive procedure in the FCFV solution algorithm is also presented., Comment: 31 pages, 24 figures

The computational cost of parametric studies currently represents the major limitation to the application of simulation-based engineering techniques in a daily industrial environment. This work presents the first nonintrusive implementation of the proper generalised decomposition (PGD) in OpenFOAM, for the approximation of parametrised laminar incompressible Navier-Stokes equations. The key feature of this approach is the seamless integration of a reduced order model (ROM) in the framework of an industrially validated computational fluid dynamics software. This is of special importance in an industrial environment because in the online phase of the PGD ROM the description of the flow for a specific set of parameters is obtained simply via interpolation of the generalised solution, without the need of any extra solution step. On the one hand, the spatial problems arising from the PGD separation of the unknowns are treated using the classical solution strategies of OpenFOAM, namely the semi-implicit method for pressure linked equations (SIMPLE) algorithm. On the other hand, the parametric iteration is solved via a collocation approach. The resulting ROM is applied to several benchmark tests of laminar incompressible Navier-Stokes flows, in two and three dimensions, with different parameters affecting the flow features. Eventually, the capability of the proposed strategy to treat industrial problems is verified by applying the methodology to a parametrised flow control in a realistic geometry of interest for the automotive industry., Comment: 35 pages, 17 figures

The main objective of this work is to describe a general and original approach for computing an off-line solution for a set of parameters describing the geometry of the domain. That is, a solution able to include information for different geometrical parameter values and also allowing to compute readily the sensitivities. Instead of problem dependent approaches, a general framework is presented for standard engineering environments where the geometry is defined by means of NURBS. The parameters controlling the geometry are now the control points characterising the NURBS curves or surfaces. The approach proposed here, valid for 2D and 3D scenarios, allows a seamless integration with CAD preprocessors. The proper generalised decomposition (PGD), which is applied here to compute explicit geometrically parametrised solutions, circumvents the curse of dimensionality. Moreover, optimal convergence rates are shown for PGD approximations of incompressible flows., Comment: 45 pages, 34 figures

A face-centred finite volume (FCFV) method is proposed for the linear elasticity equation. The FCFV is a mixed hybrid formulation, featuring a system of first-order equations, that defines the unknowns on the faces (edges in two dimensions) of the mesh elements. The symmetry of the stress tensor is strongly enforced using the well-known Voigt notation and the displacement and stress fields inside each cell are obtained element-wise by means of explicit formulas. The resulting FCFV method is robust and locking-free in the nearly incompressible limit. Numerical experiments in two and three dimensions show optimal convergence of the displacement and the stress fields without any reconstruction. Moreover, the accuracy of the FCFV method is not sensitive to mesh distortion and stretching. Classical benchmark tests including Kirch's plate and Cook's membrane problems in two dimensions as well as three dimensional problems involving shear phenomenons, pressurised thin shells and complex geometries are presented to show the capability and potential of the proposed methodology., Comment: 29 pages, 20 figures

This work proposes a superconvergent hybridizable discontinuous Galerkin (HDG) method for the approximation of the Cauchy formulation of the Stokes equation using same degree of polynomials for the primal and mixed variables. The novel formulation relies on the well-known Voigt notation to strongly enforce the symmetry of the stress tensor. The proposed strategy introduces several advantages with respect to the existing HDG formulations. First, it remedies the suboptimal behavior experienced by the classical HDG method for formulations involving the symmetric part of the gradient of the primal variable. The optimal convergence of the mixed variable is retrieved and an element-by-element post-process procedure leads to a superconvergent velocity field, even for low-order approximations. Second, no additional enrichment of the discrete spaces is required and a gain in computational efficiency follows from reducing the quantity of stored information and the size of the local problems. Eventually, the novel formulation naturally imposes physical tractions on the Neumann boundary. Numerical validation of the optimality of the method and its superconvergent properties is performed in 2D and 3D using meshes of different element types., Comment: 28 pages, 10 figures, 1 table

The first super-convergent hybridisable discontinuous Galerkin (HDG) method for linear elastic problems capable of using the same degree of approximation for both the primal and mixed variables is presented. The key feature of the method is the strong imposition of the symmetry of the stress tensor by means of the well-known and extensively used Voigt notation, circumventing the use of complex mathematical concepts to enforce the symmetry of the stress tensor either weakly or strongly. A novel procedure to construct element-by-element a super-convergent post-processed displacement is proposed. Contrary to other HDG formulations, the methodology proposed here is able to produce a super-convergent displacement field for low order approximations. The resulting method is robust and locking-free in the nearly-incompressible limit. An extensive set of numerical examples is utilised to provide evidence of the optimality of the method and its super-convergent properties in two and three dimensions and for different element types., Comment: 37 pages, 18 figures

The NURBS-enhanced finite element method (NEFEM) combined with a hybridisable discontinuous Galerkin (HDG) approach is presented for the first time. The proposed technique completely eliminates the uncertainty induced by a polynomial approximation of curved boundaries that is common within an isoparametric approach and, compared to other DG methods, provides a significant reduction in number of degrees of freedom. In addition, by exploiting the ability of HDG to compute a postprocessed solution and by using a local a priori error estimate valid for elliptic problems, an inexpensive, reliable and computable error estimator is devised. The proposed methodology is used to solve Stokes flow problems using automatic degree adaptation. Particular attention is paid to the importance of an accurate boundary representation when changing the degree of approximation in curved elements. Several strategies are compared and the superiority and reliability of HDG-NEFEM with degree adaptation is illustrated., Comment: 36 pages, 37 figures

This work proposes a novel finite volume paradigm, the face-centred finite volume (FCFV) method. Contrary to the popular vertex (VCFV) and cell (CCFV) centred finite volume methods, the novel FCFV defines the solution on the mesh faces (edges in 2D) to construct locally-conservative numerical schemes. The idea of the FCFV method stems from a hybridisable discontinuous Galerkin (HDG) formulation with constant degree of approximation, thus inheriting the convergence properties of the classical HDG. The resulting FCFV features a global problem in terms of a piecewise constant function defined on the faces of the mesh. The solution and its gradient in each element are then recovered by solving a set of independent element-by-element problems. The mathematical formulation of FCFV for Poisson and Stokes equation is derived and numerical evidence of optimal convergence in 2D and 3D is provided. Numerical examples are presented to illustrate the accuracy, efficiency and robustness of the proposed methodology. The results show that, contrary to other FV methods, the accuracy of the FCFV method is not sensitive to mesh distortion and stretching. In addition, the FCFV method shows its better performance, accuracy and robustness using simplicial elements, facilitating its application to problems involving complex geometries in 3D., Comment: 43 pages, 33 figures, 3 tables