119 results found
Abolghasemi M, Piggott MD, Spinneken J, et al., Simulating tidal turbines with mesh optimisation and RANS turbulence models, 2015 European Wave and Tidal Energy Conference
Avdis A, Jacobs CT, Hill J, et al., Shoreline and Bathymetry Approximation in Mesh Generation for Tidal Renewable Simulations
Due to the fractal nature of the domain geometry in geophysical flowsimulations, a completely accurate description of the domain in terms of acomputational mesh is frequently deemed infeasible. Shoreline and bathymetrysimplification methods are used to remove small scale details in the geometry,particularly in areas away from the region of interest. To that end, a novelmethod for shoreline and bathymetry simplification is presented. Existingshoreline simplification methods typically remove points if the resultantgeometry satisfies particular geometric criteria. Bathymetry is usuallysimplified using traditional filtering techniques, that remove unwanted Fouriermodes. Principal Component Analysis (PCA) has been used in other fields toisolate small-scale structures from larger scale coherent features in a robustway, underpinned by a rigorous but simple mathematical framework. Here wepresent a method based on principal component analysis aimed towardssimplification of shorelines and bathymetry. We present the algorithm in detailand show simplified shorelines and bathymetry in the wider region around theNorth Sea. Finally, the methods are used in the context of unstructured meshgeneration aimed at tidal resource assessment simulations in the coastalregions around the UK.
Barral N, Knepley MG, Lange M, et al., Anisotropic mesh adaptation in Firedrake with PETSc DMPlex
Despite decades of research in this area, mesh adaptation capabilities arestill rarely found in numerical simulation software. We postulate that theprimary reason for this is lack of usability. Integrating mesh adaptation intoexisting software is difficult as non-trivial operators, such as error metricsand interpolation operators, are required, and integrating available adaptiveremeshers is not straightforward. Our approach presented here is to firstintegrate Pragmatic, an anisotropic mesh adaptation library, into DMPlex, aPETSc object that manages unstructured meshes and their interactions withPETSc's solvers and I/O routines. As PETSc is already widely used, this willmake anisotropic mesh adaptation available to a much larger community. As ademonstration of this we describe the integration of anisotropic meshadaptation into Firedrake, an automated Finite Element based system for theportable solution of partial differential equations which already uses PETScsolvers and I/O via DMPlex. We present a proof of concept of this integrationwith a three-dimensional advection test case.
Collins DS, Avdis A, Allison PA, et al., Tidal dynamics and mangrove carbon sequestration during the Oligo–Miocene in the South China Sea, Nature Communications, ISSN: 2041-1723
Modern mangroves are among the most carbon-rich biomes on Earth, but their long-term (≥106 yr) impact on the global carbon cycle is unknown. The extent, productivity and preservation of mangroves are controlled by the interplay of tectonics, global sea level and sedimentation, including tide, wave and fluvial processes. The impact of these processes on mangrove-bearing successions in the Oligo–Miocene of the South China Sea (SCS) is evaluated herein. Palaeogeographic reconstructions, palaeotidal modelling, and facies analysis suggest that elevated tidal range and bed shear stress optimised mangrove development along tide-influenced tropical coastlines. Preservation of mangrove organic carbon (OC) was promoted by high tectonic subsidence and fluvial sediment supply. Lithospheric storage of OC in peripheral SCS basins potentially exceeded 4000 Gt (equivalent to 2000 ppm of atmospheric CO2). These results highlight the crucial impact of tectonic and oceanographic processes on mangrove OC sequestration within the global carbon cycle on geological timescales.
Funke SW, Kramer SC, Piggott MD, Design optimisation and resource assessment for tidal-stream renewable energy farms using a new continuous turbine approach, Renewable Energy, ISSN: 1879-0682
This paper presents a new approach for optimising the design of tidal stream turbine farms. In this approach, the turbine farm is represented by a turbine density function that specifies the number of turbines per unit area and an associated continuous locally-enhanced bottom friction field. The farm design question is formulated as a mathematical optimisation problem constrained by the shallow water equations and solved with efficient, gradient-based optimisation methods. The resulting method is accurate, computationally efficient, allows complex installation constraints, and supports different goal quantities such as to maximise power or profit. The outputs of the optimisation are the optimal number of turbines, their location within the farm, the overall farm profit, the farm's power extraction, and the installation cost.We demonstrate the capabilities of the method on a validated numerical model of the Pentland Firth, Scotland. We optimise the design of four tidal farms simultaneously, as well as individually, and study how farms in close proximity may impact upon one another.
Funke SW, Farrell PE, Piggott MD, 2017, Reconstructing wave profiles from inundation data, Computer Methods in Applied Mechanics and Engineering, Vol: 322, Pages: 167-186, ISSN: 0045-7825
© 2017 Elsevier B.V.This paper applies variational data assimilation to inundation problems governed by the shallow water equations with wetting and drying. The objective of the assimilation is to recover an unknown time-varying wave profile at an open ocean boundary from inundation observations. This problem is solved with derivative-based optimisation and an adjoint wetting and drying scheme to efficiently compute sensitivity information. The capabilities of this approach are demonstrated on an idealised sloping beach setup in which the profile of an incoming wave is reconstructed from wet/dry interface observations. The method is robust to noise in the observations if a regularisation term is added to the optimisation objective. Finally, the method is applied to a laboratory experiment of the Hokkaido-Nansei-Oki tsunami, where the wave profile is reconstructed with a relative L∞ error of less than 1%.
Parkinson SD, Funke SW, Hill J, et al., 2017, Application of the adjoint approach to optimise the initial conditions of a turbidity current with the AdjointTurbidity 1.0 model, GEOSCIENTIFIC MODEL DEVELOPMENT, Vol: 10, Pages: 1051-1068, ISSN: 1991-959X
Smith R, 2017, Numerical modelling of tsunami generated by deformable submarine slides
Submarine slides can generate tsunami waves that cause significant damage and loss of life. Numerical modelling of submarine slide generated waves is complex and computationally challenging, but is useful to understand the nature of the waves that are generated, and identify the important factors in determining wave characteristics which in turn are used in risk assessments. In this work, the open-source, finite-element, unstructured mesh fluid dynamics framework Fluidity is used to simulate submarine slide tsunami using a number of different numerical approaches. First, three alternative approaches for simulating submarine slide acceleration, deformation and wave generation with full coupling between the slide and water in two dimensions are compared. Each approach is verified against benchmarks from experimental and other numerical studies, at different scales, for deformable submarine slides. There is good agreement to both laboratory results and other numerical models, both with a fixed mesh and a dynamically adaptive mesh, tracking important features of the slide geometry as the simulation progresses. Second, Fluidity is also used in a single-layer Bousinesq approximation in conjunction with a prescribed velocity boundary condition to model the propagation of slide tsunami in two and three dimensions. A new, efficient approach for submarine slide tsunami that accounts for slide dynamics and deformation is developed by imposing slide dynamics, derived from multi-material simulations. Two submarine slides are simulated in the Atlantic Ocean, and these generate waves up to 10 m high at the coast of the British Isles. Results indicate the largest waves are generated in the direction of slide motion. The lowest waves are generated perpendicular to the slide motion. The slide velocity and acceleration are the most important factors in determining wave height. Slides that deform generate higher waves than rigid slides, although this effect is of secondary importance f
Abolghasemi MA, Piggott MD, Spinneken J, et al., 2016, Simulating tidal turbines with multi-scale mesh optimisation techniques, JOURNAL OF FLUIDS AND STRUCTURES, Vol: 66, Pages: 69-90, ISSN: 0889-9746
Adam A, Buchan AG, Piggott MD, et al., 2016, Adaptive Haar wavelets for the angular discretisation of spectral wave models, JOURNAL OF COMPUTATIONAL PHYSICS, Vol: 305, Pages: 521-538, ISSN: 0021-9991
Avdis A, Jacobs CT, Mouradian SL, et al., 2016, Meshing ocean domains for coastal engineering applications, Pages: 480-492
As we continue to exploit and alter the coastal environment, the quantification of the potential impacts from planned coastal engineering projects, as well as the minimisation of any detrimental effects through design optimisation, are receiving increasing attention. Geophysical fluid dynamics simulations can provide valuable insight towards the mitigation and prevention of negative outcomes, and as such are routinely used for planning, operational and regulatory reasons. The ability to readily create high-quality computational meshes is critical to such modelling studies as it impacts on the accuracy, efficiency and reproducibility of the numerical results. To that end, most (coastal) ocean modelling packages offer tailored mesh generation utilities. Geographical Information Systems (GIS) offer an ideal framework within which to process data for use in the meshing of coastal regions. GIS have been designed specifically for the processing and analysis of geophysical data and are a popular tool in both the academic and industrial sectors. On the other hand Computer Aided Design (CAD) is the most appropriate tool for designing coastal structures and is usually the user interface to generic three-dimensional mesh generation frameworks. In this paper we combine GIS and CAD with a view towards mesh generation for an impact study of the proposed Swansea Bay Tidal Lagoon project within the Bristol Channel and Severn Estuary. We demonstrate in this work that GIS and CAD can be used in a complementary way to deliver unstructured mesh generation capabilities for coastal engineering applications.
Culley DM, Funke SW, Kramer SC, et al., 2016, Integration of cost modelling within the micro-siting design optimisation of tidal turbine arrays, RENEWABLE ENERGY, Vol: 85, Pages: 215-227, ISSN: 0960-1481
Funke SW, Kramer SC, Piggott MD, 2016, Design optimisation and resource assessment for tidal-stream renewable energy farms using a new continuous turbine approach, Publisher: PERGAMON-ELSEVIER SCIENCE LTD
Jacobs CT, Piggott MD, Kramer SC, et al., 2016, On the validity of tidal turbine array configurations obtained from steady-state adjoint optimisation, ECCOMAS Congress 2016 - Proceedings of the 7th European Congress on Computational Methods in Applied Sciences and Engineering, Vol: 4, Pages: 8247-8261
Extracting the optimal amount of power from an array of tidal turbines requires an intricate understanding of tidal dynamics and the effects of turbine placement on the local and regional scale flow. Numerical models have contributed significantly towards this understanding, and more recently, adjoint-based modelling has been employed to optimise the positioning of the turbines in an array in an automated way and improve on simple man-made configurations (e.g. structured grids of turbines) . Adjoint-based optimisation of high-resolution and ideally 3D transient models is generally a very computationally expensive problem. Multiple approaches are therefore used in practice to obtain feasible runtimes: using high viscosity values to obtain a steady-state solution, or a sequence of steady-state solutions for "time-varying" setups; limiting the number of adjoint computations; or reformulating the problem to allow for coarser mesh resolution to make it feasible for resources assessment (e.g. , ). However, such compromises may affect the reliability of the modelled turbines, their wakes and interactions, and thus bring into question the validity of the computed optimal turbine positions. This work considers a suite of idealised simulations of flow past tidal turbine arrays in a two-dimensional channel. It compares four regular array configurations, detailed by Divett et al. , with the configuration found through adjoint optimisation in a steady-state, high-viscosity setup. The optimised configuration produces considerably more power than the other configurations (approximately 40% more than the best man-made configuration). The same configurations are then used to produce a suite of transient simulations that do not use constant high-viscosity, and instead use large eddy simulation (LES) to parameterise the resulting turbulent structures. All simulations are performed using OpenTidalFarm . It is shown that the 'low background viscosity'/LES simula
Kramer SC, Piggott MD, 2016, A correction to the enhanced bottom drag parameterisation of tidal turbines, Publisher: PERGAMON-ELSEVIER SCIENCE LTD
Kramer SC, Piggott MD, 2016, A correction to the enhanced bottom drag parameterisation of tidal turbines, Renewable Energy, Vol: 92, Pages: 385-396, ISSN: 1879-0682
Hydrodynamic modelling is an important tool for the development of tidalstream energy projects. Many hydrodynamic models incorporate the effect oftidal turbines through an enhanced bottom drag. In this paper we show thatalthough for coarse grid resolutions (kilometre scale) the resulting force exertedon the flow agrees well with the theoretical value, the force starts decreasingwith decreasing grid sizes when these become smaller than the length scale ofthe wake recovery. This is because the assumption that the upstream velocitycan be approximated by the local model velocity, is no longer valid. Using linearmomentum actuator disc theory however, we derive a relationship between thesetwo velocities and formulate a correction to the enhanced bottom drag formulationthat consistently applies a force that remains close to the theoretical value,for all grid sizes down to the turbine scale. In addition, a better understandingof the relation between the model, upstream, and actual turbine velocity, aspredicted by actuator disc theory, leads to an improved estimate of the usefullyextractable energy. We show how the corrections can be applied (demonstratedhere for the models MIKE 21 and Fluidity) by a simple modification of the dragcoefficient.
Mouradian, Avdis A, Piggott M, et al., 2016, TELEMAC model archive: Integrating open-source tools for the management and visualisation of model data, 23rd TELEMAC-MASCARET User Conference (TUC-2016)
Nunez Rattia JM, Percival JR, Yeager B, et al., 2016, Numerical simulation of scour below pipelines using flexible mesh methods, Pages: 101-108
© 2016 Taylor & Francis Group, London.Evaluating bed morphological structure and evolution (specifically the scoured bed level) accurately using numerical models is critical for analyses of the stability of many marine structures. This paper discusses the performance of an implementation within Fluidity, an open source, general purpose, Computational Fluid Dynamics (CFD) code, capable of handling arbitrary multi-scale unstructured tetrahedral meshes and including algorithms to perform dynamic anisotropic mesh adaptivity. The flexibility over mesh structure and resolution that these capabilities provide makes it potentially highly suitable for coupling the structural scale with larger scale ocean dynamics. In this very preliminary study the solver approach is demonstrated for an idealised scenario. Discontinuous Galerkin finite-element (DG-FEM) based discretisation methods have been used for the hydrodynamics and morphological calculations, and automatic mesh deformation has been utilised to account for bed evolution changes while preserving the validity and quality of the mesh. In future work, the solver will be used in three-dimensional impinging jet and other industrial and environmental scour studies.
Piggott MD, 2016, Thetis
Finite element flow solver for simulating coastal and estuarine flows.
Quattrocchi G, Gorman GJ, Piggott MD, et al., 2016, M2, overtides and compound tides generation in the Strait of Messina: the response of a non-hydrostatic, finite-element ocean model, JOURNAL OF COASTAL RESEARCH, Pages: 657-661, ISSN: 0749-0208
Smith RC, Hill J, Collins GS, et al., 2016, Comparing approaches for numerical modelling of tsunami generation by deformable submarine slides, OCEAN MODELLING, Vol: 100, Pages: 125-140, ISSN: 1463-5003
Vire A, Spinneken J, Piggott MD, et al., 2016, Application of the immersed-body method to simulate wave-structure interactions, EUROPEAN JOURNAL OF MECHANICS B-FLUIDS, Vol: 55, Pages: 330-339, ISSN: 0997-7546
Jacobs CT, Avdis A, Mouradian SL, et al., 2015, Integrating research data management into geographical information systems, CEUR Workshop Proceedings, Vol: 1529, Pages: 7-17, ISSN: 1613-0073
Ocean modelling requires the production of high-fidelity com-putational meshes upon which to solve the equations of motion. The production of such meshes by hand is often infeasible, considering the complexity of the bathymetry and coastlines. The use of Geographical Information Systems (GIS) is therefore a key component to discretising the region of interest and producing a mesh appropriate to resolve the dynamics. However, all data abociated with the production of a mesh must be provided in order to contribute to the overall recomputability of the subsequent simulation. This work presents the integration of re-search data management in QMesh, a tool for generating meshes using GIS. The tool uses the PyRDM library to provide a quick and easy way for scientists to publish meshes, and all data required to regenerate them, to persistent online repositories. These repositories are abigned unique identifiers to enable proper citation of the meshes in journal articles.
Jacobs CT, Avdis A, Mouradian SL, et al., 2015, Integrating Research Data Management into Geographical Information Systems, http://ceur-ws.org/Vol-1529/, 5th International Workshop on Semantic Digital Archives (SDA 2015), Pages: 7-17
Ocean modelling requires the production of high-fidelity computational meshes upon which to solve the equations of motion. The production of such meshes by hand is often infeasible, considering the complexity of the bathymetry and coastlines. The use of Geographical Information Systems (GIS) is therefore a key component to discretising the region of interest and producing a mesh appropriate to resolve the dynamics. However, all data associated with the production of a mesh must be provided in order to contribute to the overall recomputability of the subsequent simulation. This work presents the integration of research data management in QMesh, a tool for generating meshes using GIS. The tool uses the PyRDM library to provide a quick and easy way for scientists to publish meshes, and all data required to regenerate them, to persistent online repositories. These repositories are assigned unique identifiers to enable proper citation of the meshes in journal articles.
Jacobs CT, Goldin TJ, Collins GS, et al., 2015, An improved quantitative measure of the tendency for volcanic ash plumes to form in water: implications for the deposition of marine ash beds, JOURNAL OF VOLCANOLOGY AND GEOTHERMAL RESEARCH, Vol: 290, Pages: 114-124, ISSN: 0377-0273
Jacobs CT, Piggott MD, 2015, Firedrake-Fluids v0.1: numerical modelling of shallow water flows using an automated solution framework, GEOSCIENTIFIC MODEL DEVELOPMENT, Vol: 8, Pages: 533-547, ISSN: 1991-959X
Jordan JR, Kimura S, Holland PR, et al., 2015, On the Conditional Frazil Ice Instability in Seawater, JOURNAL OF PHYSICAL OCEANOGRAPHY, Vol: 45, Pages: 1121-1138, ISSN: 0022-3670
Martin-Short R, Hill J, Kramer SC, et al., 2015, .Tidal resource extraction in the Pentland Firth, UK: Potential impacts on flow regime and sediment transport in the Inner Sound of Stroma, RENEWABLE ENERGY, Vol: 76, Pages: 596-607, ISSN: 0960-1481
Piggott MD, 2015, OpenTidalFarm
OpenTidalFarm is an open-source software for simulating and optimising tidal turbine farms.The positioning of the turbines in a tidal farm is a crucial decision. Simulations show that the optimal positioning can increase the power generation of the farm by up to 50% and can therefore determine the viability of a project. However, finding the optimal layout is a difficult process due to the complex flow interactions. OpenTidalFarm solves this problem by applying an efficient optimisation algorithm onto a accurate flow prediction model.
Barnett GL, Funke SW, Piggott MD, 2014, Hybrid global-local optimisation algorithms for the layout design of tidal turbine arrays
Tidal stream power generation represents a promising source of renewableenergy. In order to extract an economically useful amount of power, tens tohundreds of tidal turbines need to be placed within an array. The layout ofthese turbines can have a significant impact on the power extracted and henceon the viability of the site. Funke et al. formulated the question of the bestturbine layout as an optimisation problem constrained by the shallow waterequations and solved it using a local, gradient-based optimisation algorithm.Given the local nature of this approach, the question arises of how optimal thelayouts actually are. This becomes particularly important for scenarios withcomplex bathymetry and layout constraints, both of which typically introducelocally optimal layouts. Optimisation algorithms which find the global optimagenerally require orders of magnitude more iterations than local optimisationalgorithms and are thus infeasible in combination with an expensive flow model.This paper presents an analytical wake model to act as an efficient proxy tothe shallow water model. Based upon this, a hybrid global-local two-stageoptimisation approach is presented in which turbine layouts are first optimisedwith the analytical wake model via a global optimisation algorithm, and furtheroptimised with the shallow water model via a local gradient-based optimisationalgorithm. This procedure is applied to a number of idealised cases and a morerealistic case with complex bathymetry in the Pentland Firth, Scotland. It isshown that in cases where bathymetry is considered, the two-stage optimisationprocedure is able to improve the power extracted from the array by as much as25% compared to local optimisation for idealised scenarios and by as much as12% for the more realistic Pentland Firth scenario whilst in many casesreducing the overall computation time by approximately 35%.
This data is extracted from the Web of Science and reproduced under a licence from Thomson Reuters. You may not copy or re-distribute this data in whole or in part without the written consent of the Science business of Thomson Reuters.