In the past 15 years the field of general purpose computing on graphics processing units, or GPUs, has become well developed and the practice is becoming more mainstream. Computational demands of simulation software are continuously increasing. As such for many applications traditionally computed on the central processing unit the question arises of whether moving to GPU computing is a possible cost effective way of meeting these demands.
The fundamental nature of GPU architecture that makes it so cost effective at doing bulk computation also poses restrictions on which applications are suitable for it. The shallow water equations are a simplified form of the Navier-Stokes equations and describe water levels and flow currents in suitably shallow water such as rivers, estuaries and the North sea. The main research goal of this thesis project was to determine whether the shallow water equations are suitable for implementation on a GPU. Two options exist, the equations may be solved with either an explicit or implicit time integration method. First, a literature study was conducted to familiarize with the tools required to build explicit and implicit shallow water models on a GPU. Then both an explicit and implicit shallow water solver were developed first in the MATLAB programming language and later in CUDA C++ on both CPU and GPU. The main findings are that both explicit and implicit methods are well suited for GPU implementation. Both methods proved to be compatible with a wetting and drying mechanism of numerical cells. The Cuda C++ implementation was in the order of 10 times as fast as a MATLAB implementation for both CPU and GPU. For the benchmark cases tested, the Cuda C++ GPU implementation was in the order of 50 times faster than the equivalent multithreaded CPU implementation.
The implicit implementation was benchmarked using the conjugate gradient method to solve the linear system. Various preconditioners were tested and a Repeated Red Black preconditioner was found to be the most effective. The computation time of the RRB preconditioned implicit method was compared with the explicit method and it was found that the two methods reached parity in computation time when the implicit time step was taken roughly 50 times as large as the explicit time step. For implicit time steps smaller than that the explicit method was faster and when the implicit time step was larger the implicit method was faster. For the benchmark cases tested, the implicit method using a time step 50 times larger than the explicit method was found to be less accurate and less stable than the explicit method. The conclusion is that for cases similar to the benchmark cases an explicit method is the fastest, most stable and most accurate method and thus the preferred choice.

The aim of this study is to evaluate the performance of two methods for determining the speed of a ship relative to the water during a so called speed trial.
The first method is called the Mean of means method and has been traditionally. The second is called the Iterative method and has been recently developed and is already used for speed trial analysis.
The methods have been evaluated by simulation.
A relationship between ship power and velocity was chosen based on a theoretical analysis, and this relationship is used to generate simulated speed trial data. Using the knowledge of the ’true’ relationship the performance of the methods could then be evaluated after implementation.
The findings are that the Iterative method needs a minimum of 8 measurements obtained by sailing back and forth 4 times in order to find the true relationship with acceptable
accuracy.
The Mean of means method needs a minimum of 12 measurements,
or 6 double runs.
The conclusion is that the Iterative method is more time efficient
for speed trial analysis.
In a test case assuming realistic measurement noise, 99:3% of found speeds had an error smaller than the simulated error margin of 0.1 knots, out of a sample size of 1000.
When comparing the Iterative to the Mean of means method under the same conditions, the Mean of means method made smaller errors for all observed cases.
From the results of this study has been concluded that the Iterative method is usable method for determining ship speed during a sea trial that has both advantages and disadvantages over the traditional Mean of means method.
Advantages are that the Iterative method requires less measurements than the Mean of means method, and it remains accurate even if measurements have a large or varying spacing in time.
A disadvantage is that it is more sensitive to measurement noise than the Mean of means method, making it prone to making larger errors.
A small adjustment to the current function the method assumes and uses for calculating current effects has been proposed.
The adjustment was to normalize the component of the current function that scales linearly with time.
This has been found to have a positive impact on the performance of the method.