Abstract

Brash ice is the accumulation of floating ice made up of blocks no larger than two meters across. Navigation in brash ice is becoming more usual as new navigation routes are being opened in the Artic regions. This navigation brings new concerns regarding the interaction of ice blocks with the ship. This work presents recent advances towards the development of a computational model for simulation of this navigation condition including the interaction among the ship and the ice blocks. The computational tool developed in this work is based on the coupling of a Semi-Lagrangian Particle Finite Element Method (SL-PFEM) with a multi rigid-body dynamics tool. The Particle Finite Element Method (Idelsohn et al. 2004) is a versatile framework for the analysis of fluid-structure interaction problems. The PFEM combines Lagrangian particle-based techniques with the advantage of the integral formulation of the Finite Element Method (FEM). It has been shown (Idelsohn et al. 2014 and  Nadukandi et al. 2017) to successfully simulate a wide variety of complex engineering problems, e.g. free-surface/multi-fluid flows with violent interface motions, multi-fluid mixing and buoyancy-driven segregation problems etc. The latest development within the framework of the PFEM is the X-IVAS (eXplicit Integration along the Velocity and Acceleration Streamlines) scheme. It is a semi-implicit scheme built over a Semi-Lagrangian (SL) formulation of the PFEM. In this work, the SL-PFEM model has been coupled with a multibody dynamics solver, able to handle the interactions between thousands of bodies, representing the different ice blocks. The interaction between the fluid flow and the ice blocks is taking into account by enriching the finite element space at the boundaries of the different blocks. This work is part of the research project NICESHIP sponsored by the U.S. Office of Naval Research under Grant N62909-16-1-2236.

Abstract

In this work a finite element method (FEM) is proposed to solve the problem of estimating the added resistance of a ship in waves in the time domain and using unstructured meshes. Two different schemes are used to integrate the corresponding free surface kinematic and dynamic boundary conditions: the first one based on streamlines integration (SLI); and the second one based on the Streamline-Upwind Petrov-Galerking (SUPG) stabilization. The proposed numerical schemes have been validated in different test cases, including towing tank tests with monochromatic waves. The results obtained in this work show the suitability of the present method to estimate added resistance in waves in a computationally affordable manner.

Abstract

Brash ice is the accumulation of floating ice made up of blocks no larger than two meters across. Navigation in brash ice is becoming more usual as new navigation routes are being opened in the Artic regions. This navigation brings new concerns regarding the interaction of ice blocks with the ship. This work presents recent advances towards the development of a computational model for simulation of this navigation condition including the interaction among the ship and the ice blocks. The computational tool developed in this work is based on the coupling of a Semi-Lagrangian Particle Finite Element Method (SL-PFEM) with a multi rigid-body dynamics tool. The Particle Finite Element Method (Idelsohn et al. 2004) is a versatile framework for the analysis of fluid-structure interaction problems. The PFEM combines Lagrangian particle-based techniques with the advantage of the integral formulation of the Finite Element Method (FEM). It has been shown (Idelsohn et al. 2014 and  Nadukandi et al. 2017) to successfully simulate a wide variety of complex engineering problems, e.g. free-surface/multi-fluid flows with violent interface motions, multi-fluid mixing and buoyancy-driven segregation problems etc. The latest development within the framework of the PFEM is the X-IVAS (eXplicit Integration along the Velocity and Acceleration Streamlines) scheme. It is a semi-implicit scheme built over a Semi-Lagrangian (SL) formulation of the PFEM. In this work, the SL-PFEM model has been coupled with a multibody dynamics solver, able to handle the interactions between thousands of bodies, representing the different ice blocks. The interaction between the fluid flow and the ice blocks is taking into account by enriching the finite element space at the boundaries of the different blocks. This work is part of the research project NICESHIP sponsored by the U.S. Office of Naval Research under Grant N62909-16-1-2236.

Abstract

The objective of this thesis is the research on numerical algorithms to develop numerical tools to simulate seakeeping problems as well as wave resistance problems of ships and floating structures. The first tool developed is a wave diffraction-radiation solver. It is based on the finite element method (FEM) in order to solve the Laplace equation, as well as numerical schemes based on FEM, streamline integration, and finite difference method tailored for solving the free surface boundary condition. It has been developed numerical tools to solve solid body dynamics of multibody systems with body links across them. This tool has been integrated with the wave diffraction-radiation solver to solve wave-body interaction problems. Also it has been tailored coupling algorithms with other numerical tools in order to solve multi-physics problems. In particular, it has been performed coupling with a MEF structural solver to solve fluid-structure interaction problems, with a mooring solver, and with a solver capable of simulating internal flows in tanks to solve couple seakeeping-sloshing problems. Numerical simulations have been carried out to validate and verify the developed algorithms, as well as to analyze case studies in the areas of marine engineering, offshore engineering, and offshore renewable energy.

Abstract

Wave-structure interaction is a topic of great interest in naval and offshore engineering. This interest is growing in the last years due to the boost given by the development of the marine renewable energy field. In this context the development of an efficient time-domain coupled waves-structure solver is a main request from the industry. Up to date the numerical seakeeping simulation has been mostly carried out using the frequency domain. The reason might be that the computational cost of time domain simulations were too high and computational time was too large. Moreover assumptions like linear waves and the harmonic nature of water waves made the frequency domain to be the right choice. However nowadays computing capabilities make possible to carry out numerical simulations in the time domain in a reasonable time, with the advantage of making easier the introduction of non-linearities into the algorithm and therefore coupling with other phenomena. This presentation shows the work of the authors in developing a time-domain unstructured Finite Element Method (FEM) algorithm for analysis of coupled wave-structure interaction. For this purpose, a new diffraction-radiation solver using the FEM was developed. The solver has been implemented in GPU, using CUDA architecture. The speed-up obtained ranges from 5 to 10 times compare to the implementation in a standard CPU with a conjugate gradient and ILU preconditioner. The seakeeping analysis tool has been integrated within a coupled waves-structure analysis tool. The coupling algorithm is based on a partitioned iterative algorithm, using an interpolation library able to communicate pressure forces and displacements of the structure at memory level. Furthermore, an innovative wetting and drying scheme able to improve the evaluation of the water action on the structure. The accuracy of the new seakeeping formulation for analysis of waves and floating structures interaction has been verified in different validation cases and practical applications. 

Abstract

This presentation shows the recent work of the authors in the development of a time-domain FEM model for evaluation of the seal dynamics of a surface effect ship. The fluid solver developed for this purpose, uses a potential flow approach along with a stream-line integration of the free surface. The presentation focuses on the free surface-structure algorithm that has been developed to allow the simulation of the complex and highly dynamic behavior of the seals in the interface between the air cushion, and the water. The developed fluid-structure interaction solver is based, on one side, on an implicit iteration algorithm, communicating pressure forces and displacements of the seals at memory level and, on the other side, on an innovative wetting and drying scheme able to predict the water action on the seals. The stability of the iterative scheme is improved by means of relaxation, and the convergence is accelerated using Aitken’s method. Several validations against experimental results have been carried out to demonstrate the developed algorithm.

Abstract

A Wave Adaptive Modular Vessel (WAM-V®) is a new class of ship that uses inflatable flexible hulls to conform to the surface of the water. It is similar in design to a catamaran, in that it has a twin hull design and no keel. However, the superstructure is not rigidly attached to the hulls; it uses shock absorbers and ball joints to articulate the vessel, which allows WAM-V to conform to the surface of the water while mitigating the stresses transmitted to the structure. Moreover, the inflatable hulls help to absorb the high frequency wave-loads. These features allow WAM-V to travel efficiently with low wave resistance in rough seas, by surfing on top of the waves rather than cut through them. The objective of the WAM-V is to be a lightweight watercraft capable of moving fast and efficiently on the surface of the sea. WAM-Vs are designed to allow for a variety of applications for either manned or unmanned operations and can be built in different lengths to match specific services. This presentation shows part of the work done in the project ‘Advanced Numerical Simulation and Performance Evaluation of WAM-V ® in Spray Generating Conditions’ developed by the International Center for Numerical Methods in Engineering (CIMNE) under Navy Grant N62909-12-1-7101 issued by the Office of Naval Research Global. The scope of that project included the performance analysis of the WAM-V in waves, taking into account the flexibility of the ship hulls, using fluid-structure interaction computational models. However, the focus of this paper is one of the primary concerns of that project; the development of a computational model for simulation of the WAM-V under spray generating conditions. In this regards, the final goal was to develop and demonstrate a computational engineering solver that could be used to design strategies to reduce the spray generation of the vessel.

Abstract

The aim of this work is to be able to cope with complex sloshing-seakeeping problems in an effective manner. For this purpose, two different solvers will be integrated into one coupled tool to take advantage of their characteristics. The Particle Finite Element Method is a versatile framework for the analysis of fluid-structure interaction problems. The latest development within the framework of the PFEM is the X-IVAS (eXplicit Integration along the Velocity and Acceleration Streamlines) scheme. It is a semi-implicit scheme built over a Semi-Lagrangian (SL) formulation. This new scheme was named PFEM-2 and will be used in this work to solve the fluid dynamics (sloshing) inside the tanks. The PFEM-2 will be coupled in the time domain with SeaFEM, a solver developed for seakeeping problems. SeaFEM solves the second-order diffraction-radiation equations by using the Finite Element Method (FEM) and will be used in this work to simulate the interaction between waves and a floating body. The coupling of the two tools will be accomplished by using an effective coupling algorithm and a communication technique that allows simulations to be computed without affecting the global compute time and the accuracy of the solvers. This new tool has been validated against experimental benchmarks carried out in a model basin at model scale.

Abstract

A finite element method for the solution of the up-to-second-order wave diffraction-radiation problem in the time-domain is proposed. The solver has been validated against experimental data available for the HiPRWind semisubmersible platform (designed for floating wind turbines). To perform the validation, the wave diffraction-radiation solver is coupled to a body dynamics and mooring solvers in the time-domain. The HiPRWind movements and mooring forces have been compared for a large number of test cases, including decay tests, bichromatic, and irregular waves. Good agreement has been found for both, body movements and mooring forces.

Abstract

Being capable of predicting seakeeping capabilities in the time domain is of great interest for the marine and offshore industry. However, most computer programs used in the industry work in the frequency domain, with the subsequent limitation in the accuracy of their model predictions. The main objective of this work is the development of a time domain solver based on the finite element method capable of solving multi-body seakeeping problems on unstructured meshes. In order to achieve such an objective, several techniques are combined: the use of an efficient algorithm for the free surface boundary conditions, the use of deflated conjugate gradients, and the use of a graphic processing unit for speeding up the linear solver. The results obtained results by the developed model showed good agreement with analytical solutions, experimental data for an actual offshore structure model, as well as numerical solutions obtained by other numerical method. Also, a simulation with sixteen floating bodies was carried out in a affordable CPU time, showing the potential of this approach for multi-body simulation.