We present a high-order nodal spectral element method for the two-dimensional simulation of nonlinear water waves. The model is based on the mixed Eulerian–Lagrangian (MEL) method. Wave interaction with fixed truncated structures is handled using unstructured meshes consisting of high-order iso-parametric quadrilateral/triangular elements to represent the body surfaces as well as the free surface elevation. A numerical eigenvalue analysis highlights that using a thin top layer of quadrilateral elements circumvents the general instability problem associated with the use of asymmetric mesh topology.We demonstrate how to obtain a robust MEL scheme for highly nonlinear waves using an efficient combination of (i) global L2 projection without quadrature errors, (ii) mild modal filtering and (iii) a combination of local and global re-meshing techniques. Numerical experiments for strongly nonlinear waves are presented. The experiments demonstrate that the spectral element model provides excellent accuracy in prediction of nonlinear and dispersive wave propagation. The model is also shown to accurately capture the interaction between solitary waves and fixed submerged and surface-piercing bodies. The wave motion and the wave-induced loads compare well to experimental and computational results from the literature.
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.
We present recent progress on the development of a new fully nonlinear potential flow (FNPF) model for estimation of nonlinear wave-body interactions based on a stabilized unstructured spectral element method (SEM). We introduce new proof-of-concepts for forced nonlinear wave-body interaction in two spatial dimensions to establish the methodology in the SEM setting utilizing dynamically adapted unstructured meshes. The numerical method behind the proposed methodology is described in some detail and numerical experiments on the forced motion of (i) surface piercing and (ii) submerged bodies are presented.
This book introduces a novel model to explain how the co-design and co-delivery of ocean science knowledge and solutions is influenced by ocean stakeholders with asymmetric power and resources, policy incentives and ocean conflict, ocean narratives, different knowledge systems, security concerns, principles, formal and informal rules, and communication competencies. Using the International Collaboration in Ocean Science model as a basis, the book advances with three lines of inquiry: ontological security of ocean science participants, the Ocean Decade and human well-being, and strategic narratives about international collaboration in ocean science. Through these, Carolijn van Noort shows the enabling and constraining conditions of co-creating ocean knowledge and solutions. Theoretically novel, the book provides a compelling framework for scholars to study ocean science collaboration.
The depth-integrated shallow water equations are frequently used for simulating geophysical flows, such as storm-surges, tsunamis and river flooding. In this paper a parallel shallow water solver using an unstructured high-order discontinuous Galerkin method is presented. The spatial discretization of the model is based on the Nektar++ spectral/hp library and the model is numerically shown to exhibit the expected exponential convergence. The parallelism of the model has been achieved within the Cactus Framework. The model has so far been executed successfully on up to 128 cores and it is shown that both weak and strong scaling are largely independent of the spatial order of the scheme. Results are also presented for the wave flume interaction with five upright cylinders.
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, blindly simulate 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.
Temperature data 1900–2010 from meteorological stations across the world have been analyzed and it has been found that all land areas generally have two different valid temperature trends. Coastal stations and hill stations facing ocean winds are normally more warm-trended than the valley stations that are sheltered from dominant oceans winds.
Thus, we found that in any area with variation in the topography, we can divide the stations into the more warm trended ocean air-affected stations, and the more cold-trended ocean air-sheltered stations. We find that the distinction between ocean air-affected and ocean air-sheltered stations can be used to identify the influence of the oceans on land surface. We can then use this knowledge as a tool to better study climate variability on the land surface without the moderating effects of the ocean.
We find a lack of warming in the ocean air sheltered temperature data – with less impact of ocean temperature trends – after 1950. The lack of warming in the ocean air sheltered temperature trends after 1950 should be considered when evaluating the climatic effects of changes in the Earth’s atmospheric trace amounts of greenhouse gasses as well as variations in solar conditions.
Continuous inspection and mapping of the seabed allows for monitoring the impact of anthropogenic activities on benthic ecosystems. Compared to traditional manual assessment methods which are impractical at scale, computer vision holds great potential for widespread and long-term monitoring.
We deploy an underwater remotely operated vehicle (ROV) in Jammer Bay, a heavily fished area in the Greater North Sea, and capture videos of the seabed for habitat classification. The collected JAMBO dataset is inherently ambiguous: water in the bay is typically turbid which degrades visibility and makes habitats more difficult to identify. To capture the uncertainties involved in manual visual inspection, we employ multiple annotators to classify the same set of images and analyze time spent per annotation, the extent to which annotators agree, and more.
We then evaluate the potential of vision foundation models (DINO, OpenCLIP, BioCLIP) for automating image-based benthic habitat classification. We find that despite ambiguity in the dataset, a well chosen pre-trained feature extractor with linear probing can match the performance of manual annotators when evaluated in known locations. However, generalization across time and place is an important challenge.
A moving static pressure distribution is commonly used to simulate a travelling ship. However, the ship movement changes the fluid velocity around the hull, inducing pressures on the hull surface that are no longer equal to the static pressure. Therefore, we introduce a dynamic pressure correction strategy, which can accurately simulate the impact of the ship movement on the hull-surface pressure and preserve the desired hull shape under both stationary and transient conditions. The strategy is applied to a high-order spectral model and used to investigate ship-induced waves and wave resistance over a both flat and variable topography. We explore various parameters in our study, including the average water depth to ship draft ratio (h(0)/d), the channel width to ship width ratio (W/B), the Froude number (Fr-0 = U/root gh(0)) and variations in bathymetric slope. Compared with experiments on a flat bottom, the numerical results with dynamic correction show better accuracy in the simulation of ship-induced waves and wave resistance than those obtained using a static pressure distribution. The correlation coefficient for wake waves between the numerical and experimental results is improved by approximately 0.25 with the dynamic correction strategy. The amplitude and wavelength of ship-induced mini-tsunamis over a variable topography are found to be reduced when employing a dynamic correction compared with a static pressure distribution, and this effect becomes more pronounced with higher Froude number. The static pressure approach is shown to allow large deformations of the desired hull shape and changes in ship volume which are responsible for the different wave patterns from the two approaches.
Simulating the free decay motion and wave radiation from a heaving semi-submerged sphere poses significant computational challenges due to its three-dimensional complexity. By leveraging axisymmetry, we reduce the problem to a two-dimensional simulation, significantly decreasing computational demands while maintaining accuracy. In this paper, we exploit axisymmetry to perform a large ensemble of Computational Fluid Dynamics (CFDs) simulations, aiming to evaluate and maximize both accuracy and efficiency, using the Reynolds Averaged Navier–Stokes (RANS) solver interFOAM, in the opensource finite volume CFD software OpenFOAM. Validated against highly accurate experimental data, extensive parametric studies are conducted, previously limited by computational constraints, which facilitate the refinement of simulation setups. More than 50 iterations of the same heaving sphere simulation are performed, informing efficient trade-offs between computational cost and accuracy across various simulation parameters and mesh configurations. Ultimately, by employing axisymmetry, this research contributes to the development of more accurate and efficient numerical modeling in ocean engineering.