An adaptive spectral/hp discontinuous Galerkin method for the two-dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non-conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p-1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h-type refinement, the parent element is subdivided into four similar sibling elements. The time-stepping is performed using a third-order Runge-Kutta scheme. The performance of the hp-adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p-adaptivity is more efficient than h-adaptivity with respect to degrees of freedom and computational time.
Computational fluid dynamics (CFD) is becoming an increasingly popular tool in the wave energy sector, and over the last five years we have seen many studies using CFD. While the focus of the CFD studies have been on the validation phase, comparing numerically obtained results against experimental tests, the uncertainties associated with the numerical solution has so far been more or less overlooked. There is a need to increase the reliability of the numerical solutions in order to perform simulation based optimization at early stages of development. In this paper we introduce a well-established verification and validation (V&V) technique. We focus on the solution verification stage and how to estimate spatial discretization errors for simulations where no exact solutions are available. The technique is applied to the cases of a 2D heaving box and a 3D moored cylinder. The uncertainties are typically acceptable with a few percent for the 2D box, while the 3D cylinder case show double digit uncertainties. The uncertainties are discussed with regard to physical features of the flow and numerical techniques.