861 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
eXtended hybridizable discontinuous Galerkin for incompressible flow problems with unfitted meshes and interfaces
The eXtended hybridizable discontinuous Galerkin (X-HDG) method is developed for the solution of Stokes problems with void or material interfaces. X-HDG is a novel method that combines the hybridizable discontinuous Galerkin (HDG) method with an eXtended finite element strategy, resulting in a high-order, unfitted, superconvergent method, with an explicit definition of the interface geometry by means of a level-set function. For elements not cut by the interface, the standard HDG formulation is applied, whereas a modified weak form for the local problem is proposed for cut elements. Heaviside enrichment is considered on cut faces and in cut elements in the case of bimaterial problems. Two-dimensional numerical examples demonstrate that the applicability, accuracy, and superconvergence properties of HDG are inherited in X-HDG, with the freedom of computational meshes that do not fit the interfacesPeer ReviewedPostprint (author's final draft
Superconvergent interpolatory HDG methods for reaction diffusion equations I: An HDG method
In our earlier work [8], we approximated solutions of a general class of
scalar parabolic semilinear PDEs by an interpolatory hybridizable discontinuous
Galerkin (Interpolatory HDG) method. This method reduces the computational cost
compared to standard HDG since the HDG matrices are assembled once before the
time integration. Interpolatory HDG also achieves optimal convergence rates;
however, we did not observe superconvergence after an element-by-element
postprocessing. In this work, we revisit the Interpolatory HDG method for
reaction diffusion problems, and use the postprocessed approximate solution to
evaluate the nonlinear term. We prove this simple change restores the
superconvergence and keeps the computational advantages of the Interpolatory
HDG method. We present numerical results to illustrate the convergence theory
and the performance of the method
Non-modal analysis of spectral element methods: Towards accurate and robust large-eddy simulations
We introduce a \textit{non-modal} analysis technique that characterizes the
diffusion properties of spectral element methods for linear
convection-diffusion systems. While strictly speaking only valid for linear
problems, the analysis is devised so that it can give critical insights on two
questions: (i) Why do spectral element methods suffer from stability issues in
under-resolved computations of nonlinear problems? And, (ii) why do they
successfully predict under-resolved turbulent flows even without a
subgrid-scale model? The answer to these two questions can in turn provide
crucial guidelines to construct more robust and accurate schemes for complex
under-resolved flows, commonly found in industrial applications. For
illustration purposes, this analysis technique is applied to the hybridized
discontinuous Galerkin methods as representatives of spectral element methods.
The effect of the polynomial order, the upwinding parameter and the P\'eclet
number on the so-called \textit{short-term diffusion} of the scheme are
investigated. From a purely non-modal analysis point of view, polynomial orders
between and with standard upwinding are well suited for under-resolved
turbulence simulations. For lower polynomial orders, diffusion is introduced in
scales that are much larger than the grid resolution. For higher polynomial
orders, as well as for strong under/over-upwinding, robustness issues can be
expected. The non-modal analysis results are then tested against under-resolved
turbulence simulations of the Burgers, Euler and Navier-Stokes equations. While
devised in the linear setting, our non-modal analysis succeeds to predict the
behavior of the scheme in the nonlinear problems considered
A hybridizable discontinuous Galerkin method for electromagnetics with a view on subsurface applications
Two Hybridizable Discontinuous Galerkin (HDG) schemes for the solution of
Maxwell's equations in the time domain are presented. The first method is based
on an electromagnetic diffusion equation, while the second is based on
Faraday's and Maxwell--Amp\`ere's laws. Both formulations include the diffusive
term depending on the conductivity of the medium. The three-dimensional
formulation of the electromagnetic diffusion equation in the framework of HDG
methods, the introduction of the conduction current term and the choice of the
electric field as hybrid variable in a mixed formulation are the key points of
the current study. Numerical results are provided for validation purposes and
convergence studies of spatial and temporal discretizations are carried out.
The test cases include both simulation in dielectric and conductive media
- …