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:

No abstract

We present an arbitrary-order spectral element method for general-purpose simulation of non-overturning water waves, described by fully nonlinear potential theory. The method can be viewed as a high-order extension of the classical finite element method proposed and irregular dispersive wave propagation. The benefit of using a high-order -possibly adapted -spatial discretisation for accurate water wave propagation over long times and distances is particularly attractive for marine hydrodynamics applications.1

SUMMARYWe present a spectral/hp element discontinuous Galerkin model for simulating shallow water flows on unstructured triangular meshes. The model uses an orthogonal modal expansion basis of arbitrary order for the spatial discretisation and a third-order Runge-Kutta scheme to advance in time. The local elements are coupled together by numerical fluxes, evaluated using the HLLC Riemann solver. We apply the model to test cases involving smooth flows and demonstrate the exponentially fast convergence with regard to polynomial order. We also illustrate that even for results of "engineering

Results from Blind Test Series 1, part of the Collaborative Computational Project in Wave Structure Interaction (CCP-WSI), are presented. Participants, with a range of numerical methods, simulate blindly the interaction between a fixed structure and focused waves ranging in steepness and direction. Numerical results are compared against corresponding physical data. The predictive capability of each method is assessed based on pressure and run-up measurements. In general, all methods perform well in the cases considered, however, there is notable variation in the results (even between similar methods). Recommendations are made for appropriate considerations and analysis in future comparative studies.

We present spectral/hp discontinuous Galerkin methods for modelling weakly nonlinear and dispersive water waves, described by a set of depth-integrated Boussinesq equations, on unstructured triangular meshes. When solving the equations two different formulations are considered: directly solving the coupled momentum equations and the 'scalar method', in which a wave continuity equation is solved as an intermediate step. We demonstrate that the approaches are fully equivalent and give identical results in terms of accuracy, convergence and restriction on the time step. However, the scalar method is shown to be more CPU efficient for high order expansions, in addition to requiring less storage.

Abstract: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 a 0 -C continuous expansion. 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 of 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.

This paper analyses the nonlinear forces on a moored point-absorbing wave energy converter (WEC) in resonance at prototype scale (1:1) and at model scale (1:16). Three simulation types were used: Reynolds Averaged Navier-Stokes (RANS), Euler and the linear radiation-diffraction method (linear). Results show that when the wave steepness is doubled, the response reduction is: (i) 3% due to the nonlinear mooring response and the Froude-Krylov force; (ii) 1-4% due to viscous forces; and (iii) 18-19% due to induced drag and non-linear added mass and radiation forces. The effect of the induced drag is shown to be largely scale-independent. It is caused by local pressure variations due to vortex generation below the body, which reduce the total pressure force on the hull. Euler simulations are shown to be scale-independent and the scale effects of the WEC are limited by the purely viscous contribution (1-4%) for the two waves studied. We recommend that experimental model scale test campaigns of WECs should be accompanied by RANS simulations, and the analysis complemented by scale-independent Euler simulations to quantify the scale-dependent part of the nonlinear effects.