a b s t r a c tNektar++ is an open-source software framework designed to support the development of highperformance scalable solvers for partial differential equations using the spectral/hp element method. High-order methods are gaining prominence in several engineering and biomedical applications due to their improved accuracy over low-order techniques at reduced computational cost for a given number of degrees of freedom. However, their proliferation is often limited by their complexity, which makes these methods challenging to implement and use. Nektar++ is an initiative to overcome this limitation by encapsulating the mathematical complexities of the underlying method within an efficient C++ framework, making the techniques more accessible to the broader scientific and industrial communities. The software supports a variety of discretisation techniques and implementation strategies, supporting methods research as well as application-focused computation, and the multi-layered structure of the framework allows the user to embrace as much or as little of the complexity as they need. The libraries capture the mathematical constructs of spectral/hp element methods, while the associated collection of pre-written PDE solvers provides out-of-the-box application-level functionality and a template for users who wish to develop solutions for addressing questions in their own scientific domains. Program summaryProgram title: Nektar++ Catalogue identifier: AEVV_v1_0Program summary URL:

We present an alternative "encapsulated" formulation of the Selective Frequency Damping method for finding unstable equilibria of dynamical systems, which is particularly useful when analysing the stability of fluid flows. The formulation makes use of splitting methods, which means that it can be wrapped around an existing time-stepping code as a "black box". The method is first applied to a scalar problem in order to analyse its stability and highlight the roles of the control coefficient χ and the filter width ∆ in the convergence (or not) towards the steady-state. Then the steady-state of the incompressible flow past a two-dimensional cylinder at Re = 100, obtained with a code which implements the spectral/hp element method, is presented.

The selective frequency damping (SFD) method is an alternative to classical Newton's method to obtain unstable steady-state solutions of dynamical systems. However this method has two main limitations: it does not converge for arbitrary control parameters; and when it does converge, the time necessary to reach the steady-state solution may be very long. In this paper we present an adaptive algorithm to address these two issues. We show that by evaluating the dominant eigenvalue of a "partially converged" steady flow, we can select a control coefficient and a filter width that ensure an optimum convergence of the SFD method. We apply this adaptive method to several classical test cases of computational fluid dynamics and we show that a steady-state solution can be obtained without any a priori knowledge of the flow stability properties

Accurate visualization of high-order meshes and flow fields is a fundamental tool for the verification, validation, analysis and interpretation of high-order flow simulations. Standard visualization tools based on piecewise linear approximations can be used for the display of high-order fields but their accuracy is restricted by computer memory and processing time. More often than not, the accurate visualization of complex flows using this strategy requires computational resources beyond the reach of most users. This chapter describes ElVis, a truly high-order and interactive visualization system created for the accurate and interactive visualization of scalar fields produced by high-order spectral/hp finite element simulations. We show some examples that motivate the need for such a visualization system and illustrate some of its features for the display and analysis of simulation data.

The aim of this work is to optimize an actuator design so that the flow profile at its exit section is as close as possible to a target profile. The method is founded on the penalization and level-set methods to solve direct and inverse problems on Cartesian meshes The optimization process is written and applied both for Stokes and Navier-Stokes flows. The results show that the method can be successfully applied to the non linear problem to improve the flow profile of an actuator even if the target cannot be totally reached.