Oceanic tides have the potential to yield a vast amount of renewable energy. Tidal stream generators are one of the key technologies for extracting and harnessing this potential. In order to extract an economically useful amount of power, hundreds of tidal turbines must typically be deployed in an array. This naturally leads to the question of how these turbines should be configured to extract the maximum possible power: the positioning and the individual tuning of the turbines could significantly influence the extracted power, and hence is of major economic interest. However, manual optimisation is difficult due to legal site constraints, nonlinear interactions of the turbine wakes, and the cubic dependence of the power on the flow speed. The novel contribution of this paper is the formulation of this problem as an optimisation problem constrained by a physical model, which is then solved using an efficient gradient-based optimisation algorithm. In each optimisation iteration, a two-dimensional finite element shallow water model predicts the flow and the performance of the current array configuration. The gradient of the power extracted with respect to the turbine positions and their tuning parameters is then computed in a fraction of the time taken for a flow solution by solving the associated adjoint equations. These equations propagate causality backwards through the computation, from the power extracted back to the turbine positions and the tuning parameters. This yields the gradient at a cost almost independent of the number of turbines, which is crucial for any practical application. The utility of the approach is demonstrated by optimising turbine arrays in four idealised scenarios and a more realistic case with up to 256 turbines in the Inner Sound of the Pentland Firth, Scotland.

Nonlinear systems of partial differential equations (PDEs) may permit several distinct solutions. The typical current approach to finding distinct solutions is to start Newton's method with many different initial guesses, hoping to find starting points that lie in different basins of attraction. In this paper, we present an infinite-dimensional deflation algorithm for systematically modifying the residual of a nonlinear PDE problem to eliminate known solutions from consideration. This enables the Newton-Kantorovitch iteration to converge to several different solutions, even starting from the same initial guess. The deflated Jacobian is dense, but an efficient preconditioning strategy is devised, and the number of Krylov iterations is observed not to grow as solutions are deflated. The power of the approach is demonstrated on several problems from special functions, phase separation, differential geometry, and fluid mechanics that permit distinct solutions.

The recently proposed glymphatic system suggests that bulk flow is important for clearing waste from the brain, and as such may underlie the development of e.g. Alzheimer's disease. The glymphatic hypothesis is still controversial and several biomechanical modeling studies at the micro-level have questioned the system and its assumptions. In contrast, at the macro-level, there are many experimental findings in support of bulk flow. Here, we will investigate to what extent the CSF tracer distributions seen in novel magnetic resonance imaging (MRI) investigations over hours and days are suggestive of bulk flow as an additional component to diffusion. In order to include the complex geometry of the brain, the heterogeneous CSF flow around the brain, and the transport over the timescale of days, we employed the methods of partial differential constrained optimization to identify the apparent diffusion coefficient (ADC) that would correspond best to the MRI findings. We found that the computed ADC in the cortical grey matter was 5-26% larger than the ADC estimated with DTI, which suggests that diffusion may not be the only mechanism governing transport.

No abstract

No abstract