The MFEM (Modular Finite Element Methods) library is a high-performance C++ library for finite element discretizations. MFEM supports numerous types of finite element methods and is the discretization engine powering many computational physics and engineering applications across a number of domains. This paper describes some of the recent research and development in MFEM, focusing on performance portability across leadership-class supercomputing facilities, including exascale supercomputers, as well as new capabilities and functionality, enabling a wider range of applications. Much of this work was undertaken as part of the Department of Energy’s Exascale Computing Project (ECP) in collaboration with the Center for Efficient Exascale Discretizations (CEED).

We investigate the flow over smooth (non-ribletted) shark skin denticles in an open-channel flow using direct numerical simulation (DNS) and two Reynolds averaged Navier–Stokes (RANS) closures. Large peaks in pressure and viscous drag are observed at the denticle crown edges, where they are exposed to high-speed fluid which penetrates between individual denticles, increasing shear and turbulence. Strong lift forces lead to a positive spanwise torque acting on individual denticles, potentially encouraging bristling if the denticles were not fixed. However, DNS predicts that denticles ultimately increase drag by 58% compared to a flat plate. Good predictions of drag distributions are obtained by RANS models, although an underestimation of turbulent kinetic energy production leads to an underprediction of drag. Nevertheless, RANS methods correctly predict trends in the drag data and the regions contributing most to viscous and pressure drag. Subsequently, RANS models are used to investigate the dependence of drag on the flow blockage ratio (boundary layer to roughness height ratio), finding that the drag increase due to denticles is halved when the blockage ratio δ / h is increased from 14 to 45. Our results provide an integrated understanding of the drag over non-ribletted denticles, enabling existing diverse drag data to be explained.