Search

En la actualidad las más avanzadas técnicas para diseño naval ya no se restringen sólo a proyectos de alto coste; si no que también ahora son de aplicación en el campo de las embarcaciones de recreo y turísticas, dando con ello un valor añadido al proyecto final. La utilización de un código CFD para el mayor conocimiento del comportamiento hidrodinámico de las embarcaciones de recreo y turísticas es de una gran utilidad debido a las complejidades que en el diseño de estas se presentan. Estas herramientas nos van a proporcionar una serie de datos como sustentación lateral en las embarcaciones de vela o interacción entre cascos en catamaranes con deriva lateral, que de otra forma nunca pudieran conocerse. Los resultados hidrodinámicos obtenidos nos sirven para comparar las diferentes alternativas, y para además obtener unos valores que directamente son usados para la predicción de velocidad de la embarcación, que más tarde se verán refrendadas por las pruebas de mar. La utilización de técnicas CFD, basadas en el método de los elementos finitos, para la resolución de las ecuaciones de Reynolds (RANSE), incluyendo efectos de superficie libre, hace que todo esto sea posible. En este trabajo se presentan diversos resultados de aplicación del código CFD SHYNE a varias embarcaciones deportivas y de recreo. SHYNE está basado en la técnica del Cálculo Finitesimal para la integración de las ecuaciones de Reynolds.

A stabilized semi-implicit fractional step algorithm based on the finite element method for solving ship wave problems using unstructured meshes is presented. The stabilized governing equations for the viscous incompressible fluid and the free surface are derived at a differential level via a finite calculus procedure. This allows us to obtain a stabilized numerical solution scheme. Some particular aspects of the problem solution, such as the mesh updating procedure and the transom stern treatment, are presented. Examples of the efficiency of the semi-implicit algorithm for the analysis of ship hydrodynamics problems are presented.

We present a general formulation for incompressible fluid flow analysis using the finite element method (FEM). The standard Eulerian formulation is described first. The necessary stabilization for dealing with convective effects and the incompressibility condition are introduced via the finite calculus (FIC) method. A simple extension of the fluid flow equations to an arbitrary Lagrangian–Eulerian (ALE) frame adequate for treating fluid–structure interaction problems is briefly presented. A fully Lagrangian formulation called the particle finite element method (PFEM) is also described. The PFEM is particularly attractive for fluid–structure interaction problems involving large motions of the free surface and breaking waves. Examples of application of the Eulerian, the ALE and the fully lagrangian PFEM formulations are presented.

A stabilized semi-implicit fractional step finite element method for solving coupled fluid-structure interaction problems involving free surface waves is presented. The stabilized governing equations for the viscous fluid and the free surface are derived at a differential level via a finite increment calculus procedure. A mesh updating technique based on solving a fictitious elastic problem on the moving mesh is described. Examples of the efficiency of the stabilized semi-implicit algorithm for the analysis of fluid-structure interaction problems in totally or partially submerged bodies is presented.

This presentation shows the recent work of the CIMNE in the maritime transport field. It was given at the Conference on Computation and Big Data in Transport (CM3-2017) held in November 22 – 23, 2017.

Aula FNB-CIMNE was created on 2002 through an agreement between the Faculty of Nautical Studies of Barcelona (FNB) and the International Center for Numerical Methods in Engineering (CIMNE) aiming to promote the development and application of computational methods for the solution of naval & maritime engineering problems.

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.

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 showsthe 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.&nbsp

A new Petrov–Galerkin (PG) method involving two parameters, namelyα1andα2, is presented, which yields the following schemes on rectangular meshes: (i) a compact stencil obtained by the linear interpolation of the Galerkin FEM and the classical central finite difference method (FDM), should the parameters be equal, that is,α1 = α2 = α; and (ii) the nonstandard compact stencil presented in (Int. J. Numer. Meth. Engng2011; 86:18–46) for the Helmholtz equation if the parameters are distinct, that is,α1 ≠ α2. The nonstandard compact stencil is obtained by taking the linear interpolation of the diffusive terms (specified byα1) and the mass terms (specified byα2) that appear in the stencils obtained by the standard Galerkin FEM and the classical central FDM, respectively. On square meshes, these two schemes were shown to provide solutions to the Helmholtz equation that have a dispersion accuracy of fourth and sixth order, respectively (Int. J. Numer. Meth. Engng2011; 86:18–46). The objective of this paper is to study the performance of this PG method for the Helmholtz equation using nonuniform meshes and the treatment of natural boundary conditions.

A multidimensional extension of the HRPG method using the lowest order block finite elements is presented. First, we design a nondimensional element number that quantifies the characteristic layers which are found only in higher dimensions. This is done by matching the width of the characteristic layers to the width of the parabolic layers found for a fictitious 1D reaction–diffusion problem. The nondimensional element number is then defined using this fictitious reaction coefficient, the diffusion coefficient and an appropriate element size. Next, we introduce anisotropic element length vectorsliand the stabilization parametersαi,βiare calculated along theseli. Except for the modification to include the new dimensionless number that quantifies the characteristic layers, the definitions ofαi,βiare a direct extension of their counterparts in 1D. Usingαi,βiandli, objective characteristic tensors associated with the HRPG method are defined. The numerical artifacts across the characteristic layers are manifested as the Gibbs phenomenon. Hence, we treat them just like the artifacts formed across the parabolic layers in the reaction-dominant case. Several 2D examples are presented that support the design objective—stabilization with high-resolution

We propose a fourth‐order compact scheme on structured meshes for the Helmholtz equation given byR(φ):=f(x)+Δφ+ξ2φ=0. The scheme consists of taking the alpha‐interpolation of the Galerkin finite element method and the classical central finite difference method. In 1D, this scheme is identical to the alpha‐interpolation method (J. Comput. Appl. Math.1982;8(1):15–19) and in 2D making the choice α=0.5 we recover the generalized fourth‐order compact Padé approximation (J. Comput. Phys.1995;119:252–270;Comput. Meth. Appl. Mech. Engrg1998;163:343–358) (therein using the parameter γ=2). We follow (SIAM Rev.2000;42(3):451–484;Comput. Meth. Appl. Mech. Engrg1995;128:325–359) for the analysis of this scheme and its performance on square meshes is compared with that of the quasi‐stabilized FEM (Comput. Meth. Appl. Mech. Engrg1995;128:325–359). In particular, we show that the relative phase error of the numerical solution and the local truncation error of this scheme for plane wave solutions diminish at the rateO((ξℓ)4), where ξ, ℓ represent the wavenumber and the mesh size, respectively. An expression for the parameter α is given that minimizes the maximum relative phase error in a sense that will be explained in Section 4.5. Convergence studies of the error in theL2norm, theH1semi‐norm and thel∞Euclidean norm are done and the pollution effect is found to be small. Copyright © 2010 John Wiley & Sons, Ltd.

We present the design of a high-resolution Petrov–Galerkin (HRPG) method using linear finite elements for the problem defined by the residual R (phi):= partial phi /partial t + upartial phi /partial x - kpartial^2 phi /partial x^2 +s phi -f where k; s > o =0. The structure of the method in 1D is identical to the consistent approximate upwind Petrov– Galerkin (CAU/PG) method [A.C. Galeão, E.G. Dutra do Carmo, A consistent approximate upwind Petrov– Galerkin method for the convection-dominated problems, Comput. Methods Appl. Mech. Engrg. 68 (1988) 83–95] except for the definitions of the stabilization parameters. Such a structure may also be attained via the finite-calculus (FIC) procedure [E. Oñate, Derivation of stabilized equations for numerical solution of advective–diffusive transport and fluid flow problems, Comput. Methods Appl. Mech. Engrg. 151 (1998) 233–265; E. Oñate, J. Miquel, G. Hauke, Stabilized formulation for the advection–diffusion– absorption equation using finite-calculus and linear finite elements, Comput. Methods Appl. Mech. Engrg. 195 (2006) 3926–3946] by an appropriate definition of the characteristic length. The prefix ‘high-resolution’ is used here in the sense popularized by Harten, i.e. second order accuracy for smooth/regular regimes and good shock-capturing in nonregular regimes. The design procedure embarks on the problem of circumventing the Gibbs phenomenon observed in L2-projections. Next we study the conditions on the stabilization parameters to circumvent the global oscillations due to the convective term. A conjuncture of the two results is made to deal with the problem at hand that is usually plagued by Gibbs, global and dispersive oscillations in the numerical solution. It is shown that the method indeed reproduces stabilized high-resolution numerical solutions for a wide range of values of u; k; s and f . Finally, some remarks are made on the extension of the HRPG method to multidimensions.

In the present work a novel approach has been developed for the resolution of this problem of analysis of the movement of roll of a ship. The methodology used for this is based on the modification of the differential equations of the fluid dynamics (RANSE equations), including the movement of the free surface, by applying the finite calculus method. The modified equations are solved using an implicit predictor-corrector scheme and the finite element method (FEM). This resolution scheme is considered optimal for these types of problems, both in accuracy and in calculation time. As an example, it should be noted that for typical problems (more than 1,000,000 finite elements) less than 4 hours of CPU are required to solve several cycles of ship movement (Pentium IV). This allows the analysis of an array of tests in a few days or even hours, having a computer network or sufficient computing power.

El “Innovative Naval Prototype Transformable Craft” (T-Craft) es un nuevo concepto de buque desarrollado por la marina de los Estados Unidos de América para la conexión de una base en alta mar con tierra. El TCraft puede operar en múltiple modos, usándose para conectar una base intermedia de soporte con la principal, y luego ser utilizado como conexión de alta velocidad con la costa, transportando vehículos y otros equipos y carga pesada hasta la playa. El T-Craft se ha concebido como un buque de efecto superficie (Surface-Effect Ship, SES) con un colchón de aire activo entre dos cascos rígidos, que permite al buque operar a desplazamiento completo, con soporte parcial del colchón o con el colchón plenamente activo. Este artículo presenta el trabajo reciente de los autores en el desarrollo de un modelo de evaluación del comportamiento en la mar y la maniobrabilidad del T-Craft, incluyendo la interacción del faldón. Este modelo está basado en la resolución de las ecuaciones de flujo potencial en el dominio del tiempo, mediante el método de los elementos finitos, integrando la superficie libre a lo largo de las líneas de corriente. Este artículo revisa el algoritmo de interacción fluido estructura desarrollado para la simulación de la compleja dinámica de los sellos, en la interfaz entre el colchón de aire y el agua. El algoritmo se basa, por una parte, en un algoritmo explícito alternado que permite la comunicación entre los solvers fluido-dinámico y estructural mediante TCP-IP sockets, intercambiando fuerzas de presión y desplazamientos en los sellos durante la ejecución, y por otra parte, en un innovador esquema de secado-mojado que permite predecir la interacción del agua con los sellos. Se incluyen varios casos de aplicación de la metodología desarrollada sobre el modelo l XR-1B.