La versión p-adaptable del Método de los Elementos Finitos y del Método de los Elementos de Contorno: aplicaciones en elasticidad bidimensional

Este trabajo discute la bondad de las técnicas numéricas p-autoadaptables, comparando el Método de los Elementos Finitos (MEF) y el Método de los Elementos de Contorno(MEC). Se presenta un breve resumen de las herramientas matemáticas necesarias para gobernar el proceso de refinamiento en ambos métodos. Finalmente, se presenta un ejemplo ilustrativo de relevancia práctica en ingeniería, el cual pone de manifiesto la potencia y versatilidad de las técnicas p-adaptables frente a situaciones reales.

Propagación de ondas en sistemas de barras

La acción de una fuerza aplicada bruscamente a un sistema no se transmite simultáneamente a todas las zonas del mismo; en el instante inicial las zonas alejadas de la aplicación de la carga no se alteran, propagándose las deformaciones a través del sólido en forma de ondas elásticas. Si las dimensiones del sistema son grandes, el tiempo que las ondas tardan en atravesarlo es lo suficientemente grande como para ser considerados a efec tos prácticos. El objetivo de este artículo consiste en el estudio de la respuesta de sistemas elásticos formados por barras ante cargas aplicadas brúscamente o fuerzas variables en el tiempo por el método de las ondas características. El método exige unas determinadas relaciones entre las propiedades geométricas y las características del material que permite garantizar la periodicidad del movimiento y seguir, en los casos de ondas de corte y longitudinales, la forma del impulso sin distorsionarlo.

Los Métodos de las Características y Elementos Finitos en propagación de ondas

Se presentan a continuación algunas consideraciones acerca de la propagación de ondas en medios unidimensionales elásticos y viscoplásticos. El objetivo básico es la comparación entre diferentes modelizaciones del sistema establecidas siguiendo las técnicas de Elementos Finitos. Como criterio comparativo se utilizan resultados obtenidos mediante el método de las características, que, para el caso descrito proporciona la solución más aproximada. La justificación del trabajo radica en que las técnicas de Elementos Finitos son fácilmente generalizables a problemas bi o tridimensionales,mientras que las de características pierden rápidamente la ventaja de su simplicidad.

Método interactivo para la extrapolación de ondas sísmicas. Aplicación a la prospección geofísica (Parte 1)

En los análisis sísmicos por reflexión existe una creciente demanda de estimaciones geológicas y geométricas más precisas y sofisticadas. En último término se pretende extraer estimaciones fiables de las secciones de reflexión, estimaciones que afectan a parámetros tales como la velocidad de propagación, densidad y, por supuesto, contenido de hidrocarburos. Ello exige gran precisión en las técnicas de procesado y representación gráfica, pudiendo decirse que los métodos gráficos interactivos son clave en las tendencias futuras de prospección. Un área de investigación actual es el procesado y representación de datos representados por "horizontes". Un horizonte sísmico se define como la superficie que separa dos capas diferentes de roca. Estas interfases están asociadas a reflexiones que se pueden identificar claramente en las secciones sísmicas. Usando una pantalla gráfica adecuada estos horizontes aparecen coloreados en la sección,pudiendo tornarse puntos selectos para definirlas. Si se utilizan secciones concentradas ("stacked") es preciso determinar la posición correcta de las interfases. Es el proceso llamado "emigración". Tradicionalmente se han aplicado fórmulas de representación o transformaciones integrales para realizar este proceso. Sin embargo es suficiente hacer emigrar las curvas seleccionadas (superficies en el caso 3-D) que representan menos del 1% del total de los datos, con lo que el tiempo de CPU se reduce considerablemente. Si se desprecia la información sobre amplitudes, fases y frecuencias la teoría de rayos es el método más rápido y simple para la conversión en profundidad, poniéndose todo el énfasis en el proceso interactivo y en la representación gráfica. Puesto que en estructuras complejas, como las del Mar del Norte o las del Norte de Alemania, las funciones de velocidad nunca se conocen con precisión, es preciso repetir la emigración mediante diferentes velocidades especificadas por el usuario en forma interactiva. En este artículo se presenta la aplicación del procedimiento para horizontes de separación nula (zero-offset sections) y al final se indica la posibilidad de extensión del estudio a casos con separación constante y al análisis de velocidad.

Método interactivo para la extrapolación de ondas sísmicas. Aplicación a la prospección geofísica (Parte 2)

Este artículo es continuación del publicado en el número 2 de la revista "Anales de Ingeniería Mecánica"(diciembre 1982), sobre los métodos de análisis sísmico por reflexión. Como allí se indicó, con estos métodos se pretende determinar la geometría de los horizontes sísmicos en una determinada sección así como las características físicas de las interfases que separa. En el presente artículo se introducen mejoras en lo que se refiere a tiempo compucional y tipos de campos de velocidades. En este nuevo método propuesto para la resolución del proceso de emigración, mediante un método gráfico interactivo, pueden señalarse como ventjas fundamentales, las siguientes: a) Tiempo de resolución acorde con la complejidad del campo de velocidades, produciéndose economía de dicho tiempo, sobre todo en los casos sencillos; b) El número de elementos de la malla y las dimensiones de los mismos son función de la complejidad del campo de velocidades; c) El estudio en cada elemento se realiza de forma general e independiente. General, en cuanto a que los resultados que se obtienen en cada elemento son valores globales, referidos a la malla y al problema en conjunto, e independiente, puesto que los datos necesarios en cada elemento se obtienen a partir de valores nodales; d) Posibilidad de resolución del problema para un campo de velocidades cualesquiera, y por lo tanto, incluyendo los casos de discontinuidades de velocidad, que se evalúan mediante la ley de Snell; e) Posibilidad de extender el método al caso de tres dimensiones, sin más que elegir como elementos de la malla, tetraedros en vez de triángulos. En conjunto, resulta un método flexible, de fácil utilización e interpretación de resultados.

Utilización de la técnica de los elementos discretos en tomografía sísmica

La tomografia es un sistema de obtención de imágenes internas de sólidos a los que no se tiene acceso directo. Ésta técnica ha sido ya de enorme fecundidad en el campo del diagnóstico médico. Hoy día el sistema sigue perfeccionandose y el presente trabajo se enmarca dentro de esta linea, con aplicaciones en el campo de la Ingeniería Civil, utilizando técnicas de discretización del medio investigado, así como la teoría de rayos. Los resultados obtenidos son esperanzadores, y es posible que en un corto espacio de tiempo esta técnica goce de la fiabilidad necesaria en el campo de la Ingeniería.

Modelo de fisuración. Aplicación a estructuras axisimétricas de hormigón armado

El comportamiento mecánico del hormigón sigue siendo un problema no resuelto en cuanto a la modelización según las teorias de la mecánica del continuo, a pesar de tratarse de un material ampliamente extendido en su uso. Una de las facetas más caracteristicas de su comportamiento, la fisuración, representa tal vez el aspecto más dificil de tratar en los modelos bi y tridimensionales. En este articulo se describe un modelo propuesto de comportamiento que salva el compromiso entre complejidad y aplicabilidad: La fisuración es considerada como un deterioro del material -daño, tal vez- y su orientación no es registrada en la historia; al igual que ocurre en la Teoria de la Plasticidad, existe un trabajo de fisuración que sirve para medir el deterioro existente. El modelo es además compatible con la plasticidad como tal, tipica del comportamiento a compresión del hormigón. Dicho modelo, junto con una versión modificada del método de longitud de arco para resolución de problemas no lineales, es implementado y utilizado para la resolución de algunos ejemplos.

Localización de reflectores mediante la utilización de la teoría de rayos a nivel dinámico

A partir de los años 60, comienza un desarrollo importante de los métodos que utilizan las características de reflexión del medio para la determinación de la estructura del subsuelo. Se han desarrollado técnicas muy precisas en el campo de la detección, fundamentalmente de recursos energéticos, pero aplicadas únicamente a nivel cinemática. Se intenta en este trabajo mostrar la posibilidad de utilizar la teoría de rayos a un segundo nivel, dinámico, para extraer más información del medio en menos tiempo, y, por lo tanto, con menor costo. Se ha utilizado el método de los elementos discretos triangulares = Since the sixties an important developmen of the methods which use the characteristics of environment reflexion for the determination of subsoil structure has taken place. Some very accurate technics have been developed in the detection field, principally energetic resources technics, but they have been only applied to a kinematic level. In this work we try to show the possibility of using ray theory in a second level: dinamic level, to get more information about the environment in less time, and so with less cost. The method of triangular discret elements has been used.

Aplicación de la teoría de fiabilidad al establecimiento de niveles de seguridad en túneles

La perforación de túneles constituye una tarea clave en la mejora de las comunicaciones y en el desarrollo de la industria minera. Debido a su complejidad y al alto número de factores que influyen en su desdarrollo, este tipo de proyectos han sido considerados tradicionalmente obras de alto e indefinido riesgo. La difusión del llamado Nuevo Método Austriaco, ha supuesto un gran impulso en la tecnificación de este tipo de proyectos. Paradójicamente, aún en la actualidad, el problema de los márgenes de seguridad en el cálculo se suele reducir en la práctica a la evaluación de un único coeficiente cuya definición no satisface los criterios actualmente aceptados en este campo. En este trabajo se ha pretendido, de acuerdo a las tendencias actuales, introducir las técnicas de Fiabilidad Estructural en el modelo de cálculo del sostenimiento propuesto por el NATM. En el desarrollo se ha tomado como punto de partida el método de cálculo de las Curvas Características, para lo cual se han formulado las correspondientes leyes de comportamiento del terreno y del sostenimiento. Se han establecido 3 Estados Límite de referencia y se ha seleccionado un conjunto de variables básicas representativas del problema estudiado. Finalmente la aplicación de los conceptos de Fiabilidad Estructural en Nivel II junto con el algoritmo de estimación puntual de E. Rosenblueth, han permitido el desarrollo de modelos de cálculo que permiten la determinación de los Indices de Fiabilidad Estructural respecto a los 3 Estados Límite considerados. El esquema descrito ha sido implementado en programas de cálculo. Estos programas de cálculo han sido utilizados en el estudio de diversos ejemplos correspondientes a situaciones reales de proyecto y se han desarrollado diversos estudios paramétricos. Todo ello ha permitido comprobar la validez del procedimiento desarrollado y llegar al establecimiento de un propuesta de calibración de factores parciales de ponderación que podrían ser incorporados en recomendaciones prácticas de cálculo.

Optimización del acondicionamiento sísmico en puentes mediante el aislamiento en la base del tablero

Una solución para conseguir minimizar los daños que los terremotos causan en los puentes consiste en el aislamiento sísmico en la base, utilizando apoyos deslizantes entre pilas y tablero que transmitan únicamente las cargas verticales, mientras que los esfuerzos horizontales son transmitidos por unos dispositivos especiales que proporcionan flexibilidad horizontal, fuerzas de reacción y disipación de energía. En este trabajo se desarrolla un método para definir cuales serían las características ideales que deberían tener esos dispositivos de aislamiento para conseguir que las fuerzas horizontales transmitidas a las pilas sean inferiores a su límite elástico, minimizando los desplazamientos permitidos del tablero, y consiguiendo idénticas deformaciones en todos los puntos de la estructura en cada momento. Este método se desarrollará para todo tipo de puentes y se particularizará para tipos propuestos en los estudios realizados en la elaboración del Eurocódigo Sísmico analizándose las variables que intervienen, y las respuestas dinámicas esperadas.

Influencia de la carga aplicada sobre la microdureza en polímeros

En este trabajo se presenta un modelo mecánico que permite relacionar la microdureza de polímeros con la carga aplicada. Los resultados de las medidas experimentales realizadas sobre polietilenos de diferentes tipos se han ajustado al modelo y se discuten los coeficientes resultantes del ajuste. A partir de los valores de estos coeficientes se propone un criterio para la selección de la carga a aplicar.

Caracterización de vehículos ferroviarios para analizar la sensibilidad de sus parámetros dinámicos

Este trabajo forma parte de un proyecto más amplio, cuya finalidad es realizar un estudio de sensibilidad que permita determinar qué grado de exactitud se requiere en la definición de los principales parámetros que influyen en la dinámica de los vehículos ferroviarios. Como paso previo a este análisis de sensibilidad, es preciso seleccionar un vehículo de referencia en el que basar dicho estudio. Dicho vehículo se empleará como punto de partida, para después variar el valor de los parámetros que lo caracterizan, a fin de determinar su influencia sobre la dinámica de marcha. Para alcanzar este objetivo se ha construido una base de datos de vehículos ferroviarios, que almacena las características de los parámetros que afectan al comportamiento dinámico de los vehículos. El valor de referencia correspondiente a cada parámetro puede asociarse al promedio de los valores almacenados en la base de datos, y sus rangos de variación pueden obtenerse a partir de las medidas de dispersión de estos datos. Se ha realizado un test de normalidad sobre los datos almacenados en la base de datos, encontrándose que la hipótesis de normalidad resulta, en general, inapropiada, lo que requiere el uso de estimadores más robustos de los habituales. Aunque el estudio de sensibilidad podría basarse en el vehículo ficticio definido por los valores promedio de cada parámetro, cabe la posibilidad de que este vehículo promedio presente un comportamiento dinámico distinto al que tendría un vehículo real. Para evitar este posible contratiempo, se han utilizado técnicas de escalado multidimensional para seleccionar los registros de la base de datos más próximos al vehículo promedio. Para determinar cuál de estos vehículos resulta más apropiado, se ha comparado su comportamiento dinámico empleando técnicas de simulación de sistemas multicuerpo. A partir de los resultados obtenidos en estas simulaciones, se ha valorado el comportamiento de cada modelo desde el punto de vista de la seguridad, de la agresión a la vía y del confort, siguiendo las indicaciones de la norma EN-14363. Finalmente, se ha seleccionado como vehículo de referencia aquél que, presentando un comportamiento adecuado, requiere el menor esfuerzo computacional para llevar a cabo las simulaciones.

Verificación experimental de las correlaciones de transferencia de calor por ebullición en película, en piscina, en torno a esferas

La ebullición en película es el mecanismo de transferencia de calor básico que acopla térmicamente un líquido y una superficie caliente cuando existe una gran diferencia de temperatura entre ambos. El conocimiento preciso del comportamiento de este mecanismo térmico en torno a esferas es necesario para el análisis de seguridad de escenarios industriales en los que exista contacto entre un líquido y un material fundido fragmentado, generalmente en piezas esferoidales, y resulta esencial para garantizar la seguridad de los reactores nucleares ante escenarios accidentales de muy baja probabilidad, pero de gran severidad, en los que se postule la rotura en guillotina del circuito de refrigeración del reactor simultánea al fallo total activo del sistema de refrigeración de emergencia del núcleo. En tal hipotético escenario, se produciría la fusión del combustible en el plazo de algunas horas, con dispersión y relocalización de fragmentos sólidos esferoidales a muy alta temperatura. Para hacer frente a este escenario se precisa desarrollar procedimientos automáticos y manuales de operación de emergencia, resultando imprescindible disponer de modelos térmicos confiables, con un calificado soporte experimental, que permitan analizar de forma realista la refrigeración por ebullición en película, en modo ebullición en piscina, de las geometrías esferoidales sólidas resultantes. En el presente trabajo UNET-UPM abordan la verificación de las correlaciones más conocidas para ebullición en película en piscina en torno a esferas, mediante la comparación de las mismas con los resultados experimentales obtenidos por Liu-Theofanous. Algunos de los aspectos resaltantes son la limitada aplicabilidad de la correlación empírica de Frederking-Clark, el buen ajuste que muestran las correlaciones de Tou- Tso y de Grigoriew frente a los datos experimentales usados, y las deficiencias que muestran las correlaciones que toman en cuenta el subenfriamiento del líquido, sobre todo a altos niveles de subenfriamiento. En una futura segunda fase de esta investigación, de tipo analítica-numérica-experimental, se abordará en el desarrollo de nuevas correlaciones semi-empíricas, de mejor ajuste, que permitan una mejor capacidad predictiva en los modelos.

Revisión general de la transferencia de calor por ebullición en película en piscina

La ebullición en película es el mecanismo de transferencia de calor básico que acopla térmicamente un líquido, saturado o subenfriado, y una superficie caliente cuando existe una gran diferencia de temperatura entre ambos. Dicho mecanismo presenta la complejidad física asociada a la existencia de una capa límite convectiva de vapor en torno a la superficie caliente, así como al acoplamiento térmico radiante entre la superficie y los dos medios participativos que la rodean (vapor y líquido). Este mecanismo tiene aplicaciones tecnológicas de interés en propulsión aeroespacial, en sistemas criogénicos, en procesos industriales metalúrgicos, y en otras áreas de la Ingeniería en donde existe la necesidad de enfriar mediante ebullición superficies a alta temperatura. Por otra parte, la refrigeración por ebullición en película es un mecanismo que resulta de importancia esencial en el análisis de la seguridad de procesos industriales en los que sea previsible la interacción de agua y materiales fundidos, entre los que cabe destacar el caso de los reactores nucleares ante escenarios de riesgo en los que el acoplamiento térmico convectivo-radiante entre el combustible nuclear y el refrigerante, e incluso entre la superficie exterior de la vasija y el refrigerante, garantiza el cumplimiento de los criterios de aceptación en caso de fallo mecánico importante del sistema de refrigeración del reactor. En el presente trabajo, y para identificar el marco en que UNET-UPM desarrollará sus investigaciones sobre ebullición en película en piscina en torno a esferas, se realiza una revisión de los estudios más relevantes desarrollados en esta temática, en la que las primeras investigaciones teóricas y experimentales, de tipo fundamental, se remontan al periodo 1950-65 en que se desarrollaron los primeros proyectos de navegación espacial, y se diseñaron y construyeron los primeros reactores nucleares de aplicación comercial, promoviéndose posteriormente numerosas investigaciones que aportaron correlaciones "ad-hoc" de aplicación restringida a ciertas tecnologías, así como otras de propósito general, de menor precisión pero con mayor soporte teórico.

Hydrodynamic forces on heave plates for offshore systems oscillating close to the seabed or the free surface

La energía eólica marina es uno de los recursos energéticos con mayor proyección pudiendo contribuir a reducir el consumo de combustibles fósiles y a cubrir la demanda de energía en todo el mundo. El concepto de aerogenerador marino está basado en estructuras fijas como jackets o en plataformas flotantes, ya sea una semisumergible o una TLP. Se espera que la energía eólica offshore juegue un papel importante en el perfil de producción energética de los próximos años; por tanto, las turbinas eólicas deben hacerse más fables y rentables para ser competitivas frente a otras fuentes de energía. Las estructuras flotantes pueden experimentar movimientos resonantes en estados de la mar con largos períodos de oleaje. Estos movimientos disminuyen su operatividad y pueden causar daños en los componentes eléctricos de las turbinas y en las palas, también en los risers y moorings. La respuesta de la componente vertical del movimiento puede reducirse mediante diferentes actuaciones: (1) aumentando la amortiguación del sistema, (2) manteniendo el período del movimiento vertical fuera del rango de la energía de la ola, y (3) reduciendo las fuerzas de excitación verticales. Un ejemplo típico para llevar a cabo esta reducción son las "Heave Plates". Las heave plates son placas que se utilizan en la industria offshore debido a sus características hidrodinámicas, ya que aumentan la masa añadida y la amortiguación del sistema. En un análisis hidrodinámico convencional, se considera una estructura sometida a un oleaje con determinadas características y se evalúan las cargas lineales usando la teoría potencial. El amortiguamiento viscoso, que juega un papel crucial en la respuesta en resonancia del sistema, es un dato de entrada para el análisis. La tesis se centra principalmente en la predicción del amortiguamiento viscoso y de la masa añadida de las heave plates usadas en las turbinas eólicas flotantes. En los cálculos, las fuerzas hidrodinámicas se han obtenido con el f n de estudiar cómo los coeficientes hidrodinámicos de masa añadida5 y amortiguamiento varían con el número de KC, que caracteriza la amplitud del movimiento respecto al diámetro del disco. Por otra parte, se ha investigado la influencia de la distancia media de la ‘heave plate’ a la superficie libre o al fondo del mar, sobre los coeficientes hidrodinámicos. En este proceso, un nuevo modelo que describe el trabajo realizado por la amortiguación en función de la enstrofía, es descrito en el presente documento. Este nuevo enfoque es capaz de proporcionar una correlación directa entre el desprendimiento local de vorticidad y la fuerza de amortiguación global. El análisis también incluye el estudio de los efectos de la geometría de la heave plate, y examina la sensibilidad de los coeficientes hidrodinámicos al incluir porosidad en ésta. Un diseño novedoso de una heave plate, basado en la teoría fractal, también fue analizado experimentalmente y comparado con datos experimentales obtenidos por otros autores. Para la resolución de las ecuaciones de Navier Stokes se ha usado un solver basado en el método de volúmenes finitos. El solver usa las librerías de OpenFOAM (Open source Field Operation And Manipulation), para resolver un problema multifásico e incompresible, usando la técnica VOF (volume of fluid) que permite capturar el movimiento de la superficie libre. Los resultados numéricos han sido comparados con resultados experimentales llevados a cabo en el Canal del Ensayos Hidrodinámicos (CEHINAV) de la Universidad Politécnica de Madrid y en el Canal de Experiencias Hidrodinámicas (CEHIPAR) en Madrid, al igual que con otros experimentos realizados en la Escuela de Ingeniería Mecánica de la Universidad de Western Australia. Los principales resultados se presentan a continuación: 1. Para pequeños valores de KC, los coeficientes hidrodinámicos de masa añadida y amortiguamiento incrementan su valor a medida que el disco se aproxima al fondo marino. Para los casos cuando el disco oscila cerca de la superficie libre, la dependencia de los coeficientes hidrodinámicos es más fuerte por la influencia del movimiento de la superficie libre. 2. Los casos analizados muestran la existencia de un valor crítico de KC, donde la tendencia de los coeficientes hidrodinámicos se ve alterada. Dicho valor crítico depende de la distancia al fondo marino o a la superficie libre. 3. El comportamiento físico del flujo, para valores de KC cercanos a su valor crítico ha sido estudiado mediante el análisis del campo de vorticidad. 4. Introducir porosidad al disco, reduce la masa añadida para los valores de KC estudiados, pero se ha encontrado que la porosidad incrementa el valor del coeficiente de amortiguamiento cuando se incrementa la amplitud del movimiento, logrando un máximo de damping para un disco con 10% de porosidad. 5. Los resultados numéricos y experimentales para los discos con faldón, muestran que usar este tipo de geometrías incrementa la masa añadida cuando se compara con el disco sólido, pero reduce considerablemente el coeficiente de amortiguamiento. 6. Un diseño novedoso de heave plate basado en la teoría fractal ha sido experimentalmente estudiado a diferentes calados y comparado con datos experimentales obtenidos por otro autores. Los resultados muestran un comportamiento incierto de los coeficientes y por tanto este diseño debería ser estudiado más a fondo. ABSTRACT Offshore wind energy is one of the promising resources which can reduce the fossil fuel energy consumption and cover worldwide energy demands. Offshore wind turbine concepts are based on either a fixed structure as a jacket or a floating offshore platform like a semisubmersible, spar or tension leg platform. Floating offshore wind turbines have the potential to be an important part of the energy production profile in the coming years. In order to accomplish this wind integration, these wind turbines need to be made more reliable and cost efficient to be competitive with other sources of energy. Floating offshore artifacts, such oil rings and wind turbines, may experience resonant heave motions in sea states with long peak periods. These heave resonances may increase the system downtime and cause damage on the system components and as well as on risers and mooring systems. The heave resonant response may be reduced by different means: (1) increasing the damping of the system, (2) keeping the natural heave period outside the range of the wave energy, and (3) reducing the heave excitation forces. A typical example to accomplish this reduction are “Heave Plates”. Heave plates are used in the offshore industry due to their hydrodynamic characteristics, i.e., increased added mass and damping. Conventional offshore hydrodynamic analysis considers a structure in waves, and evaluates the linear and nonlinear loads using potential theory. Viscous damping, which is expected to play a crucial role in the resonant response, is an empirical input to the analysis, and is not explicitly calculated. The present research has been mainly focused on the prediction of viscous damping and added mass of floating offshore wind turbine heave plates. In the calculations, the hydrodynamic forces have been measured in order to compute how the hydrodynamic coefficients of added mass1 and damping vary with the KC number, which characterises the amplitude of heave motion relative to the diameter of the disc. In addition, the influence on the hydrodynamic coefficients when the heave plate is oscillating close to the free surface or the seabed has been investigated. In this process, a new model describing the work done by damping in terms of the flow enstrophy, is described herein. This new approach is able to provide a direct correlation between the local vortex shedding processes and the global damping force. The analysis also includes the study of different edges geometry, and examines the sensitivity of the damping and added mass coefficients to the porosity of the plate. A novel porous heave plate based on fractal theory has also been proposed, tested experimentally and compared with experimental data obtained by other authors for plates with similar porosity. A numerical solver of Navier Stokes equations, based on the finite volume technique has been applied. It uses the open-source libraries of OpenFOAM (Open source Field Operation And Manipulation), to solve 2 incompressible, isothermal immiscible fluids using a VOF (volume of fluid) phase-fraction based interface capturing approach, with optional mesh motion and mesh topology changes including adaptive re-meshing. Numerical results have been compared with experiments conducted at Technical University of Madrid (CEHINAV) and CEHIPAR model basins in Madrid and with others performed at School of Mechanical Engineering in The University of Western Australia. A brief summary of main results are presented below: 1. At low KC numbers, a systematic increase in added mass and damping, corresponding to an increase in the seabed proximity, is observed. Specifically, for the cases when the heave plate is oscillating closer to the free surface, the dependence of the hydrodynamic coefficients is strongly influenced by the free surface. 2. As seen in experiments, a critical KC, where the linear trend of the hydrodynamic coefficients with KC is disrupted and that depends on the seabed or free surface distance, has been found. 3. The physical behavior of the flow around the critical KC has been explained through an analysis of the flow vorticity field. 4. The porosity of the heave plates reduces the added mass for the studied porosity at all KC numbers, but the porous heave plates are found to increase the damping coefficient with increasing amplitude of oscillation, achieving a maximum damping coefficient for the heave plate with 10% porosity in the entire KC range. 5. Another concept taken into account in this work has been the heave plates with flaps. Numerical and experimental results show that using discs with flaps will increase added mass when compared to the plain plate but may also significantly reduce damping. 6. A novel heave plate design based on fractal theory has tested experimentally for different submergences and compared with experimental data obtained by other authors for porous plates. Results show an unclear behavior in the coefficients and should be studied further. Future work is necessary in order to address a series of open questions focusing on 3D effects, optimization of the heave plates shapes, etc.