Search

Given an arbitrary function in |$\boldsymbol{H}({\operatorname{div}})$|⁠ , we show that the error attained by the global-best approximation by |$\boldsymbol{H}({\operatorname{div}})$| -conforming piecewise polynomial Raviart–Thomas–Nédélec elements under additional constraints on the divergence and normal flux on the boundary is, up to a generic constant, equivalent to the sum of independent local-best approximation errors over individual mesh elements, without constraints on the divergence or normal fluxes. The generic constant only depends on the shape-regularity of the underlying simplicial mesh, the space dimension, and the polynomial degree of the approximations. The analysis also gives rise to a stable, local, commuting projector in |$\boldsymbol{H}({\operatorname{div}})$|⁠ , delivering an approximation error that is equivalent to the local-best approximation. We next present a variant of the equivalence result, where robustness of the constant with respect to the polynomial degree is attained for unbalanced approximations. These two results together further enable us to derive rates of convergence of global-best approximations that are fully optimal in both the mesh size |$h$| and the polynomial degree |$p$|⁠ , for vector fields that only feature elementwise the minimal necessary Sobolev regularity. We finally show how to apply our findings to derive optimal a priori |$hp$| -error estimates for mixed and least-squares finite element methods applied to a model diffusion problem. [ABSTRACT FROM AUTHOR]

A variational discrete element method is applied to simulate quasi-static crack propagation. Cracks are considered to propagate between the mesh cells through the mesh facets. The elastic behaviour is parametrized by the continuous mechanical parameters (Young modulus and Poisson ratio). A discrete energetic cracking criterion coupled to a discrete kinking criterion guide the cracking process. Two-dimensional numerical examples are presented to illustrate the robustness and versatility of the method. [ABSTRACT FROM AUTHOR]

We investigate the computational performance of hybrid high-order methods applied to flow simulations in extremely large discrete fracture networks (over one million of fractures). We study the choice of basis functions, the trade-off between increasing the polynomial order and refining the mesh, and how to take advantage of polygonal cells to reduce the number of degrees of freedom. [ABSTRACT FROM AUTHOR]

We investigate artificial compressibility (AC) techniques for the time discretization of the incompressible Navier–Stokes equations. The space discretization is based on a lowest-order face-based scheme supporting polytopal meshes, namely discrete velocities are attached to the mesh faces and cells, whereas discrete pressures are attached to the mesh cells. This face-based scheme can be embedded into the framework of hybrid mixed mimetic schemes and gradient schemes, and has close links to the lowest-order version of hybrid high-order methods devised for the steady incompressible Navier–Stokes equations. The AC time-stepping uncouples at each time step the velocity update from the pressure update. The performances of this approach are compared against those of the more traditional monolithic approach which maintains the velocity-pressure coupling at each time step. We consider both first-order and second-order time schemes and either an implicit or an explicit treatment of the nonlinear convection term. We investigate numerically the CFL stability restriction resulting from an explicit treatment, both on Cartesian and polytopal meshes. Finally, numerical tests on large 3D polytopal meshes highlight the efficiency of the AC approach and the benefits of using second-order schemes whenever accurate discrete solutions are to be attained. [ABSTRACT FROM AUTHOR]

We design and analyze a hybrid high-order method on unfitted meshes to approximate the Stokes interface problem. The interface can cut through the mesh cells in a very general fashion. A cell-agglomeration procedure prevents the appearance of small cut cells. Our main results are inf-sup stability and a priori error estimates with optimal convergence rates in the energy norm. Numerical simulations corroborate these results. [ABSTRACT FROM AUTHOR]