25 research outputs found
Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity
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).Peer ReviewedPostprint (author's final draft
Robust stabilised finite element solvers for generalised Newtonian fluid flows
Various materials and solid-fluid mixtures of engineering and biomedical interest can be modelled as generalised Newtonian fluids, as their apparent viscosity depends locally on the flow field. Despite the particular features of such models, it is common practice to combine them with numerical techniques originally conceived for Newtonian fluids, which can bring several issues such as spurious pressure boundary layers, unsuitable natural boundary conditions and coupling terms spoiling the efficiency of nonlinear solvers and preconditioners. In this work, we present a finite element framework dealing with such issues while maintaining low computational cost and simple implementation. The building blocks of our algorithm are (i) an equal-order stabilisation method preserving consistency even for lowest-order discretisations, (ii) robust extrapolation of velocities in the time-dependent case to decouple the rheological law from the overall system, (iii) adaptive time step selection and (iv) a fast physics-based preconditioned Krylov subspace solver, to tackle the relevant range of discretisation parameters including highly varying viscosity. Selected numerical experiments are provided demonstrating the potential of our approach in terms of robustness, accuracy and efficiency for problems of practical interest
The face-centered finite volume method (FCFV) for steady-state incompressible Navier-Stokes equations
The face-centered finite volumes (FCFV) has been proposed [2,4]. FCFV may be derived as a hybridizable discontinuous Galerkin (HDG) method with constant degree of approximation for all the variables [1,3,5]. Contrary to CCFV and VCFV approaches, the proposed FCFV method provides LBB-stable discretizations and achieves first-order convergence of velocity, pressure and gradient of velocity using unstructured meshes, without the need to perform flux reconstruction. Numerical experiments show that the method is robust in presence of distorted and stretched elements [2,4]
HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB
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
Interpolatory HDG Method for Parabolic Semilinear PDEs
We propose the interpolatory hybridizable discontinuous Galerkin
(Interpolatory HDG) method for a class of scalar parabolic semilinear PDEs. The
Interpolatory HDG method uses an interpolation procedure to efficiently and
accurately approximate the nonlinear term. This procedure avoids the numerical
quadrature typically required for the assembly of the global matrix at each
iteration in each time step, which is a computationally costly component of the
standard HDG method for nonlinear PDEs. Furthermore, the Interpolatory HDG
interpolation procedure yields simple explicit expressions for the nonlinear
term and Jacobian matrix, which leads to a simple unified implementation for a
variety of nonlinear PDEs. For a globally-Lipschitz nonlinearity, we prove that
the Interpolatory HDG method does not result in a reduction of the order of
convergence. We display 2D and 3D numerical experiments to demonstrate the
performance of the method
Spectral/hp element methods: recent developments, applications, and perspectives
The spectral/hp element method combines the geometric flexibility of the
classical h-type finite element technique with the desirable numerical
properties of spectral methods, employing high-degree piecewise polynomial
basis functions on coarse finite element-type meshes. The spatial approximation
is based upon orthogonal polynomials, such as Legendre or Chebychev
polynomials, modified to accommodate C0-continuous expansions. Computationally
and theoretically, by increasing the polynomial order p, high-precision
solutions and fast convergence can be obtained and, in particular, under
certain regularity assumptions an exponential reduction in approximation error
between numerical and exact solutions can be achieved. This method has now been
applied in many simulation studies of both fundamental and practical
engineering flows. This paper briefly describes the formulation of the
spectral/hp element method and provides an overview of its application to
computational fluid dynamics. In particular, it focuses on the use the
spectral/hp element method in transitional flows and ocean engineering.
Finally, some of the major challenges to be overcome in order to use the
spectral/hp element method in more complex science and engineering applications
are discussed
Recommended from our members
Advanced Computational Engineering
The finite element method is the established simulation tool for the numerical solution of partial differential equations in many engineering problems with many mathematical developments such as mixed finite element methods (FEMs) and other nonstandard FEMs like least-squares, nonconforming, and discontinuous Galerkin (dG) FEMs. Various aspects on this plus related topics ranging from order-reduction methods to isogeometric analysis has been discussed amongst the pariticpants form mathematics and engineering for a large range of applications
Monolithic multiphysics simulation of hypersonic aerothermoelasticity using a hybridized discontinuous Galerkin method
This work presents implementation of a hybridized discontinuous Galerkin (DG) method for robust simulation of the hypersonic aerothermoelastic multiphysics system. Simulation of hypersonic vehicles requires accurate resolution of complex multiphysics interactions including the effects of high-speed turbulent flow, extreme heating, and vehicle deformation due to considerable pressure loads and thermal stresses. However, the state-of-the-art procedures for hypersonic aerothermoelasticity are comprised of low-fidelity approaches and partitioned coupling schemes. These approaches preclude robust design and analysis of hypersonic vehicles for a number of reasons. First, low-fidelity approaches limit their application to simple geometries and lack the ability to capture small scale flow features (e.g. turbulence, shocks, and boundary layers) which greatly degrades modeling robustness and solution accuracy. Second, partitioned coupling approaches can introduce considerable temporal and spatial inaccuracies which are not trivially remedied. In light of these barriers, we propose development of a monolithically-coupled hybridized DG approach to enable robust design and analysis of hypersonic vehicles with arbitrary geometries. Monolithic coupling methods implement a coupled multiphysics system as a single, or monolithic, equation system to be resolved by a single simulation approach. Further, monolithic approaches are free from the physical inaccuracies and instabilities imposed by partitioned approaches and enable time-accurate evolution of the coupled physics system. In this work, a DG method is considered due to its ability to accurately resolve second-order partial differential equations (PDEs) of all classes. We note that the hypersonic aerothermoelastic system is composed of PDEs of all three classes. Hybridized DG methods are specifically considered due to their exceptional computational efficiency compared to traditional DG methods. It is expected that our monolithic hybridized DG implementation of the hypersonic aerothermoelastic system will 1) provide the physical accuracy necessary to capture complex physical features, 2) be free from any spatial and temporal inaccuracies or instabilities inherent to partitioned coupling procedures, 3) represent a transition to high-fidelity simulation methods for hypersonic aerothermoelasticity, and 4) enable efficient analysis of hypersonic aerothermoelastic effects on arbitrary geometries