Universidad Politécnica de Cartagena: Escuela Técnica Superior de Ingeniería Naval y Oceánica
Universidad Politécnica de Cartagena: Escuela Técnica Superior de Ingeniería Naval y Oceánica
Universidad Politécnica de Cartagena: Escuela Técnica Superior de Ingeniería Naval y Oceánica
Muchísimas gracias
ÍNDICE GENERAL
LISTADO DE SÍMBOLOS............................................................................................9
PRÓLOGO.....................................................................................................................21
5
5.1 Introducción ..................................................................................................... 68
5.2 Metodología y aplicación................................................................................. 73
5.3 Resultados ........................................................................................................ 75
5.3.1 Curvas hidrostáticas para los distintos trimados. ..................................... 82
6
1.5.2.2 Ecuaciones de Navier-Stokes.............................................................. 145
1.5.2.2.1 La suposición de incompresibilidad .............................................. 146
1.5.3 Capa límite.............................................................................................. 146
1.5.3.1 Desprendimiento de la capa límite...................................................... 148
1.6 Modelos físico-matemáticos empleados en CFD’s........................................ 149
1.6.1 Modelos de flujo potencial ..................................................................... 149
1.6.2 Modelos para flujo ideal ......................................................................... 150
1.6.3 Modelos para flujos ideales, estacionarios y rotacionales ...................... 150
1.6.4 Modelos parabólicos de las ecuaciones de Navier-Stokes ..................... 150
1.6.5 Solución directa de las ecuaciones de Navier -Stokes............................ 151
1.6.6 Modelos de flujo incompresible ............................................................. 151
1.6.7 Modelos para la resolución de las turbulencias ...................................... 152
1.6.7.1 Simulación directa de las ecuaciones (Direct Simulation, DS) .......... 152
1.6.7.2 Simulación de grandes vórtices (Large Eddy Simulation, LES) ........ 152
1.6.7.3 Modelos URANS (Unsteady Reynolds Averaged Navier-Stokes)..... 153
1.6.7.4 Modelos que promedian las ecuaciones de Navier-Stokes (RANS) ... 153
1.7 Modelos de turbulencia .................................................................................. 154
1.7.1 Modelos algebraicos (de cero ecuaciones de cierre) .............................. 155
1.7.2 Modelos de ecuación de energía cinética turbulenta .............................. 156
1.7.2.1 La ecuación de energía cinética turbulenta ......................................... 156
1.7.3 Modelos de una ecuación de cierre......................................................... 157
1.7.4 Modelos de dos ecuaciones de cierre ..................................................... 158
2. METODOLOGÍA DEL ENSAYO DE REMOLQUE DE UNA CARENA ........ 161
2.1 Proceso de cálculo de un código CFD ........................................................... 161
2.1.1 Pre-proceso ............................................................................................. 162
2.1.1.1 Creación o importación de la geometría ............................................. 162
2.1.1.1.1 Creación del volumen de control .................................................. 163
2.1.1.2 Asignación de condiciones de contorno ............................................. 163
2.1.1.3 Materiales............................................................................................ 172
2.1.1.4 Datos del problema ............................................................................. 173
2.1.1.4.1 Discretización espacial .................................................................. 174
Clasificación de los métodos de discretización espacial ............................ 174
Tipos de mallado ........................................................................................ 177
7
Generación de la malla ............................................................................... 179
Tamaño del elemento de malla .................................................................. 180
2.1.1.4.2 Discretización temporal ................................................................ 181
2.1.1.5 Selección del modelo de turbulencia adecuado .................................. 182
2.1.2 Cálculo .................................................................................................... 182
2.1.3 Post-proceso ........................................................................................... 183
2.2 Resumen de parámetros ................................................................................. 183
3. APLICACIÓN EN UN CASO PRÁCTICO ......................................................... 184
3.1 Definición de la geometría ............................................................................. 184
3.1.1 Generación del volumen de control ........................................................ 185
3.2 Asignación de las condiciones de contorno ................................................... 186
3.2.1 Campo de presiones y de velocidades .................................................... 186
3.2.2 Condiciones de la capa límite ................................................................. 187
3.2.3 Mallado ................................................................................................... 189
3.2.4 Modelo de turbulencia ............................................................................ 191
3.3 Visión, análisis y comparación de resultados ................................................ 199
4. CONCLUSIONES ................................................................................................ 224
8
LISTADO DE SÍMBOLOS
9
Tamaño del elemento de malla
Puntal a la cubierta principal
Puntal a la cubierta superior
Semiángulo de entrada en la flotación
Energía cinética turbulenta; Conductividad térmica; Velocidad de escala
Rugosidad del casco
Altura del centro de carena sobre la base
Altura metacéntrica longitudinal
Altura metacéntrica transversal
Longitud de escala turbulenta
Eslora del bulbo
Longitud de mezcla
Eslora
Longitud de la cuchara
Eslora de desplazamiento
Eslora en la flotación
Longitud de onda
Eslora entre perpendiculares
Posición longitudinal del centro de carena
Posición longitudinal del centro de flotación
Masa
Momento para alterar el trimado un centímetro
Vector normal a una superficie
Número mínimo de intervalos temporales
Número final de intervalos temporales
Presión
Presión absoluta
Perpendicular de popa
Perpendicular de proa
Deposición de calor por unidad de masa
Factor de crecimiento asignado
Resistencia debida al coeficiente de correlación modelo-buque
Resistencia de los apéndices
Resistencia de presión producida por el bulbo
Resistencia de fricción
10
Constante del gas
Resistencia total
Resistencia debida a la inmersión del espejo
Resistencia viscosa
Resistencia por formación de olas
Número de Reynolds
Superficie mojada del buque
Superficie mojada de los apéndices
Superficie mojada del bulbo
Superficie mojada del buque con apéndices
Superficie mojada de la tobera
Superficie mojada del timón y del codaste cerrado
Tiempo
Calado; Temperatura
Calado medio
Calado en la perpendicular de popa
Calado en la perpendicular de proa
Tiempo de simulación
Toneladas por centímetro de inmersión
Nivel de intensidad de turbulencia
Componente de la velocidad; Componente fluctuante de la velocidad
Valor medio de la velocidad
Velocidad de fricción en la pared
Velocidad promediada en el tiempo
Velocidad promediada en el tiempo adimensional
Componente de la velocidad; Volumen específico
′ Velocidad local
Velocidad media del fluido
Velocidad el buque; Volumen
Componente de la velocidad
Longitud del flujo laminar
Distancia normal a la pared
Distancia adimensional a la pared
∆ Desplazamiento
Volumen de carena o de desplazamiento
11
Operador de Laplace
Espesor de la capa límite
Incremento de tiempo
Incremento de tiempo final
Disipación de energía cinética turbulenta
Longitud de onda
Esfuerzo cortante; Tensiones de Reynolds
Esfuerzo cortante normal
Esfuerzo cortante tangencial
Esfuerzo cortante en la pared
Viscosidad dinámica
Viscosidad dinámica turbulenta
Viscosidad cinemática
Rango específico de disipación de energía
Densidad
⨂ Sección media
12
ÍNDICE DE TABLAS
13
Tabla 4.8. Valores de la resistencia a carga parcial obtenidos mediante el método de
Sander M. Calisal y Dan McGreer. ....................................................................... 125
Tabla 4.9. Comparación de los valores de las resistencias obtenidos mediante los
distintos métodos de predicción estadística. ......................................................... 127
Tabla 5.1. Valores de TIL en función del número de Froude. .................................... 164
Tabla 5.2. Coordenadas de los puntos que forman el volumen de control. ................. 185
Tabla 5.3. Valores de la longitud del flujo laminar ‘x’ para las distintas velocidades. 188
Tabla 5.4. Valores de los parámetros de turbulencia iniciales para 20 nudos. ............ 188
Tabla 5.5. Características de los distintos mallados empleados en el análisis de
sensibilidad. ........................................................................................................... 189
Tabla 5.6. Valores de los parámetros de turbulencia para las distintas velocidades. .. 197
Tabla 5.7. Valores iniciales de los parámetros referentes a la discretización espacial
para las distintas velocidades. ............................................................................... 198
Tabla 5.8. Valores de los parámetros referentes a la discretización temporal para las
distintas velocidades.............................................................................................. 198
Tabla 5.9. Valores de las distancias de los cortes para la obtención del perfil de olas
respecto a crujía..................................................................................................... 198
Tabla 5.10. Valores de la resistencia para los distintos modelos de turbulencia
seleccionados......................................................................................................... 199
Tabla 5.11. Valores de la resistencia viscosa y por formación de olas respecto a la
resistencia total. ..................................................................................................... 201
Tabla 6.1. Comparación de los valores de las resistencias obtenidos mediante los
distintos métodos de predicción de resistencia. .................................................... 229
Tabla 6.2. Comparación de métodos de predicción de resistencia respecto al de Calisal
(plena carga) para distintas velocidades. ............................................................... 233
14
ÍNDICE DE FIGURAS
15
Figura 5.22. Ejemplo de mallado multibloque. ........................................................... 179
Figura 5.23. Ejemplo de convergencia de las fuerzas de presión para la velocidad de 20
nudos. .................................................................................................................... 183
Figura 5.24. Parámetros iniciales a introducir en Tdyn. ............................................. 183
Figura 5.25. Casco del buque dividido a una altura de 4 metros................................. 184
Figura 5.26. Volumen de control................................................................................. 185
Figura 5.27. Identificación de las entidades geométricas en capas. ............................ 186
Figura 5.28. Campo de velocidades. ........................................................................... 187
Figura 5.29. Campo de presiones. ............................................................................... 187
Figura 5.30. Comparación entre los diferentes mallados. ........................................... 190
Figura 5.31. Distribución de presiones en el casco a 3 nudos. .................................... 203
Figura 5.32. Distribución de presiones en la superficie libre a 3 nudos...................... 203
Figura 5.33. Distribución de presiones en el volumen de control a 3 nudos............... 203
Figura 5.34. Distribución de velocidades en el casco a 3 nudos. ................................ 204
Figura 5.35. Distribución de velocidades en la superficie libre a 3 nudos. ................. 204
Figura 5.36. Distribución de velocidades en el volumen de control a 3 nudos. .......... 204
Figura 5.37. Mapa de olas en la superficie libre a 3 nudos. ........................................ 205
Figura 5.38. Mapa de olas en la superficie libre a 3 nudos. ........................................ 205
Figura 5.39. Mapa de olas en el volumen de control a 3 nudos. ................................. 205
Figura 5.40. Distribución de presiones en el casco a 5 nudos. .................................... 207
Figura 5.41. Distribución de presiones en la superficie libre a 5 nudos...................... 207
Figura 5.42. Distribución de presiones en el volumen de control a 5 nudos............... 207
Figura 5.43. Distribución de velocidades en el casco a 5 nudos. ................................ 208
Figura 5.44. Distribución de velocidades en la superficie libre a 5 nudos. ................. 208
Figura 5.45. Distribución de velocidades en el volumen de control a 5 nudos. .......... 208
Figura 5.46. Mapa de olas en la superficie libre a 5 nudos. ........................................ 209
Figura 5.47. Mapa de olas en la superficie libre a 5 nudos. ........................................ 209
Figura 5.48. Mapa de olas en el volumen de control a 5 nudos. ................................. 209
Figura 5.49. Distribución de presiones en el casco a 10 nudos. .................................. 211
Figura 5.50. Distribución de presiones en la superficie libre a 10 nudos.................... 211
Figura 5.51. Distribución de presiones en el volumen de control a 10 nudos............. 211
Figura 5.52. Distribución de velocidades en el casco a 10 nudos. .............................. 212
Figura 5.53. Distribución de velocidades en la superficie libre a 10 nudos. ............... 212
Figura 5.54. Distribución de velocidades en el volumen de control a 10 nudos. ........ 212
Figura 5.55. Mapa de olas en la superficie libre a 10 nudos. ...................................... 213
Figura 5.56. Mapa de olas en la superficie libre a 10 nudos. ...................................... 213
Figura 5.57. Mapa de olas en el volumen de control a 10 nudos. ............................... 213
Figura 5.58. Distribución de presiones en el casco a 15 nudos. .................................. 215
Figura 5.59. Distribución de presiones en la superficie libre a 15 nudos.................... 215
Figura 5.60. Distribución de presiones en el volumen de control a 15 nudos............. 215
Figura 5.61. Distribución de velocidades en el casco a 15 nudos. .............................. 216
Figura 5.62. Distribución de velocidades en la superficie libre a 15 nudos. ............... 216
Figura 5.63. Distribución de velocidades en el volumen de control a 15 nudos. ........ 216
Figura 5.64. Mapa de olas en la superficie libre a 15 nudos. ...................................... 217
16
Figura 5.65. Mapa de olas en la superficie libre a 15 nudos. ...................................... 217
Figura 5.66. Mapa de olas en el volumen de control a 15 nudos. ............................... 217
Figura 5.67. Distribución de presiones en el casco a 20 nudos. .................................. 219
Figura 5.68. Distribución de presiones en la superficie libre a 20 nudos.................... 219
Figura 5.69. Distribución de presiones en el volumen de control a 20 nudos............. 219
Figura 5.70. Distribución de velocidades en el casco a 20 nudos. .............................. 220
Figura 5.71. Distribución de velocidades en la superficie libre a 20 nudos. ............... 220
Figura 5.72. Distribución de velocidades en el volumen de control a 20 nudos. ........ 220
Figura 5.73. Mapa de olas en la superficie libre a 20 nudos. ...................................... 221
Figura 5.74. Mapa de olas en la superficie libre a 20 nudos. ...................................... 221
Figura 5.75. Mapa de olas en el volumen de control a 20 nudos. ............................... 221
Figura 5.76. Distribución de presiones (izquierda) y velocidades (derecha) en la zona
de proa. .................................................................................................................. 223
17
ÍNDICE DE GRÁFICAS
18
Gráfica 3.27. Curvas de la posición longitudinal del centro de carena y del centro de la
flotación para un trimado de +0,5 m. ...................................................................... 95
Gráfica 3.28. Curva de la altura del centro de carena para un trimado de +0,5 m. ....... 96
Gráfica 3.29. Curva del radio metacéntrico transversal para un trimado de +0,5 m. .... 96
Gráfica 3.30. Curva del radio metacéntrico longitudinal para un trimado de +0,5 m. .. 97
Gráfica 3.31. Curva de la distancia del metacentro transversal sobre la base para un
trimado de +0,5 m. .................................................................................................. 97
Gráfica 3.32. Curva de la distancia del metacentro longitudinal sobre la base para un
trimado de +0,5 m. .................................................................................................. 98
Gráfica 3.33. Curva de las toneladas por centímetro de inmersión para un trimado de
+0,5 m. .................................................................................................................... 98
Gráfica 3.34. Curva del momento para alterar el trimado un centímetro para un trimado
de +0,5 m. ................................................................................................................ 99
Gráfica 3.35. Curva del desplazamiento para un trimado de +1 m. ............................ 100
Gráfica 3.36. Curvas de la superficie mojada y del área de la flotación para un trimado
de +1 m. ................................................................................................................. 100
Gráfica 3.37. Curvas de los coeficientes de formas para un trimado de +1 m. ........... 101
Gráfica 3.38. Curvas de la posición longitudinal del centro de carena y del centro de la
flotación para un trimado de +1 m. ....................................................................... 101
Gráfica 3.39. Curva de la altura del centro de carena para un trimado de +1 m. ........ 102
Gráfica 3.40. Curva del radio metacéntrico transversal para un trimado de +1 m. ..... 102
Gráfica 3.41. Curva del radio metacéntrico longitudinal para un trimado de +1 m. ... 103
Gráfica 3.42. Curva de la distancia del metacentro transversal sobre la base para un
trimado de +1 m. ................................................................................................... 103
Gráfica 3.43. Curva de la distancia del metacentro longitudinal sobre la base para un
trimado de +1 m. ................................................................................................... 104
Gráfica 3.44. Curva de las toneladas por centímetro de inmersión para un trimado de
+1 m. ..................................................................................................................... 104
Gráfica 3.45. Curva del momento para alterar el trimado un centímetro para un trimado
de +1 m. ................................................................................................................. 105
Gráfica 3.46. Curva del desplazamiento para un trimado de +1,5 m. ......................... 106
Gráfica 3.47. Curvas de la superficie mojada y del área de la flotación para un trimado
de +1,5 m. .............................................................................................................. 106
Gráfica 3.48. Curvas de los coeficientes de formas para un trimado de +1,5 m. ........ 107
Gráfica 3.49. Curvas de la posición longitudinal del centro de carena y del centro de la
flotación para un trimado de +1,5 m. .................................................................... 107
Gráfica 3.50. Curva de la altura del centro de carena para un trimado de +1,5 m. ..... 108
Gráfica 3.51. Curva del radio metacéntrico transversal para un trimado de +1,5 m. .. 108
Gráfica 3.52. Curva del radio metacéntrico longitudinal para un trimado de +1,5 m. 109
Gráfica 3.53. Curva de la distancia del metacentro transversal sobre la base para un
trimado de +1,5 m. ................................................................................................ 109
Gráfica 3.54. Curva de la distancia del metacentro longitudinal sobre la base para un
trimado de +1,5 m. ................................................................................................ 110
19
Gráfica 3.55. Curva de las toneladas por centímetro de inmersión para un trimado de
+1,5 m. .................................................................................................................. 110
Gráfica 3.56. Curva del momento para alterar el trimado un centímetro para un trimado
de +1,5 m. .............................................................................................................. 111
Gráfica 4.1. Curvas de las componentes de la resistencia obtenidas mediante el método
de J. Holtrop y G. G. J. Mennen.............................................................................115
Gráfica 4.2. Curva de resistencia obtenida mediante el método de J. Holtrop y G. G. J.
Mennen. ................................................................................................................. 116
Gráfica 4.3. Curvas de la resistencia viscosa y por formación de olas respecto a la
resistencia total. ..................................................................................................... 118
Gráfica 4.4. Curvas de resistencia del casco con bulbo y sin bulbo obtenida mediante el
método de Amadeo García Gómez. ...................................................................... 121
Gráfica 4.5. Curva de resistencia obtenida mediante el método de Van Oortmersen. 123
Gráfica 4.6. Curvas de resistencia a plena carga y a carga parcial obtenidas mediante el
método de Sander M. Calisal y Dan McGreer. ..................................................... 126
Gráfica 4.7. Comparación de las curvas de resistencia obtenidas mediante los distintos
métodos de predicción estadística. ........................................................................ 128
Gráfica 5.1. Comparación de los distintos modelos de turbulencia empleados en el
análisis de sensibilidad. ......................................................................................... 191
Gráfica 5.2. Curvas de resistencia del modelo de turbulencia K-Epsilon High Reynolds.
............................................................................................................................... 192
Gráfica 5.3. Curvas de resistencia del modelo de turbulencia K-Epsilon Low Reynolds.
............................................................................................................................... 193
Gráfica 5.4. Curvas de resistencia del modelo de turbulencia K-Epsilon Two Layers.
............................................................................................................................... 193
Gráfica 5.5. Curvas de resistencia del modelo de turbulencia K-Omega. ................... 194
Gráfica 5.6. Curvas de resistencia del modelo de turbulencia K-Omega SST. ........... 195
Gráfica 5.7. Curvas de resistencia del modelo de turbulencia Smagorinsky. ............. 195
Gráfica 5.8. Curvas de resistencia del modelo de turbulencia K-Energy Two Layers.196
Gráfica 5.9. Comparación de las resistencias obtenidas por los distintos modelos de
turbulencia seleccionados...................................................................................... 200
Gráfica 5.10. Curvas de la resistencia viscosa y por formación de olas respecto a la
resistencia total. ..................................................................................................... 202
Gráfica 5.11. Perfil de olas a una velocidad de 3 nudos.............................................. 206
Gráfica 5.12. Perfil de olas a una velocidad de 5 nudos.............................................. 210
Gráfica 5.13. Perfil de olas a una velocidad de 10 nudos............................................ 214
Gráfica 5.14. Perfil de olas a una velocidad de 15 nudos............................................ 218
Gráfica 5.15. Perfil de olas a una velocidad de 20 nudos............................................ 222
Gráfica 5.16. Perfil de olas a una velocidad de 20 y 10 nudos.................................... 223
Gráfica 6.1. Comparación de las curvas de resistencia obtenidas mediante los distintos
métodos de predicción de resistencia. ................................................................... 228
Gráfica 6.2. Comparación de las curvas de resistencia obtenidas mediante los distintos
métodos de predicción de resistencia. ................................................................... 230
20
PRÓLOGO
21
CAPÍTULO I
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
1. INTRODUCCIÓN
23
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
da nunca en los buques mercantes, de pasaje o de servicio, en las que el máximo
peso muerto se da a la salida de puerto.
La estiba de las mercancías se produce, asimismo, en la mar, y frecuentemente
en condiciones adversas, lo que implica que la misma no está sometida a control
alguno por personal ajeno al propio buque, lo que puede conllevar una
disminución peligrosa de la estabilidad.
Aunque ya se ha mencionado anteriormente, es importante destacar que durante
una gran parte de su permanencia en la mar, el buque de pesca tiene su
maniobrabilidad restringida, especialmente los buques de arrastre. En ello se
parecen a los remolcadores, cableros, etc., pero se diferencian de los restantes.
Un factor de gran importancia, y también diferenciador con respecto a otro tipo
de embarcaciones, es que la tripulación suele participar, de una manera u otra,
del negocio, lo que implica que se cuide más su habilidad como ‘productor’ que
como ‘navegante’. Es por ello, que los conocimientos marineros y de navegación
son, con frecuencia, inferiores a los de los tripulantes de buques mercantes,
influyendo de manera decisiva en la seguridad del buque.
El tiempo de permanencia en la mar, sin tocar puerto, puede ser muy alto, de
varios meses para muchos pesqueros de gran altura. Ello debe ser considerado a
la hora de diseñar el buque, teniéndose que prestar una especial atención a la
habitabilidad, capacidad de combustible y otros consumibles, reparaciones y
facilidad de mantenimiento.
La propulsión del buque de pesca debe ser estudiada para diferentes formas de
operación, principalmente en los buques de arrastre, lo que obliga al proyectista
a buscar un compromiso entre las diferentes alternativas.
Por último, destacar que en cierto tipo de buques se da la característica de que
actúan en flotillas, particularidad que sólo es usual en buques de guerra.
24
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
2.1 Arrastrero
25
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
Los ramperos son los más importantes en número y suelen ser barcos con bastante
puntal y gran obra viva, necesaria para que las hélices puedan ejercer fuerza en el agua
para poder arrastrar los artes por el fondo. En el diseño del casco de un arrastrero,
priman las características hidrodinámicas, en particular la adecuada relación
eslora/manga, que posibiliten el desplazamiento, al, o desde el caladero en el menor
tiempo posible.
La caseta de gobierno se sitúa a proa para que la popa quede libre para la
maniobra, precisando menos personal y evitando riesgos de accidentes. Con esta
disposición el área sobre la que se realiza la faena además de estar protegida esta menos
sujeta al cabeceo. Estas y otras ventajas, como el que la embarcación pueda mantenerse
pescando cuando va recuperando la red, ha hecho que en los últimos años el número de
arrastreros por el costado no haya crecido, cediendo el primer puesto a los que lo hacen
por la popa.
Figura 1.1. Maniobra de arrastre por popa (Fuente: Micó Barba, F., 1978, Introducción a la pesca
marítima, publicación del ETSIN, Madrid).
26
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
27
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
2.2 Cerquero
Figura 1.3. Maniobra de cerco (Fuente:
http://estaticos.ipmedios.com/media/1/3/360bbc1a7232998b4718617309e70610-original.jpg).
28
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
1
copo en donde queda atrapado el cardumen ; después, se va cobrando el arte por uno o
varios extremos, ayudándose por medio de plumas, hasta que los peces capturados
quedan en un espacio mínimo; éstos se suben a bordo con un ‘brailer’ o salabardo2,
descargándolos sobre una tolva montada en cubierta, de la que por medio de canaletas
se dirigen hacia la escotilla de la bodega correspondiente.
Este tipo de pesca se utiliza para aquellas especies que navegan en bancos como el
atún, el salmón o la anchoa.
Figura 1.4. Buque atunero ‘Euskadi Alai’ (Fuente: http://i.ytimg.com/vi/IW5PS8p-
Ejc/maxresdefault.jpg).
1
Concentración grande de peces, generalmente de la misma especie, que se desplazan juntos.
2
Pequeño saco de red cuya boca va sujeta a un aro unido a un mango, que se emplea para extraer la pesca de las
redes grandes.
29
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
2.3 Palangrego
Figura 1.5. Cobrado de la línea por popa (izquierda) y por el costado (derecha) (Fuente: Micó Barba, F.,
1978, Introducción a la pesca marítima, publicación del ETSIN, Madrid).
30
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
sotavento. Si los sedales están amarrados a la línea principal, se desenganchan los peces
a medida que llegan a bordo. Cuando el buque cobre el aparejo por la popa, la maniobra
es similar.
Figura 1.6. Buque palangrero ‘Novo Airiños’ (Fuente:
http://www.insenaval.com/images/productos/buques_pesca/palangrero_50_metros_l.jpg).
31
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
2.4 Nasero
Figura 1.7. Buque nasero ‘Northwestern’ (Fuente:
http://r.ddmcdn.com/w_2594/s_f/o_1/cx_0/cy_0/cw_2594/ch_1729/DSC/uploads/2014/05/01-
northwestern-2600.jpg).
3
Máquina empleada para subir los aparejos de pesca a bordo de la embarcación.
32
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
33
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
El hecho de que el precio de realización de un ensayo sea independiente del
tamaño y naturaleza del barco, implica un encarecimiento relativo de esta herramienta
de diseño, a medida que el coste total del mismo disminuye, esto es, en proporción a su
tamaño y a la sofisticación del equipo instalado.
Por lo general, todos los elementos que constituyen un conjunto están de alguna
forma ligados entre sí; y, por lo tanto, no es extraño pensar que cuanto menor sea la
inversión precisa para construir un buque, tanto menor será el presupuesto disponible
para su diseño, y por ende el ‘cuidado’ y ‘dedicación’ que se ponga en el mismo.
Siguiendo el mismo hilo de razonamiento, es lógico suponer que aquellas Oficinas
Técnicas que se dedican a proyectar buques de un cierto porte (petroleros, graneleros,
etc.), o con una gran dosis de sofisticación (frigoríficos, portacontenedores, pasaje, etc.),
han de poseer una mayor cantidad de ‘experiencia’ debido fundamentalmente al hecho
de que están dotadas de una estructura suficiente como para permitirles mantener
Técnicos especializados en diferentes áreas.
Teniendo en cuenta que únicamente en muy contadas ocasiones el tamaño o el
grado de sofisticación del equipo de un buque de pesca justifica un coste total del
mismo lo suficientemente elevado como para que la experimentación física resulte
‘rentable’, en lo que concierne a la importancia relativa de tal partida en el coste total
del proyecto; la herramienta habitualmente más utilizada por los proyectistas de formas
de esta clase de barcos, es la ‘experiencia’, que en muchas ocasiones se traduce como el
empleo de una razonable cantidad de ‘intuición’ y ‘buen ojo’ para ‘deformar’ las líneas
de un buque base, a fin de adaptar su diseño a unas nuevas dimensiones.
En el caso de que las características intrínsecas del proyecto permitieran la
realización de ensayos de canal, bien sea porque el coste de producción del barco sea
elevado, o bien porque se pretenda construir una serie de buques semejantes; las
posibilidades de utilizar la experimentación física como medio para optimizar la carena
siguen siendo relativamente pequeñas. En efecto, no es difícil pensar que la realización
de un completo estudio hidrodinámico del comportamiento en aguas tranquilas de un
buque, encaminado a la optimización de sus formas, implica, por lo general, llevar a
cabo ensayos con al menos tres modelos de carena alternativos; representando su coste
una inversión suficientemente importante como para que en muy pocas ocasiones pueda
ser asumida. La evaluación durante la fase de diseño de otros aspectos relativos a las
prestaciones en servicio del buque, tan importantes como su comportamiento en la mar
o su maniobrabilidad, se efectúa únicamente en muy raras ocasiones mediante la
realización de ensayos de canal.
A modo de resumen, se puede concluir que las formas de los buques de pesca se
suelen diseñar básicamente mediante la adaptación de proyectos anteriores a unas
nuevas dimensiones, utilizando generalmente la experimentación física tan sólo para
comprobar sus prestaciones propulsivas en aguas tranquilas, en muchas ocasiones
únicamente por motivos legales, sin que a veces sus resultados puedan ser utilizados
para mejorar el diseño, ya sea porque se conocen en una fase muy avanzada del mismo
como para retocar las formas, o bien porque la información obtenida sea incompleta
debido a la realización de un exiguo número de ensayos.
La resistencia al avance de un barco se compone fundamentalmente de dos partes,
una debida a la viscosidad del fluido que, al contacto con la superficie exterior de la
obra viva, genera una fuerza que opone al movimiento y otra que es debida a la
gravedad y que se conoce con el nombre de resistencia por formación de olas. Esta
última se pone de manifiesto, de ahí le viene el nombre, por el sistema de olas que
34
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
acompaña al barco. Su valor puede determinarse, lo que no es tan fácil, calculando la
energía gastada en generar el sistema de olas.
La magnitud relativa de cada una de las componentes de la resistencia al avance
de un barco depende de la velocidad a la que se navega y de las formas de la carena.
Para valores pequeños de la velocidad la resistencia ocasionada por la viscosidad
predomina sobre la resistencia por formación de olas. La resistencia viscosa puede
llegar a ser hasta un 80% de la resistencia total. A medida que la velocidad crece la
resistencia por formación de olas también lo hace, mientras que la que tiene su origen en
la viscosidad lo hace en menor proporción. Cuando la velocidad llega a ser
suficientemente alta, el valor que toma la resistencia debida a la formación de olas es
tan alto que supone una barrera infranqueable para carenas de formas convencionales.
Cuando se trata de estudiar la influencia de las formas geométricas de la carena
sobre la resistencia es necesario tener presente que el sistema de olas se produce por que
el barco, al pasar, se abre camino empujando al agua y perturbando el equilibrio de la
superficie libre. Podría decirse que cuando la embarcación navega en la superficie del
mar amontona el agua; como consecuencia la acción de la aceleración gravitatoria la
superficie libre tiende a su estado de equilibrio inicial, lo que se traduce en una
generación de olas. Cuanto mayor es la perturbación que el barco produce en la
superficie libre mayor será el sistema de olas y por consiguiente mayor la resistencia.
Matemáticamente el fenómeno se estudia suponiendo que lo que el barco hace al
cortar la superficie es dar, de forma continuada especialmente con la proa (aunque en
teoría lo hace con toda su obra viva), impulsos de presión. La integración de los
sistemas de olas que cada uno de estos impulsos genera da el sistema de olas generado
por el barco.
La base del proyecto hidrodinámico de una embarcación rápida está en hacer que
la perturbación que, al navegar, produce sobre la superficie del agua sea lo menor
posible. Por esta razón los buques rápidos tienen proas y cuerpos de entradas muy finos,
de manera que el área de las secciones cambie suavemente y la separación del agua por
parte del casco esté producida principalmente por la velocidad. El lanzamiento de la
proa, característica común en los barcos rápidos, tiene la misma razón de ser aunque en
este caso es para evitar que la proa altere el campo de presión simultáneamente en todas
las líneas de agua. Con esto la proa corta al agua más suavemente y la perturbación es
menor. La finura de la proa hace que en los barcos rápidos el desplazamiento del cuerpo
de proa sea menor que el del cuerpo de popa y por tanto que el centro de empuje está a
popa de la sección media. La importancia de la resistencia por formación de olas a altas
velocidades y su dependencia de la perturbación creada en la superficie libre pone de
manifiesto la necesidad de cuidar al máximo el proyecto de las formas en este tipo de
barcos. Un error en las líneas del casco de un barco rápido se paga muy caro, al menos
desde el punto de vista hidrodinámico.
En los barcos lentos predomina la resistencia debida a la viscosidad por lo que es
bueno proyectar las formas de manera que esa resistencia sea lo menor posible. Por
tanto conviene que la superficie del casco que está en contacto con el agua, la que se
conoce como superficie mojada, sea la menor posible. Otro fenómeno que tiene su
origen en la viscosidad y que lleva consigo un fuerte aumento de resistencia es la
separación de la capa límite. Se produce casi siempre en la popa y se pone de manifiesto
por una serie de torbellinos claramente visibles cuando se mira la estela que va dejando
el barco. Para evitar la separación es conveniente hacer más suave el flujo del agua en la
popa aumentando la longitud del cuerpo de salida. En este caso el desplazamiento se
35
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
lleva hacia proa y el centro de carena se sitúa a proa de la cuaderna media: lo contrario
de lo que pasaba cuando el barco era rápido.
Aunque la variedad de buques de pesca es muy grande, por lo que no es posible
dar reglas generales, lo que les caracteriza desde el punto de vista de la resistencia al
avance es que todos ellos tienen la necesidad, o al menos eso es lo que quieren sus
armadores, de alcanzar el caladero y regresar de él lo antes posible. Por consiguiente
son barcos a los que se le exige navegar con rapidez. Conviene recordar que lo que
define si una embarcación es rápida no es el valor absoluto de la velocidad sino una
cierta relación entre ella y la eslora, el número de Froude ()ݎܨ.
Un barco de pesca se mueve a números de Froude próximos a 0,4 por lo que la
componente de la resistencia al avance debida a las olas tiene un valor significativo. Son
por tanto, como antes se ha dicho, barcos de los considerados semi-rápidos y tiene, o al
menos deben tener, el centro de carena a popa de la cuaderna maestra. Hay que repetir
sin embargo que la variedad de los barcos de pesca impide dar normas indiscutibles
sobre lo que es mejor, ya que cada caso es distinto. Es también cierto que, debido a
tener que permanecer durante largos periodos de tiempo en la mar, necesitan un
excelente comportamiento en mares adversos. Por eso se admite que el centro de carena
pueda estar un poco a proa de la posición óptima desde el punto de vista de la
resistencia y que se les dé más manga de la que les correspondería si se tuvieran en
cuenta únicamente criterios relacionados con la resistencia al avance.
Muchos han sido los intentos que a lo largo de los años se han hecho para mejorar
las características hidrodinámicas de los barcos de pesca. Según José A. Aláez Zazurca
en su publicación “La Hidrodinámica del buque de pesca” [2], el aumento progresivo
de la velocidad a la que navegaban hizo que las formas evolucionaran, abandonándose
en el cuerpo de proa las antiguas formas en ‘U’ y adaptándose las formas en ‘V’, a la
vez que se afinaba el ángulo de entrada en la flotación. Para hacer el cabeceo más suave
a las cuadernas de proa se les da un abanico apreciable lo que permite mantener seca la
cubierta aunque se navegue con mal tiempo.
Figura 1.8. Evolución de las formas de las cuadernas en las embarcaciones de pesca. (Fuente:
http://bertan.gipuzkoakultura.net/23/argazkiak/g/161.jpg).
36
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
sencillez geométrica y el relativamente bajo coste de construcción. Como por razón de
la velocidad, es conveniente que el desplazamiento del cuerpo de proa sea menor que el
del cuerpo de popa, los pesqueros tienen la cuaderna maestra, que es la que tiene mayor
área, situada ligeramente a popa de la media. Suelen tener también astilla muerta cuyo
valor debe elegirse acorde con el radio de pantoque para evitar que se formen
torbellinos en esa zona.
Figura 1.9. Evolución de las proas y popas de las embarcaciones de pesca. (Fuente:
http://bertan.gipuzkoakultura.net/23/argazkiak/g/231.jpg).
4
Golpe o impacto en la embarcación cuando el extremo de la misma entra en el agua (tras haberse elevado la proa o
popa previamente por encima del agua) que produce un súbito cambio en la aceleración vertical del buque.
37
CAPÍTULO I - INTRODUCCIÓN A LOS BUQUES DE PESCA
DESCRIPCIÓN, CARACTERÍSTICAS GENERALES, FORMAS TÍPICAS Y OTRAS PARTICULARIDADES
Figura 1.10. Representación del “slamming” de un buque. (Fuente: Aláez Zazurca, J. A., 30 de marzo de
1995, Comportamiento del buque en la mar, Madrid).
38
CAPÍTULO II
1. INTRODUCCIÓN
40
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
modelo a ensayar, por lo que se suele recurrir a modelos matemáticos (Computational
Fluid Dynamics, CFD), el cual ya se verá más adelante, o métodos analíticos para hacer
las primeras aproximaciones de la resistencia al avance.
Para aplicar uno de estos métodos se deberá elegir entre la información
disponible, poniendo especial interés en que las características de nuestro buque
respondan a las características de los buques basándose en los cuales se desarrolló el
citado método.
Estos métodos no sustituyen, ni mucho menos, a los métodos experimentales ya
que estos son más precisos en el cálculo de la resistencia, y además se dispondrá de
otros muchos datos para el correcto diseño de la propulsión, de la disposición de
apéndices, etc. pudiéndose hacer pequeñas modificaciones de las formas que mejoren el
comportamiento hidrodinámico de la carena.
Entre los métodos analíticos podemos distinguir las series sistemáticas de formas
y los métodos numéricos.
a) Serie sistemática
Una serie sistemática se desarrolla para optimizar las formas de la carena de un
determinado tipo de buque. Se seleccionan unas formas ‘patrón’ y, a partir de ellas, se
realizan variaciones sistemáticas, tanto de sus dimensiones principales como de sus
coeficientes geométricos adimensionales. Paralelamente se han de seleccionar tres o
cuatro coeficientes de carena, o de relaciones entre sus dimensiones principales, que
tengan mayor influencia en la resistencia al avance. Se construyen una serie de modelos
de la serie y se realizan los ensayos para medir la resistencia al avance, que pueden
verse complementados por ensayos de carenas con apéndices, ensayos de líneas de
corriente, ensayos de autopropulsión, etc.
Los resultados de los ensayos se publican de forma que se facilite la interpretación
y aplicación a los futuros usuarios. La idoneidad de una serie sistemática será máxima
cuando el buque estudiado prevea formas similares a las estudiadas en la serie y las
relaciones adimensionales de las dimensiones principales del buque estén dentro de los
rangos de variación cubiertos por la serie. Por el contrario, cuanto más diferencia haya
ira disminuyendo la fiabilidad de los resultados.
b) Métodos numéricos
Por otro lado están los métodos numéricos, los cuales la mayoría de ellos para
predecir la resistencia al avance tienen base estadística.
Para ello primero se elige una base de datos apropiada evitando que contenga
casos singulares o información demasiado dispar. A continuación se eligen la variables
independientes con mayor influencia en la magnitud que se está tratando de estimar
(resistencia). Posteriormente se selecciona el tipo de función de dependencia que
relacione la magnitud a estudiar con las variables independientes elegidas y, finalmente,
se realizan ajustes de la curvas para la obtención de los coeficientes de las funciones
empleando las técnicas más adecuadas para cada tipo de función.
La base de datos puede formarse con resultados de buques sin ninguna relación
entre sí, o con buques pertenecientes a un determinado tipo. La predicción basada en
una regresión será tanto más buena como lo sea la base de datos empleada, dependiendo
41
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
la precisión del grado de semejanza entre los buques de la muestra y la del buque
considerado.
42
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
2. MÉTODO DE J. HOLTROP Y G. G. J. MENNEN
43
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
0,06
1 , (2.5)
4 1
en esta fórmula es el coeficiente prismático referido a la eslora en la flotación,
es la posición longitudinal del centro de carena a partir de la sección media como
porcentaje de ( será positivo hacia proa).
El coeficiente pretende introducir la forma de popa en el factor de forma:
1 0,011 (2.6)
Cstern
Popas tipo góndola -25
Cuadernas en V -10
Cuadernas normales 0
Cuadernas en U con popa Hogner 10
Tabla 2.2. Valores del coeficiente Cstern.
44
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
Tipo de apéndice (1 + k2)
Henchimientos protectores 3,0
Henchimientos integrados 2,0
Ejes 2,0 – 4,0
Aletas estabilizadoras 2,8
Domos 2,7
Quillas de balance 1,4
Tabla 2.3. Valores de (1 + k2) para varios tipos de apéndices.
- Para 0,4
(2.10)
,
, , (2.11)
2223105 90 ,
/ 0,11 ,
0,229577 /
0,11 / 0,25 /
/ 0,25 0,5 0,0625 /
Tabla 2.4. Valores del coeficiente C7.
0,9 (2.13)
, (2.14)
,
0,56
, (2.15)
0,31
siendo el área transversal del bulbo en la vertical de la intersección de la flotación
con el perfil de proa (m2) y la altura sobre la base del c.d.g. de la sección del bulbo
de proa en la perpendicular de proa (m).
0,8
1 (2.16)
45
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
con:
, (2.18)
0,4
con:
/ 512 1,69385
512 / 1727 / / 8,0
1,69385
2,36
/ 1727 0
Tabla 2.6. Valores del coeficiente C15.
/B 12 1,446 0,03
/B 12 1,446
Tabla 2.7. Valores del coeficiente λ.
- Para 0,55
(2.19)
, ,
, (2.20)
6919,3 2
, ,
7,2035 (2.21)
, ,
0,4 (2.22)
46
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
,
(2.24)
0,11 ,
donde el coeficiente es una medida para la parte que emerge de la proa y es el
número de Froude referido a la inmersión:
0,56
(2.25)
1,5
(2.26)
0,25 0,15
5 0,2 1 0,2
5 0
Tabla 2.8. Valores del coeficiente C6.
2 (2.28)
47
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
/ 0,04 /
/ 0,04 0,04
Tabla 2.9. Valores del coeficiente C4.
Si se desea predecir el valor de para una rugosidad mayor de 150 μm, se puede
utilizar la fórmula:
/
0,105 0,05579
∆ /
(2.33)
48
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
3. MÉTODO DE AMADEO GARCÍA GÓMEZ
Característica Rango
Eslora entre perpendiculares, (m) 25 – 60
Número de Froude, 0,25 – 0,40
Relación / 0,095 – 0,165
Tabla 2.10. Rango de aplicación del método de Amadeo García.
(2.39)
1,24
0,265 2,151 0,298 (2.40)
49
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
(2.41)
1
La superficie mojada sin apéndices, , puede estimarse por:
0,5 (2.44)
1 (2.45)
0,07 20
La superficie mojada del timón y del codaste cerrado puede estimarse por:
0,1 , (2.46)
donde es el calado medio (m).
(2.47)
%
1
100
% (2.48)
, (2.51)
La superficie mojada sin apéndices para pesqueros con bulbo de proa, , puede
estimarse por:
(2.53)
1
50
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
4. MÉTODO DE VAN OORTMERSSEN
Característica Rango
Eslora en la flotación, (m) 8 – 80
Volumen de carena, m2 5 – 3000
Relación eslora/manga, / 3 – 6,2
Relación manga/calado, / 1,9 – 4,0
Coeficiente prismático, 0,50 – 0,725
Coeficiente de la maestra, 0,73 – 0,97
Posición longitudinal del centro de carena, -8% LD – +2,8% LD
Semiángulo de entrada en la flotación, (º) 10 – 46
Relación velocidad/eslora, / 0 – 1,79
Número de Froude, 0 – 0,50
Tabla 2.11. Rango de aplicación del método de Van Oortmersen.
(2.54)
2
La resistencia total viene dada por
, (2.55)
donde
1
(2.56)
2
0,075
(2.57)
2
1 1
(2.58)
∆
,
0,14347 (2.59)
51
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
10 , , , , , , ,
(2.60)
, , , , ,
(2.61)
i 1 2 3 4
di,0 79,32134 6714,88397 -908,44371 3012,14549
di,1 -0,09287 19,83 2,52704 2,71437
di,2 -0,00209 2,66997 -0,35794 0,25521
di,3 -246,45596 -19662,024 755,1866 -9198,8084
di,4 187,13664 14099,904 -48,93952 6886,60416
di,5 -1,42983 137,33613 9,86873 -159,92694
di,6 0,11898 -13,36938 -0,77652 16,23621
di,7 0,15727 -4,49852 3,7902 -0,82014
di,8 -0,00064 0,021 -0,01879 0,00225
di,9 -2,52862 216,44923 -9,24399 236,3797
di,10 0,50619 -35,07602 1,28571 -44,1782
di,11 1,62851 -128,72535 250,6491 207,2558
Tabla 2.12. Valores de los coeficientes di,j de la regresión.
1
(2.62)
2
Corrección 103 CA
Rugosidad 0,35
Gobierno 0,04
Quillas de balance 0,04
Aire 0,08
Tabla 2.13. Valores del coeficiente CA según distintas correcciones.
52
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
5. MÉTODO DE SANDER M. CALISAL Y DAN MCGREER
Característica Rango
Coeficiente de bloque, 0,531 – 0,615
Relación eslora/manga, / 2,6 – 4,0
Relación manga/calado, / 2–3
Tabla 2.14. Rango de aplicación del método de Sander M. Calisal y Dan McGreer.
, (2.67)
donde
, , , , , (2.68)
,
0,14347 (2.69)
Los valores de , fueron determinados mediante análisis de regresión de datos de
ensayos con modelos, los cuales vienen dados para su correspondiente coeficiente de
bloque . Para valores intermedios de dicho coeficiente será necesario realizar una
interpolación.
53
CAPÍTULO II - INTRODUCCIÓN A LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA
CB = 0,531
i 1 2 3 4
di,0 0,006056 0,166001 -0,067109 -0,047944
di,1 1,42E-06 -0,000817 0,000425 0,000326
di,2 9,99E-06 -0,005752 0,00299 0,002295
di,3 -0,001758 -0,00068 0,000302 0,000147
di,4 0,000246 -0,002688 0,001521 0,00118
Tabla 2.15. Valores de los coeficientes di,j para un coeficiente de bloque CB = 0,531.
CB = 0,615
i 1 2 3 4
di,0 0,074654 0,076958 -0,000162 -0,018642
di,1 0,001879 7,70E-05 7,42E-05 1,05E-05
di,2 -0,000701 1,80E-05 0,000907 0,000119
di,3 -0,05158 -0,005247 -0,001137 -0,001944
di,4 0,009871 0,000572 -0,000661 0,001169
Tabla 2.16. Valores de los coeficientes di,j para un coeficiente de bloque CB = 0,615.
(2.71)
Los valores de los coeficientes vienen recogidos en la siguiente tabla para plena
carga y carga parcial:
54
CAPÍTULO III
Generación de formas
CAPÍTULO III - GENERACIÓN DE FORMAS
1. INTRODUCCIÓN
En este capítulo se va a analizar las formas que se van a adoptar así como la
generación de las mismas. Como en todos los buques las formas del casco son, como ya
se ha expuesto en el Capítulo I, función del tipo de mercado al que va destinado el
barco, así como de la carga (tipo y volumen o peso), y de la zona de navegación y
puertos inicialmente establecidos como previsibles, al igual que los canales por los que
previsiblemente pudiese navegar. Otro punto importante a tener en consideración es la
meteorología de la zona donde llevará a cabo su servicio.
Hay que señalar que, como en todo proyecto, existen unas determinadas
dimensiones que se deben optimizar con el fin de:
1. Conseguir que las formas del casco tengan la menor resistencia al avance
posible, lo que se traduce en la necesidad de menor potencia y, por tanto, el
motor propulsor será de menor potencia, lo que se traduce en un menor tamaño
del motor propulsor, con lo que disminuye el consumo y el precio del mismo.
Así mismo se reduce su tamaño con lo que disminuye el empacho, dejando un
mayor espacio para albergar la carga, así como el peso en rosca del buque.
2. Lograr una buena estabilidad, suficiente para la seguridad de la embarcación,
pero no excesiva, por motivos de confortabilidad y de comportamiento en la
mar, con lo que también se aumenta la seguridad del buque y de la tripulación.
3. Dotar a la embarcación de espacios sobre cubiertas suficientes para conseguir
una disposición general óptima, característica muy importante en este tipo de
buques factoría, pues deben estar pensados no sólo para navegar y transportar
carga, sino para llevar a cabo un proceso productivo a bordo (pesca. preparación,
procesado y congelación de las capturas).
4. Aumentar en lo posible la capacidad de carga, ya que a mayor volumen de
bodegas, mayor es la cantidad de capturas que se pueden almacenar en la misma
marea, aumentando la rentabilidad del buque. La maximización de la capacidad
de carga debe realizase dentro de un rango razonable, ya que al margen de
razones constructivas se debe de tener en cuenta las limitaciones de pesca
impuestas en los caladeros, con el fin de permitir la regeneración de las especies.
5. Las formas elegidas deben favorecer la mejor maniobrabilidad posible, aspecto
muy importante durante las faenas de pesca.
56
CAPÍTULO III - GENERACIÓN DE FORMAS
Además es conveniente el uso de popas de espejo en aquellos pesqueros que
necesitan amplios espacios de trabajo en esa zona y es importante que, a ser posible, el
barco tenga un codaste amplio, para poder meter una hélice de gran diámetro, que gire a
pocas revoluciones y que tenga inmersión suficiente (es aconsejable dotar a la
embarcación de un cierto trimado de proyecto para aumentar la inmersión del
propulsor).
La posición longitudinal de la sección de área máxima tiene una gran influencia
en la eficiencia del casco. Para una resistencia óptima, la sección de área máxima en
buques pesqueros debe estar algo a popa de la sección media.
A la hora de definir las formas del barco, podemos abordar el proyecto por
diversos caminos diferentes; pero fundamentalmente se obtiene de las siguientes tres
maneras:
Para abordar el análisis de las formas se empleará este último modo, es decir, se
partirá de un buque arrastrero ya construido, aunque previamente se realizará un análisis
de sus formas.
57
CAPÍTULO III - GENERACIÓN DE FORMAS
58
CAPÍTULO III - GENERACIÓN DE FORMAS
Relación L/B
Aunque los pesqueros son buques pequeños y rápidos, con lo que en principio les
corresponderían relaciones / muy altas, a su vez deben tener mucha manga por
consideraciones de capacidad de bodegas y sobre todo de estabilidad, con lo que / se
ve comprendida entre 4,0 y 5,0. El buque base tiene una relación
50,5
4,208
12
Relación B/T
Los valores típicos de / para pesqueros están entre 2,2 y 2,5. El buque base
tiene una relación
12
2,5
4,8
Curva de Áreas
Se comprueba que la curva de áreas del buque base tenga unos hombros con radio
de curvatura amplio, para que no haya desprendimientos de flujo en el hombro de popa
ni creación de sistemas secundarios de olas en el hombro de proa. También se verifica la
concavidad en la zona de proa.
14 0,51445
0,324
9,81 50,5
Línea de flotación
La de flotación es la más importante de las líneas de agua, pues los fenómenos
hidrodinámicos más importantes se producen cerca de la flotación. Para números de
Froude comprendidos entre 0,25 y 0,40 (caso del buque base) la flotación a proa debe
ser ligeramente cóncava, y así se comprueba al observar la curva de áreas del buque
base y que será la que se adopte para el buque proyecto.
59
CAPÍTULO III - GENERACIÓN DE FORMAS
Forma de las cuadernas
La zona de proa se proyecta con unas formas en ‘U’, con el fin de amortiguar los
movimientos de cabeceo, así como de un amplio abanico para reducir el embarque de
agua en cubierta.
La astilla muerta en la cuaderna maestra es amplia para obtener un plano de deriva
suficiente, manteniendo bajo el coeficiente de bloque, de manera que se favorece la
maniobrabilidad a la vez que se reduce la resistencia al avance.
La popa será en espejo, minorando el ángulo de salida de las líneas de agua cerca
de la flotación, reduciendo así el riesgo de formación de fenómenos de flujo indeseables
(desprendimiento de flujo, turbulencias.). Además de lo anterior, una popa de espejo
proporciona una superficie útil en popa mucho mayor que las popas de crucero, por lo
que se facilitan las operaciones de pesca.
El vano de la hélice deberá ser lo mayor posible pero permitiendo, en todo
momento, una inmersión suficiente del propulsor, con el fin de favorecer el rendimiento
propulsivo.
El ángulo de salida de las líneas de agua en la zona de la hélice será lo menor
posible para facilitar la llegada de agua a la hélice, mejorar la maniobrabilidad y reducir
el coeficiente de succión.
Se adoptaran unas formas redondeadas que presentan unas condiciones marineras
globales bastantes buenas y permiten un buen aprovechamiento de los espacios
interiores. Desde otro punto de vista, las formas redondeadas presentan una menor
resistencia en olas y menores momentos flectores.
60
CAPÍTULO III - GENERACIÓN DE FORMAS
3. OBTENCIÓN DE LA GEOMETRÍA
En este apartado se realizará una descripción de todos los pasos a seguir para
llevar a cabo la generación de la geometría a partir del buque base escogido
anteriormente.
Inicialmente los datos de los que se disponen es el plano de formas del barco base,
el cual se muestra a continuación. Para la realización de la superficie del casco se utiliza
un programa de diseño CAD (concretamente Rhinoceros). Para ello se importa a
Rhinoceros cada una de las vistas del plano de formas (alzado, planta y perfil) como
‘bitmap de fondo’ de forma individual, escalados según las dimensiones del barco.
61
CAPÍTULO III - GENERACIÓN DE FORMAS
Figura 3.1. Plano de formas del buque base (Fuente: Velázquez Sánchez, J. A., octubre de 2010, Arrastrero congelador factoría de fondo y pelágico, ETSIN, Madrid).
62
CAPÍTULO III - GENERACIÓN DE FORMAS
Posteriormente se procede a ‘calcar’, es decir, a dibujar las distintas curvas que
forman el plano de formas. Una vez hecho esto se giran las curvas de las vistas
correspondientes (alzado, planta y perfil) y se llevan a su verdadera posición tanto las
cuadernas como las líneas de agua y las longitudinales, obteniendo finalmente el
‘esqueleto’ del barco.
Figura 3.2. ‘Esqueleto’ del casco del buque base (Fuente: Rhinoceros).
Figura 3.3. Superficie del casco del buque base (Fuente: Rhinoceros).
63
CAPÍTULO III - GENERACIÓN DE FORMAS
Figura 3.4. Superficie del casco afín (Fuente: Maxsurf).
Tras esta transformación afín se procede a la realización del plano de formas del
buque final con la ayuda del programa Rhinoceros, el cual se puede ver a continuación.
64
L
C
Sec. 16 Sec. 17 Sec. 18 Sec. 19 Ppr L. 8 L. 6 L. 4 L. 2 L. 2 L. 4 L. 6 L. 8
L. 9 L. 7 L. 5 L. 3 L. 1 L. 1 L. 3 L. 5 L. 7 L. 9
Ppp Sec. 1 Sec. 2 Sec. 3 Sec. 4 Sec. 5 Sec. 6 Sec. 7 Sec. 8 Sec. 9 Sec. 10 Sec. 11 Sec. 12 Sec. 13 Sec. 14 Sec. 15
L. de Agua 8
20
19
L. de Agua 7 L. de Agua 7
18
17
L. de Agua 6 L. de agua 6
16
.9
ng
15
Lo
L. de Agua 5 L. de Agua 5
2 0
8
.
.6
ng
.7
.
.1
14
ng
ng
3
Lo
ng
ng
g.
.
Lo
Lo
1
ng
n
Lo
Lo
Lo
Lo
L. de Agua 4 L. de Agua 4
.4
2
ng
Lo
3 13
L. de Agua 3 L. de Agua 3
4 12
5
L. de Agua 2 L. de Agua 2 6
11
7
8 10
L. de Agua 1 L. de Agua 1
9
L. 9 L. 7 L. 5 L. 3 L. 1 L. 1 L. 3 L. 5 L. 7 L. 9
Sec. 1 Sec. 2 Sec. 3 Sec. 4 Sec. 5 Sec. 6 Sec. 7 Sec. 8 Sec. 9 Sec. 11 Sec. 12 Sec. 13 Sec. 14 Sec. 15 Sec. 16 Sec. 17 Sec. 18 Sec. 19 L. 8 L. 6 L. 4 L. 2 L. 2 L. 4 L. 6 L. 8
Ppp Ppr
Sec. 1 Sec. 2 Sec. 3 Sec. 4 Sec. 5 Sec. 6 Sec. 7 Sec. 8 Sec. 9 Sec. 11 Sec. 12 Sec. 13 Sec. 14 Sec. 15 Sec. 16 Sec. 17 Sec. 18 Sec. 19
L
C
CARACTERÍSTICAS PRINCIPALES
Eslora total ___________________________________ 48,421 m.
Eslora entre perpendiculares _____________________ 44,490 m.
Eslora de la flotación ____________________________ 44,879 m.
Long. 9 Long. 9 Calado de proyecto _____________________________ 4,000 m.
Manga flot. proyecto ____________________________ 10,000 m.
Long. 8 Long. 8
e A. 4
L. d
Long. 7
L. d
Long. 7 CARACTERÍSTICAS DE TRAZADO
.3 e A. 5 Separación entre líneas de agua __________________ 1,000 m.
Long. 6 Long. 6
deA L. d
e A. Separación entre secciones ______________________ 2,225 m.
L. 6
Long. 5 .2 Long. 5 Separación entre longitudinales ___________________ 0,500 m.
eA .1
L. d
e
.d A. 7
Long. 4 L eA L.
de Long. 4
L .d A.
Long. 3 8 Long. 3 Nombre Firma
Long. 2 Long. 2 UNIVERSIDAD POLITÉCNICA
Alumno Daniel García López
Long. 1 Long. 1
DE CARTAGENA
L L
Director Domingo L. García López
C C Codirector José Enrique Gutiérrez Romero
Escuela Técnica Superior de Ingeniería Naval y Oceánica
4. CURVA DE ÁREAS
(3.1)
⨂
66
CAPÍTULO III - GENERACIÓN DE FORMAS
45
40
35
30
25
Área (m2)
20
15
10
0
‐2 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46
Eslora (m)
Gráfica 3.1. Curvas de áreas del buque para un calado de 4 m.
67
CAPÍTULO III - GENERACIÓN DE FORMAS
5. CURVAS HIDROSTÁTICAS
5.1 Introducción
Área de la flotación
Es la curva básica para obtener todas las demás. Para un calado dado, el área de la
flotación se obtiene integrando las distintas mangas, o semimangas (pues lo normal es
que la flotación sea simétrica), a lo largo de la eslora en la flotación correspondiente a
ese calado
2 , (3.2)
(3.4)
100
La suposición que se hace aquí es considerar que el área de la flotación permanece
constante en toda la variación de calado que se produzca. Por tanto las variaciones de
calado tendrán que ser pequeñas (a no ser que el buque sea de costados rectos) o el error
que se cometa puede ser muy grande.
68
CAPÍTULO III - GENERACIÓN DE FORMAS
⨂
⨂ (3.5)
⨂ 2 (3.6)
2 , (3.7)
Desplazamiento de trazado
El desplazamiento se calcula multiplicando el volumen de carena por el peso
específico supuesto constante
∆ (3.8)
El volumen de trazado, correspondiente a cada flotación, es el volumen que se
calcula con las lecturas del plano de formas o de la cartilla de trazado
(3.9)
(3.10)
, (3.11)
⨂ ⨂ , (3.13)
69
CAPÍTULO III - GENERACIÓN DE FORMAS
Radio metacéntrico transversal
El radio metacéntrico transversal es el radio de curvatura de la curva que describe
el centro de carena al girar el buque un ángulo infinitesimal en el plano transversal,
manteniendo constante el volumen sumergido, y se calcula mediante la expresión
, (3.14)
, (3.16)
siendo el momento de inercia del área de la flotación respecto del eje transversal que
pasa por su centro de gravedad ( ).
Como el centro de la flotación puede estar en cualquier posición resulta más
cómodo calcular primero el momento de inercia respecto del eje paralelo que pasa por la
sección media, ⨂ , y después, aplicando al teorema de Steiner, obtener
⨂ ⨂ (3.17)
Si la flotación es simétrica
⨂ 2 (3.18)
, (3.21)
100
70
CAPÍTULO III - GENERACIÓN DE FORMAS
siendo el mismo momento de inercia que se utiliza para calcular , el peso
específico supuesto para el agua de mar y la eslora entre perpendiculares.
Coeficientes de forma
Estos coeficientes son unas relaciones entre los volúmenes sumergidos o las
superficies del buque, y las de los cuerpos geométricos que ‘más o menos’ los
circunscriben.
Coeficiente de bloque ( )
Es la relación entre el volumen sumergido del buque y el volumen de un prisma
de dimensiones su eslora, manga y calado.
, (3.22)
Coeficiente de la maestra ( )
Es la relación entre el área sumergida de la sección media y el área del rectángulo
de dimensiones la manga y el calado.
⨂
, (3.23)
Coeficiente de la flotación ( )
A este coeficiente a veces se le llama de afinamiento de las líneas de agua, y se le
define como la relación entre el área correspondiente a una flotación dada y la del
rectángulo de dimensiones la eslora y la manga.
, (3.24)
71
CAPÍTULO III - GENERACIÓN DE FORMAS
Los valores de este coeficiente para el desplazamiento máximo pueden oscilar
entre 0,70, o incluso algo menos, para buques con extremos muy afinados, hasta 0,90
cuando tienen mucho cuerpo central.
(3.25)
⨂
(3.26)
Las tres variables que intervienen son las de trazado correspondientes a una
flotación dada. Los valores de este coeficiente para el desplazamiento máximo pueden
variar desde 0,70 o 0,71, para destructores o lanchas rápidas, hasta 0,96 para grandes
petroleros.
Un valor muy alto da una idea de un buque con secciones en forma de ‘U’,
mientras que si es bajo sus secciones tendrán forma de ‘V’.
Superficie mojada
La superficie mojada es el área de la superficie exterior de la carena que está en
contacto con el agua para una flotación determinada y se utiliza para calcular la
resistencia de fricción del buque durante la marcha.
72
CAPÍTULO III - GENERACIÓN DE FORMAS
5.2 Metodología y aplicación
Figura 3.5. Superficie del casco a estudiar hidrodinámicamente (Fuente: Maxsurf).
73
CAPÍTULO III - GENERACIÓN DE FORMAS
Nº NOMBRE LT LPP B T HCP HCS HCP – T
(m) (m) (m) (m) (m) (m) (m)
6 PLAYA DE LAGOS 35,25 29 8 3 3,6 5,7 0,6
7 URONDO 39 30,6 8,8 3,25 3,8 5,95 0,55
8 CACHACHO 38,85 32 8,4 3,05 3,7 5,9 0,65
9 MARDE MARES 38,85 32 8,4 3,05 3,7 5,9 0,65
10 BOGAVANTE SEGUNDO 38,4 32,4 8,4 3,5 3,9 6,05 0,4
11 MOHAMMED (5) 42 34,4 10 4 5 7 1
12 VILLA DE PITANXO 50 41 9,7 4 4,35 6,85 0,35
13 PATRICIA SOLTELO 51 42,5 9,8 4 4,3 6,75 0,3
14 RIO CAXIL 50,5 42 10 4 4,5 6,6 0,5
15 MOHAMMED (13) 52 45,8 9 4 5 7 1
16 ANA GANDÓN 58 48,75 10 4 4,5 6,7 0,5
17 FERRALEMES 59 50,65 11 4 5 7 1
18 HERMANDOS TOUZA 60 51,9 10,2 4 5 6,7 1
19 DABFAR 60 54 12 4,5 5 7,3 0,5
20 FONSECA (ex H. Garrido) 65 56,7 11 4 5 7 1
21 RIO NERVION 66,5 57 11,2 4,65 5 7,1 0,35
22 PASCOAL ATLANTICO 69 60 12 4 5 7 1
23 D. GUILLERMO GEFAELL (12) 92 80 15 6 7 10 1
Tabla 3.2. Base de datos de buques arrastreros.
74
CAPÍTULO III - GENERACIÓN DE FORMAS
5.3 Resultados
Una vez ya definidas las líneas principales sobre el casco se procede al cálculo de
las características hidrostáticas para el calado de proyecto, las cuales se pueden observar
en la siguiente tabla:
CARACTERÍSTICAS HIDROSTÁTICAS
Calado (T) 4 m
Desplazamiento (Δ) 1033 t
Eslora de la flotación (LF) 44,879 m
Manga máxima en la flotación (Bmax) 10 m
Superficie mojada (S) 572,547 m2
Área de la flotación (Af) 336,922 m2
Coeficiente prismático (CP) 0,606
Coeficiente de bloque (CB) 0,561
Coeficiente de la maestra (CM) 0,932
Coeficiente de la flotación (CF) 0,751
LCB desde PPP 21,675 m
LCF desde PPP 19,194 m
KB 2,28 m
BMt 2,204 m
BML 37,282 m
KMt 4,484 m
KML 39,562 m
TCI 3,453 t/cm
MTC 8,254 t·m/cm
Tabla 3.3. Características hidrostáticas del casco al calado de proyecto.
El que el buque posea un coeficiente de bloque relativamente alto nos indica que
es un buque de formas llenas, es decir, es ‘muy cuadrado’, que se puede observar en la
curva de áreas que se analizará en el siguiente apartado. Aunque en este caso tiene un
valor medio.
El coeficiente de la maestra es muy próximo a uno, lo que es normal debido a la
forma de la cuaderna maestra.
El coeficiente de flotación nos relaciona el área en la flotación con la manga y la
eslora del buque, por ello un coeficiente de flotación grande (próximo a uno) nos indica
que las formas en proa y popa para la flotación considerada (en este caso al calado de
diseño, 4 m) son muy llenas. Así en este caso se ve que en la zona de popa las formas
son muy llenas, como corresponde a una popa de espejo, mientras que a proa las formas
son más afinadas para que así la resistencia al avance del buque sea menor, por lo que la
potencia necesaria también lo será.
Por último la posición longitudinal del centro de carena es recomendable que este
a popa de la sección maestra, pero dentro de un margen que se ha fijado entre el 1% y el
2,5% respecto a la eslora entre perpendiculares del buque.
75
CAPÍTULO III - GENERACIÓN DE FORMAS
Así, para este caso, el valor de la posición longitudinal del centro de carena que se
obtiene es
21,675
que respecto a la maestra es
44,49
21,675
2 100 1,281 % ó
44,49
Como ya se comentó anteriormente, las características hidrostáticas obtenidas
anteriormente para el calado de proyecto se calcularán para las distintas líneas de agua
para el rango comprendido entre 3 y 5 metros en incrementos de 0,1 metros. De este
modo se obtienen una serie de valores, con los que se trazarán las distintas curvas
hidrostáticas ya explicadas.
A continuación se muestra de forma tabular y gráfica las características
hidrostáticas según los distintos calados estudiados para los diversos trimados obtenidas
mediante Maxsurf.
76
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = -0,5 m
3 702,5 462,618 293,955 0,601 0,519 0,906 0,691 22,925 20,958 1,708 2,75 37,608 4,458 39,314 3,013 5,576
3,1 732,8 472,976 297,286 0,604 0,524 0,909 0,699 22,841 20,837 1,763 2,675 37,16 4,438 38,921 3,047 5,752
3,2 763,5 483,333 300,636 0,607 0,529 0,912 0,706 22,758 20,715 1,818 2,605 36,762 4,423 38,578 3,082 5,933
3,3 794,5 493,654 304,001 0,61 0,534 0,915 0,713 22,676 20,594 1,873 2,539 36,408 4,413 38,279 3,116 6,121
3,4 825,8 503,929 307,449 0,612 0,538 0,918 0,72 22,595 20,47 1,929 2,478 36,121 4,406 38,047 3,151 6,32
3,5 857,5 514,476 311,01 0,613 0,541 0,92 0,725 22,514 20,342 1,984 2,42 35,905 4,404 37,887 3,188 6,532
3,6 889,6 525,823 314,68 0,614 0,544 0,922 0,731 22,433 20,211 2,04 2,366 35,753 4,405 37,79 3,225 6,757
3,7 922,1 533,883 318,202 0,614 0,546 0,924 0,735 22,353 20,091 2,096 2,315 35,562 4,41 37,655 3,262 6,975
3,8 954,9 545,058 322,319 0,614 0,547 0,926 0,74 22,273 19,94 2,152 2,267 35,614 4,418 37,764 3,304 7,247
3,9 988,2 556,128 326,684 0,614 0,549 0,928 0,745 22,192 19,778 2,208 2,222 35,764 4,43 37,97 3,349 7,546
4 1022 567,423 330,247 0,613 0,55 0,93 0,747 22,111 19,665 2,264 2,179 35,63 4,443 37,892 3,385 7,785
4,1 1056 579,077 334,644 0,612 0,55 0,932 0,752 22,03 19,508 2,321 2,138 35,81 4,459 38,129 3,43 8,101
4,2 1090 590,044 338,543 0,611 0,551 0,933 0,755 21,949 19,381 2,378 2,099 35,82 4,477 38,195 3,47 8,382
4,3 1125 601,286 342,538 0,611 0,552 0,935 0,759 21,867 19,252 2,435 2,062 35,871 4,497 38,303 3,511 8,678
4,4 1161 612,755 346,709 0,61 0,553 0,936 0,762 21,785 19,115 2,492 2,026 35,994 4,518 38,484 3,554 8,998
4,5 1196 624,13 350,734 0,609 0,553 0,938 0,766 21,704 18,988 2,55 1,992 36,076 4,541 38,623 3,595 9,313
4,6 1233 635,64 354,84 0,609 0,554 0,939 0,769 21,622 18,859 2,607 1,959 36,195 4,566 38,8 3,637 9,644
4,7 1269 647,308 359,055 0,608 0,554 0,941 0,772 21,54 18,726 2,665 1,926 36,359 4,591 39,021 3,68 9,993
4,8 1306 658,317 362,282 0,611 0,558 0,942 0,778 21,459 18,651 2,723 1,895 36,204 4,617 38,925 3,713 10,256
4,9 1343 668,882 364,748 0,614 0,562 0,943 0,783 21,381 18,62 2,781 1,862 35,827 4,643 38,606 3,739 10,453
5 1381 679,171 366,869 0,617 0,566 0,944 0,787 21,306 18,609 2,838 1,829 35,374 4,668 38,21 3,76 10,622
Tabla 3.4. Características hidrostáticas del casco para un trimado de -0,5 m.
77
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = 0 m
3 707,5 465,367 298,641 0,601 0,541 0,909 0,703 22,482 20,631 1,715 2,788 38,885 4,503 40,6 3,061 5,82
3,1 738,3 475,76 302,078 0,604 0,546 0,912 0,71 22,402 20,502 1,771 2,711 38,44 4,482 40,21 3,096 6,009
3,2 769,4 486,115 305,606 0,606 0,55 0,915 0,716 22,323 20,37 1,827 2,639 38,077 4,466 39,904 3,132 6,209
3,3 801 496,92 309,245 0,607 0,553 0,917 0,722 22,243 20,231 1,883 2,572 37,797 4,455 39,68 3,17 6,423
3,4 832,9 508,402 312,96 0,608 0,555 0,92 0,727 22,163 20,09 1,939 2,509 37,574 4,448 39,512 3,208 6,648
3,5 865,1 516,384 316,619 0,608 0,557 0,922 0,731 22,084 19,953 1,995 2,451 37,352 4,446 39,348 3,245 6,873
3,6 897,8 527,531 320,854 0,608 0,558 0,924 0,736 22,003 19,785 2,052 2,395 37,382 4,447 39,434 3,289 7,15
3,7 930,9 538,485 324,536 0,608 0,559 0,926 0,739 21,922 19,653 2,109 2,344 37,203 4,452 39,312 3,326 7,389
3,8 964,4 550,05 328,85 0,607 0,56 0,928 0,744 21,84 19,486 2,166 2,295 37,29 4,46 39,455 3,371 7,686
3,9 998,3 561,101 332,758 0,607 0,561 0,93 0,747 21,758 19,346 2,223 2,248 37,231 4,471 39,453 3,411 7,956
4 1033 572,547 336,922 0,606 0,561 0,932 0,751 21,675 19,194 2,28 2,204 37,282 4,484 39,562 3,453 8,254
4,1 1067 584,193 341,135 0,606 0,562 0,933 0,754 21,592 19,042 2,338 2,162 37,363 4,5 39,701 3,497 8,566
4,2 1103 595,583 345,034 0,605 0,562 0,935 0,757 21,508 18,91 2,396 2,122 37,341 4,518 39,736 3,537 8,857
4,3 1138 607,633 349,572 0,604 0,563 0,936 0,762 21,425 18,745 2,454 2,083 37,565 4,537 40,019 3,583 9,215
4,4 1174 619,102 353,517 0,604 0,563 0,938 0,765 21,341 18,617 2,512 2,046 37,587 4,557 40,099 3,624 9,528
4,5 1211 630,564 357,259 0,606 0,566 0,939 0,77 21,257 18,502 2,57 2,009 37,55 4,58 40,12 3,662 9,83
4,6 1247 641,265 359,846 0,609 0,57 0,941 0,775 21,175 18,454 2,628 1,973 37,14 4,601 39,768 3,688 10,029
4,7 1284 651,721 362,02 0,613 0,574 0,942 0,779 21,096 18,43 2,687 1,935 36,629 4,622 39,316 3,711 10,196
4,8 1322 661,898 363,934 0,616 0,578 0,943 0,783 21,021 18,423 2,745 1,898 36,073 4,643 38,818 3,73 10,344
4,9 1359 671,634 365,741 0,619 0,582 0,944 0,786 20,95 18,422 2,803 1,861 35,523 4,664 38,325 3,749 10,486
5 1397 681,122 367,377 0,622 0,585 0,945 0,789 20,882 18,432 2,86 1,825 34,955 4,686 37,815 3,766 10,615
Tabla 3.5. Características hidrostáticas del casco para un trimado de 0 m.
78
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = +0,5 m
3 713,7 468,362 304 0,596 0,508 0,906 0,713 22,021 20,278 1,729 2,826 40,414 4,555 42,14 3,116 6,118
3,1 745 479,447 307,715 0,598 0,512 0,909 0,719 21,945 20,128 1,785 2,747 40,029 4,533 41,812 3,154 6,332
3,2 776,8 487,964 311,366 0,599 0,515 0,912 0,723 21,868 19,986 1,842 2,674 39,657 4,516 41,497 3,192 6,547
3,3 808,9 499,041 315,395 0,599 0,518 0,914 0,728 21,79 19,824 1,899 2,606 39,492 4,505 41,388 3,233 6,797
3,4 841,5 510,21 319,793 0,599 0,52 0,917 0,733 21,71 19,642 1,956 2,542 39,504 4,498 41,458 3,278 7,084
3,5 874,4 521,453 323,513 0,599 0,522 0,919 0,736 21,63 19,501 2,014 2,482 39,249 4,496 41,26 3,316 7,323
3,6 907,8 533,213 327,96 0,599 0,524 0,921 0,741 21,548 19,319 2,071 2,426 39,314 4,497 41,383 3,362 7,627
3,7 941,6 544,381 331,987 0,599 0,526 0,923 0,744 21,466 19,164 2,129 2,373 39,225 4,502 41,352 3,403 7,904
3,8 975,9 555,759 336,034 0,598 0,527 0,925 0,747 21,383 19,009 2,187 2,322 39,158 4,51 41,343 3,444 8,19
3,9 1011 567,404 340,285 0,598 0,529 0,927 0,751 21,299 18,845 2,246 2,275 39,19 4,52 41,434 3,488 8,501
4 1046 579,02 344,443 0,597 0,53 0,929 0,754 21,214 18,688 2,304 2,229 39,203 4,533 41,504 3,531 8,813
4,1 1081 590,771 348,651 0,597 0,531 0,931 0,758 21,129 18,531 2,363 2,185 39,251 4,548 41,611 3,574 9,138
4,2 1117 602,402 352,634 0,598 0,533 0,932 0,762 21,043 18,39 2,422 2,143 39,23 4,564 41,649 3,614 9,451
4,3 1153 613,324 355,41 0,602 0,538 0,934 0,767 20,958 18,32 2,481 2,101 38,79 4,581 41,268 3,643 9,659
4,4 1190 623,946 357,648 0,605 0,543 0,935 0,772 20,877 18,281 2,539 2,058 38,204 4,597 40,741 3,666 9,824
4,5 1227 634,317 359,584 0,609 0,548 0,937 0,776 20,799 18,262 2,598 2,015 37,562 4,613 40,158 3,686 9,966
4,6 1264 644,259 361,377 0,613 0,552 0,938 0,779 20,725 18,252 2,656 1,974 36,922 4,63 39,576 3,704 10,101
4,7 1301 653,798 362,986 0,616 0,556 0,939 0,782 20,655 18,253 2,715 1,933 36,266 4,647 38,978 3,721 10,223
4,8 1338 663,28 364,496 0,619 0,561 0,941 0,785 20,588 18,262 2,772 1,892 35,621 4,665 38,392 3,736 10,339
4,9 1376 672,724 365,937 0,623 0,565 0,942 0,788 20,526 18,276 2,83 1,853 34,994 4,684 37,822 3,751 10,451
5 1413 682,141 367,323 0,626 0,569 0,943 0,79 20,466 18,294 2,888 1,816 34,387 4,703 37,272 3,765 10,562
Tabla 3.6. Características hidrostáticas del casco para un trimado de +0,5 m.
79
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = +1 m
3 721,3 470,357 310,582 0,586 0,475 0,907 0,72 21,538 19,919 1,751 2,864 42,535 4,614 44,275 3,183 6,529
3,1 753,3 481,742 314,602 0,588 0,48 0,91 0,725 21,465 19,71 1,808 2,784 42,14 4,592 43,937 3,225 6,761
3,2 785,8 492,892 318,458 0,588 0,483 0,913 0,729 21,388 19,541 1,866 2,71 41,759 4,575 43,615 3,264 6,996
3,3 818,7 504,557 322,829 0,589 0,486 0,916 0,733 21,31 19,353 1,924 2,64 41,667 4,564 43,58 3,309 7,282
3,4 852 515,821 326,909 0,589 0,489 0,918 0,737 21,231 19,187 1,982 2,575 41,474 4,557 43,446 3,351 7,552
3,5 885,7 527,445 331,218 0,589 0,492 0,92 0,74 21,15 19,01 2,041 2,514 41,408 4,554 43,439 3,395 7,849
3,6 919,9 539,173 335,566 0,589 0,494 0,922 0,744 21,067 18,833 2,1 2,456 41,381 4,556 43,47 3,44 8,158
3,7 954,5 550,788 339,728 0,589 0,496 0,924 0,748 20,983 18,668 2,159 2,401 41,298 4,56 43,447 3,482 8,459
3,8 989,6 562,609 344,021 0,589 0,498 0,926 0,751 20,898 18,497 2,218 2,349 41,291 4,567 43,498 3,526 8,781
3,9 1025 574,495 348,302 0,589 0,5 0,928 0,755 20,812 18,328 2,278 2,299 41,298 4,576 43,565 3,57 9,111
4 1061 585,641 351,476 0,592 0,506 0,93 0,761 20,726 18,224 2,337 2,25 40,892 4,587 43,219 3,603 9,347
4,1 1097 596,272 353,811 0,597 0,511 0,932 0,765 20,643 18,168 2,397 2,201 40,214 4,597 42,601 3,627 9,512
4,2 1133 606,58 355,771 0,601 0,516 0,933 0,769 20,563 18,134 2,456 2,152 39,458 4,607 41,904 3,647 9,649
4,3 1170 616,67 357,549 0,605 0,522 0,935 0,773 20,487 18,113 2,515 2,104 38,696 4,618 41,201 3,665 9,775
4,4 1207 626,217 359,128 0,608 0,527 0,936 0,776 20,415 18,104 2,573 2,057 37,923 4,63 40,487 3,681 9,888
4,5 1244 635,701 360,588 0,612 0,531 0,938 0,779 20,347 18,104 2,632 2,011 37,163 4,642 39,786 3,696 9,994
4,6 1281 645,246 361,972 0,616 0,536 0,939 0,781 20,282 18,109 2,69 1,966 36,428 4,655 39,109 3,71 10,096
4,7 1318 654,765 363,293 0,619 0,541 0,94 0,784 20,222 18,12 2,748 1,923 35,719 4,67 38,458 3,724 10,195
4,8 1355 664,238 364,578 0,622 0,545 0,941 0,786 20,165 18,136 2,805 1,881 35,043 4,686 37,839 3,737 10,295
4,9 1392 673,646 365,825 0,626 0,549 0,943 0,788 20,111 18,155 2,863 1,841 34,394 4,704 37,248 3,75 10,394
5 1430 683,04 367,041 0,629 0,553 0,944 0,791 20,06 18,177 2,92 1,803 33,773 4,722 36,684 3,762 10,493
Tabla 3.7. Características hidrostáticas del casco para un trimado de +1 m.
80
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = +1,5 m
3 730,4 475,082 320,325 0,568 0,44 0,909 0,721 21,018 19,593 1,781 2,904 46,567 4,684 48,321 3,283 7,274
3,1 763,4 487,31 324,068 0,572 0,446 0,912 0,727 20,951 19,324 1,84 2,822 45,75 4,661 47,564 3,322 7,472
3,2 796,8 499,177 327,441 0,576 0,453 0,914 0,734 20,878 19,074 1,899 2,746 44,802 4,643 46,675 3,356 7,639
3,3 830,5 510,856 331,203 0,578 0,457 0,917 0,738 20,801 18,864 1,958 2,675 44,228 4,632 46,161 3,395 7,866
3,4 864,7 522,78 335,551 0,579 0,461 0,919 0,742 20,721 18,663 2,018 2,608 44,049 4,625 46,042 3,439 8,166
3,5 899,3 534,707 339,92 0,579 0,464 0,921 0,745 20,638 18,475 2,078 2,545 43,936 4,621 45,989 3,484 8,481
3,6 934,3 546,517 344,14 0,579 0,467 0,924 0,749 20,554 18,301 2,138 2,485 43,795 4,621 45,908 3,527 8,794
3,7 969,8 558,082 347,868 0,582 0,472 0,926 0,754 20,469 18,158 2,198 2,427 43,478 4,624 45,651 3,566 9,072
3,8 1006 568,827 350,382 0,587 0,478 0,927 0,759 20,386 18,086 2,258 2,369 42,712 4,627 44,946 3,591 9,246
3,9 1042 579,268 352,44 0,591 0,484 0,929 0,764 20,305 18,039 2,318 2,312 41,836 4,629 44,131 3,613 9,386
4 1078 589,471 354,15 0,596 0,49 0,931 0,767 20,229 18,013 2,378 2,256 40,902 4,632 43,257 3,63 9,499
4,1 1114 599,101 355,824 0,6 0,496 0,933 0,77 20,157 17,988 2,437 2,201 40,026 4,637 42,44 3,647 9,615
4,2 1151 608,57 357,258 0,604 0,501 0,934 0,773 20,088 17,978 2,496 2,148 39,128 4,643 41,602 3,662 9,712
4,3 1187 617,989 358,588 0,608 0,507 0,936 0,776 20,023 17,975 2,555 2,096 38,256 4,65 40,789 3,676 9,804
4,4 1224 627,378 359,848 0,612 0,512 0,937 0,778 19,962 17,977 2,613 2,047 37,418 4,659 40,01 3,688 9,893
4,5 1261 636,746 361,051 0,615 0,517 0,938 0,78 19,904 17,984 2,671 1,999 36,616 4,669 39,266 3,701 9,979
4,6 1298 646,126 362,22 0,619 0,522 0,94 0,783 19,85 17,996 2,729 1,953 35,853 4,681 38,561 3,713 10,066
4,7 1335 655,501 363,36 0,622 0,526 0,941 0,785 19,799 18,011 2,786 1,909 35,125 4,694 37,891 3,724 10,153
4,8 1372 664,939 364,473 0,625 0,531 0,942 0,787 19,751 18,029 2,843 1,867 34,433 4,709 37,256 3,736 10,24
4,9 1410 674,378 365,566 0,628 0,535 0,943 0,789 19,706 18,05 2,9 1,827 33,773 4,726 36,654 3,747 10,327
5 1447 683,855 366,652 0,631 0,539 0,944 0,791 19,664 18,075 2,957 1,788 33,147 4,744 36,085 3,758 10,416
Tabla 3.8. Características hidrostáticas del casco para un trimado de +1,5 m.
81
CAPÍTULO III - GENERACIÓN DE FORMAS
5.3.1 Curvas hidrostáticas para los distintos trimados.
Trimado = -0,5 m
Desplazamiento
5,2
5
4,8
4,6
4,4
Calado (m)
4,2
4
3,8
3,6
3,4
3,2
3
2,8
600 700 800 900 1000 1100 1200 1300 1400 1500
Desplazamiento (t)
Gráfica 3.2. Curva del desplazamiento para un trimado de -0,5 m.
Superficie mojada Área de la flotación
5,2
5
4,8
4,6
4,4
Calado (m)
4,2
4
3,8
3,6
3,4
3,2
3
2,8
200 300 400 500 600 700 800
Área (m2)
Gráfica 3.3. Curvas de la superficie mojada y del área de la flotación para un trimado de -0,5 m.
82
CAPÍTULO III - GENERACIÓN DE FORMAS
CB CP CM CF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
0,5 0,55 0,6 0,65 0,7 0,75 0,8 0,85 0,9 0,95 1
Coeficiente
Gráfica 3.4. Curvas de los coeficientes de formas para un trimado de -0,5 m.
LCB LCF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
18 18,5 19 19,5 20 20,5 21 21,5 22 22,5 23
Centro longitudinal desde Ppp (m)
Gráfica 3.5. Curvas de la posición longitudinal del centro de carena y del centro de la flotación para un
trimado de -0,5 m.
83
CAPÍTULO III - GENERACIÓN DE FORMAS
KB
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,7 1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9
KB (m)
Gráfica 3.6. Curva de la altura del centro de carena para un trimado de -0,5 m.
BMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8
BMt (m)
Gráfica 3.7. Curva del radio metacéntrico transversal para un trimado de -0,5 m.
84
CAPÍTULO III - GENERACIÓN DE FORMAS
BMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
34,8 35,2 35,6 36 36,4 36,8 37,2 37,6 38
BMl (m)
Gráfica 3.8. Curva del radio metacéntrico longitudinal para un trimado de -0,5 m.
KMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
4,38 4,4 4,42 4,44 4,46 4,48 4,5 4,52 4,54 4,56 4,58 4,6 4,62 4,64 4,66 4,68 4,7
KMt (m)
Gráfica 3.9. Curva de la distancia del metacentro transversal sobre la base para un trimado de -0,5 m.
85
CAPÍTULO III - GENERACIÓN DE FORMAS
KMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
37 37,4 37,8 38,2 38,6 39 39,4 39,8
KMl (m)
Gráfica 3.10. Curva de la distancia del metacentro longitudinal sobre la base para un trimado de -0,5 m.
TCI
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
3 3,1 3,2 3,3 3,4 3,5 3,6 3,7 3,8
TCI (t/cm)
Gráfica 3.11. Curva de las toneladas por centímetro de inmersión para un trimado de -0,5 m.
86
CAPÍTULO III - GENERACIÓN DE FORMAS
MTC
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
5,5 6 6,5 7 7,5 8 8,5 9 9,5 10 10,5 11
MTC (t∙m)
Gráfica 3.12. Curva del momento para alterar el trimado un centímetro para un trimado de -0,5 m.
87
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = 0 m
Desplazamiento
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
600 700 800 900 1000 1100 1200 1300 1400 1500
Desplazamiento (t)
Gráfica 3.13. Curva del desplazamiento para un trimado de 0 m.
Superficie mojada Área de la flotación
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
200 300 400 500 600 700 800
Área (m2)
Gráfica 3.14. Curvas de la superficie mojada y del área de la flotación para un trimado de 0 m.
88
CAPÍTULO III - GENERACIÓN DE FORMAS
CB CP CM CF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
0,5 0,55 0,6 0,65 0,7 0,75 0,8 0,85 0,9 0,95 1
Coeficiente
Gráfica 3.15. Curvas de los coeficientes de formas para un trimado de 0 m.
LCB LCF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
18 18,5 19 19,5 20 20,5 21 21,5 22 22,5 23
Centro longitudinal desde Ppp (m)
Gráfica 3.16. Curvas de la posición longitudinal del centro de carena y del centro de la flotación para un
trimado de 0 m.
89
CAPÍTULO III - GENERACIÓN DE FORMAS
KB
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,7 1,8 1,9 2,0 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9
KB (m)
Gráfica 3.17. Curva de la altura del centro de carena para un trimado de 0 m.
BMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8
BMt (m)
Gráfica 3.18. Curva del radio metacéntrico transversal para un trimado de 0 m.
90
CAPÍTULO III - GENERACIÓN DE FORMAS
BMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
34,8 35,2 35,6 36 36,4 36,8 37,2 37,6 38 38,4 38,8 39,2
BMl (m)
Gráfica 3.19. Curva del radio metacéntrico longitudinal para un trimado de 0 m.
KMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
4,42 4,44 4,46 4,48 4,5 4,52 4,54 4,56 4,58 4,6 4,62 4,64 4,66 4,68 4,7
KMt (m)
Gráfica 3.20. Curva de la distancia del metacentro transversal sobre la base para un trimado de 0 m.
91
CAPÍTULO III - GENERACIÓN DE FORMAS
KMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
37,8 38,1 38,4 38,7 39 39,3 39,6 39,9 40,2 40,5 40,8
KMl (m)
Gráfica 3.21. Curva de la distancia del metacentro longitudinal sobre la base para un trimado de 0 m.
TCI
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
3 3,1 3,2 3,3 3,4 3,5 3,6 3,7 3,8
TCI (t/cm)
Gráfica 3.22. Curva de las toneladas por centímetro de inmersión para un trimado de 0 m.
92
CAPÍTULO III - GENERACIÓN DE FORMAS
MTC
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
5,5 6 6,5 7 7,5 8 8,5 9 9,5 10 10,5 11
MTC (t∙m)
Gráfica 3.23. Curva del momento para alterar el trimado un centímetro para un trimado de 0 m.
93
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = +0,5 m
Desplazamiento
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
600 700 800 900 1000 1100 1200 1300 1400 1500
Desplazamiento (t)
Gráfica 3.24. Curva del desplazamiento para un trimado de +0,5 m.
Superficie mojada Área de la flotación
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
200 300 400 500 600 700 800
Área (m2)
Gráfica 3.25. Curvas de la superficie mojada y del área de la flotación para un trimado de +0,5 m.
94
CAPÍTULO III - GENERACIÓN DE FORMAS
CB CP CM CF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
0,45 0,5 0,55 0,6 0,65 0,7 0,75 0,8 0,85 0,9 0,95 1
Coeficiente
Gráfica 3.26. Curvas de los coeficientes de formas para un trimado de +0,5 m.
LCB LCF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
18 18,5 19 19,5 20 20,5 21 21,5 22 22,5
Centro longitudinal desde Ppp (m)
Gráfica 3.27. Curvas de la posición longitudinal del centro de carena y del centro de la flotación para un
trimado de +0,5 m.
95
CAPÍTULO III - GENERACIÓN DE FORMAS
KB
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,7 1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9
KB (m)
Gráfica 3.28. Curva de la altura del centro de carena para un trimado de +0,5 m.
BMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9
BMt (m)
Gráfica 3.29. Curva del radio metacéntrico transversal para un trimado de +0,5 m.
96
CAPÍTULO III - GENERACIÓN DE FORMAS
BMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
34,2 34,8 35,4 36 36,6 37,2 37,8 38,4 39 39,6 40,2 40,8
BMl (m)
Gráfica 3.30. Curva del radio metacéntrico longitudinal para un trimado de +0,5 m.
KMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
4,48 4,5 4,52 4,54 4,56 4,58 4,6 4,62 4,64 4,66 4,68 4,7 4,72
KMt (m)
Gráfica 3.31. Curva de la distancia del metacentro transversal sobre la base para un trimado de +0,5 m.
97
CAPÍTULO III - GENERACIÓN DE FORMAS
KMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
36,6 37,2 37,8 38,4 39 39,6 40,2 40,8 41,4 42 42,6
KMl (m)
Gráfica 3.32. Curva de la distancia del metacentro longitudinal sobre la base para un trimado de +0,5 m.
TCI
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
3 3,1 3,2 3,3 3,4 3,5 3,6 3,7 3,8
TCI (t/cm)
Gráfica 3.33. Curva de las toneladas por centímetro de inmersión para un trimado de +0,5 m.
98
CAPÍTULO III - GENERACIÓN DE FORMAS
MTC
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
5,5 6 6,5 7 7,5 8 8,5 9 9,5 10 10,5 11
MTC (t∙m)
Gráfica 3.34. Curva del momento para alterar el trimado un centímetro para un trimado de +0,5 m.
99
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = +1 m
Desplazamiento
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
600 700 800 900 1000 1100 1200 1300 1400 1500
Desplazamiento (t)
Gráfica 3.35. Curva del desplazamiento para un trimado de +1 m.
Superficie mojada Área de la flotación
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
200 300 400 500 600 700 800
Área (m2)
Gráfica 3.36. Curvas de la superficie mojada y del área de la flotación para un trimado de +1 m.
100
CAPÍTULO III - GENERACIÓN DE FORMAS
CB CP CM CF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
0,45 0,5 0,55 0,6 0,65 0,7 0,75 0,8 0,85 0,9 0,95 1
Coeficiente
Gráfica 3.37. Curvas de los coeficientes de formas para un trimado de +1 m.
LCB LCF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
17,5 18 18,5 19 19,5 20 20,5 21 21,5 22
Centro longitudinal desde Ppp (m)
Gráfica 3.38. Curvas de la posición longitudinal del centro de carena y del centro de la flotación para un
trimado de +1 m.
101
CAPÍTULO III - GENERACIÓN DE FORMAS
KB
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,7 1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9 3
KB (m)
Gráfica 3.39. Curva de la altura del centro de carena para un trimado de +1 m.
BMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,7 1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9 3
BMt (m)
Gráfica 3.40. Curva del radio metacéntrico transversal para un trimado de +1 m.
102
CAPÍTULO III - GENERACIÓN DE FORMAS
BMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
33,2 33,8 34,4 35 35,6 36,2 36,8 37,4 38 38,6 39,2 39,8 40,4 41 41,6 42,2 42,8 43,4
BMl (m)
Gráfica 3.41. Curva del radio metacéntrico longitudinal para un trimado de +1 m.
KMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
4,54 4,56 4,58 4,6 4,62 4,64 4,66 4,68 4,7 4,72 4,74
KMt (m)
Gráfica 3.42. Curva de la distancia del metacentro transversal sobre la base para un trimado de +1 m.
103
CAPÍTULO III - GENERACIÓN DE FORMAS
KMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
36,5 37 37,5 38 38,5 39 39,5 40 40,5 41 41,5 42 42,5 43 43,5 44 44,5 45
KMl (m)
Gráfica 3.43. Curva de la distancia del metacentro longitudinal sobre la base para un trimado de +1 m.
TCI
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
3,1 3,2 3,3 3,4 3,5 3,6 3,7 3,8
TCI (t/cm)
Gráfica 3.44. Curva de las toneladas por centímetro de inmersión para un trimado de +1 m.
104
CAPÍTULO III - GENERACIÓN DE FORMAS
MTC
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
6 6,5 7 7,5 8 8,5 9 9,5 10 10,5 11
MTC (t∙m)
Gráfica 3.45. Curva del momento para alterar el trimado un centímetro para un trimado de +1 m.
105
CAPÍTULO III - GENERACIÓN DE FORMAS
Trimado = +1,5 m
Desplazamiento
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
600 700 800 900 1000 1100 1200 1300 1400 1500
Desplazamiento (t)
Gráfica 3.46. Curva del desplazamiento para un trimado de +1,5 m.
Superficie mojada Área de la flotación
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
200 300 400 500 600 700 800
Área (m2)
Gráfica 3.47. Curvas de la superficie mojada y del área de la flotación para un trimado de +1,5 m.
106
CAPÍTULO III - GENERACIÓN DE FORMAS
CB CP CM CF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
0,4 0,45 0,5 0,55 0,6 0,65 0,7 0,75 0,8 0,85 0,9 0,95 1
Coeficiente
Gráfica 3.48. Curvas de los coeficientes de formas para un trimado de +1,5 m.
LCB LCF
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
17,5 18 18,5 19 19,5 20 20,5 21 21,5
Centro longitudinal desde Ppp (m)
Gráfica 3.49. Curvas de la posición longitudinal del centro de carena y del centro de la flotación para un
trimado de +1,5 m.
107
CAPÍTULO III - GENERACIÓN DE FORMAS
KB
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,7 1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9 3
KB (m)
Gráfica 3.50. Curva de la altura del centro de carena para un trimado de +1,5 m.
BMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
1,7 1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9 3
BMt (m)
Gráfica 3.51. Curva del radio metacéntrico transversal para un trimado de +1,5 m.
108
CAPÍTULO III - GENERACIÓN DE FORMAS
BMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
32,8 33,8 34,8 35,8 36,8 37,8 38,8 39,8 40,8 41,8 42,8 43,8 44,8 45,8 46,8
BMl (m)
Gráfica 3.52. Curva del radio metacéntrico longitudinal para un trimado de +1,5 m.
KMt
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
4,6 4,62 4,64 4,66 4,68 4,7 4,72 4,74 4,76
KMt (m)
Gráfica 3.53. Curva de la distancia del metacentro transversal sobre la base para un trimado de +1,5 m.
109
CAPÍTULO III - GENERACIÓN DE FORMAS
KMl
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
35,7 36,7 37,7 38,7 39,7 40,7 41,7 42,7 43,7 44,7 45,7 46,7 47,7 48,7
KMl (m)
Gráfica 3.54. Curva de la distancia del metacentro longitudinal sobre la base para un trimado de +1,5 m.
TCI
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
3,2 3,3 3,4 3,5 3,6 3,7 3,8
TCI (t/cm)
Gráfica 3.55. Curva de las toneladas por centímetro de inmersión para un trimado de +1,5 m.
110
CAPÍTULO III - GENERACIÓN DE FORMAS
MTC
5,2
5
4,8
4,6
4,4
4,2
Calado (m)
4
3,8
3,6
3,4
3,2
3
2,8
6,5 7 7,5 8 8,5 9 9,5 10 10,5 11
MTC (t∙m)
Gráfica 3.56. Curva del momento para alterar el trimado un centímetro para un trimado de +1,5 m.
111
CAPÍTULO IV
1. INTRODUCCIÓN
113
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
114
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
60000
50000
40000
Resistencia (kg)
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 4.1. Curvas de las componentes de la resistencia obtenidas mediante el método de J. Holtrop y G. G. J. Mennen.
115
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
Rt
60000
50000
40000
Resistencia (kg)
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 4.2. Curva de resistencia obtenida mediante el método de J. Holtrop y G. G. J. Mennen.
116
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
Como se comentó en el Capítulo I, la magnitud relativa de cada una de las
componentes de la resistencia al avance de un barco depende de la velocidad a la que se
navega y de las formas de la carena.
Para valores pequeños de la velocidad la resistencia ocasionada por la viscosidad
predomina sobre la resistencia por formación de olas, llegando la resistencia viscosa a
ser aproximadamente un 81,5% de la resistencia total. A medida que la velocidad crece
la resistencia por formación de olas también lo hace, mientras que la que tiene su origen
en la viscosidad lo hace en menor proporción. Cuando la velocidad llega a ser
suficientemente alta, el valor que toma la resistencia debida a la formación de olas es
tan alto que llega a ser un 84% respecto de la resistencia total.
117
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
Rv Rw
90
80
70
60
Resistencia (%RT)
50
40
30
20
10
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 4.3. Curvas de la resistencia viscosa y por formación de olas respecto a la resistencia total.
118
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
V V Fr Re CF RR/RT CT RT
(nudos) (m/s) (kg)
3 1,5434 0,0739 5,834E+07 2,256E-03 0,0187 3,313E-03 236,0
4 2,0578 0,0985 7,778E+07 2,161E-03 0,0728 3,404E-03 431,2
5 2,5723 0,1231 9,723E+07 2,092E-03 0,1272 3,537E-03 700,0
6 3,0867 0,1478 1,167E+08 2,038E-03 0,1820 3,707E-03 1056,5
7 3,6012 0,1724 1,361E+08 1,993E-03 0,2370 3,917E-03 1519,3
8 4,1156 0,1970 1,556E+08 1,956E-03 0,2924 4,171E-03 2113,0
9 4,6301 0,2216 1,750E+08 1,924E-03 0,3481 4,478E-03 2871,4
10 5,1445 0,2463 1,945E+08 1,896E-03 0,4041 4,852E-03 3841,1
11 5,6590 0,2709 2,139E+08 1,872E-03 0,4605 5,313E-03 5089,2
12 6,1734 0,2955 2,333E+08 1,850E-03 0,5171 5,891E-03 6715,1
13 6,6879 0,3201 2,528E+08 1,829E-03 0,5741 6,632E-03 8872,4
14 7,2023 0,3448 2,722E+08 1,811E-03 0,6314 7,613E-03 11812,9
15 7,7168 0,3694 2,917E+08 1,794E-03 0,6890 8,971E-03 15978,0
Tabla 4.4. Valores de la resistencia del casco sin bulbo obtenidos mediante el método de Amadeo García
Gómez.
119
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
120
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
Rt (CB) Rt (SB)
20000
18000
16000
14000
12000
Resistencia (kg)
10000
8000
6000
4000
2000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Velocidad (nudos)
Gráfica 4.4. Curvas de resistencia del casco con bulbo y sin bulbo obtenida mediante el método de Amadeo García Gómez.
121
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
V V Fr Re CF CF+CA RF RR/∆ RR RA RT
(nudos) (m/s) (kg) (kg) (kg) (kg)
3 1,54335 0,0737 5,808E+07 2,257E-03 2,757E-03 160,831 3,114E-07 0,322 35,623 196,8
4 2,0578 0,0983 7,744E+07 2,163E-03 2,663E-03 273,919 1,497E-05 15,469 63,330 352,7
5 2,57225 0,1229 9,681E+07 2,093E-03 2,593E-03 414,252 8,992E-05 92,890 98,954 606,1
6 3,0867 0,1474 1,162E+08 2,039E-03 2,539E-03 581,049 2,381E-04 245,961 142,493 969,5
7 3,60115 0,1720 1,355E+08 1,995E-03 2,495E-03 773,698 4,284E-04 442,545 193,949 1410,2
8 4,1156 0,1966 1,549E+08 1,957E-03 2,457E-03 991,698 6,305E-04 651,344 253,322 1896,4
9 4,63005 0,2211 1,743E+08 1,925E-03 2,425E-03 1234,628 8,469E-04 874,840 320,610 2430,1
10 5,1445 0,2457 1,936E+08 1,898E-03 2,398E-03 1502,126 1,197E-03 1236,839 395,815 3134,8
11 5,65895 0,2703 2,130E+08 1,873E-03 2,373E-03 1793,873 1,779E-03 1837,897 478,936 4110,7
12 6,1734 0,2949 2,323E+08 1,851E-03 2,351E-03 2109,587 3,416E-03 3528,690 569,974 6208,3
13 6,68785 0,3194 2,517E+08 1,831E-03 2,331E-03 2449,016 5,074E-03 5240,995 668,928 8358,9
14 7,2023 0,3440 2,711E+08 1,812E-03 2,312E-03 2811,931 7,162E-03 7398,428 775,798 10986,2
15 7,71675 0,3686 2,904E+08 1,796E-03 2,296E-03 3198,121 1,169E-02 12077,984 890,584 16166,7
16 8,2312 0,3931 3,098E+08 1,780E-03 2,280E-03 3607,394 1,898E-02 19608,956 1013,287 24229,6
17 8,74565 0,4177 3,291E+08 1,766E-03 2,266E-03 4039,572 2,762E-02 28531,308 1143,906 33714,8
18 9,2601 0,4423 3,485E+08 1,752E-03 2,252E-03 4494,491 3,604E-02 37228,551 1282,441 43005,5
19 9,77455 0,4669 3,679E+08 1,740E-03 2,240E-03 4971,996 4,340E-02 44828,815 1428,892 51229,7
20 10,289 0,4914 3,872E+08 1,728E-03 2,228E-03 5471,943 4,957E-02 51208,791 1583,260 58264,0
Tabla 4.6. Valores de la resistencia obtenidos mediante el método de Van Oortmersen.
122
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
Rt
60000
50000
40000
Resistencia (kg)
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 4.5. Curva de resistencia obtenida mediante el método de Van Oortmersen.
123
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
A continuación se presentan los resultados del método de Sander M. Calisal y Dan McGreer.
124
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
125
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
Rt (PC) Rt (CP)
90000
80000
70000
60000
Resistencia (kg)
50000
40000
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 4.6. Curvas de resistencia a plena carga y a carga parcial obtenidas mediante el método de Sander M. Calisal y Dan McGreer.
126
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
6. COMPARACIÓN DE RESULTADOS
Velocidad Van Oortmersen Holtrop Calisal (CP) Calisal (PC) Amadeo (CB) Amadeo (SB)
V RT RT RT RT RT RT
(nudos) (kg) (kg) (kg) (kg) (kg) (kg)
3 196,8 252,3 - - 239,3 236,0
4 352,7 433,2 183,8 - 427,7 431,2
5 606,1 659,5 428,9 - 688,8 700,0
6 969,5 931,8 724,4 - 1048,4 1056,5
7 1410,2 1258,5 1081,9 - 1543,5 1519,3
8 1896,4 1665,5 1519,0 1424,6 2212,9 2113,0
9 2430,1 2208,5 2062,3 2693,5 3081,8 2871,4
10 3134,8 2961,6 2752,3 4112,3 4155,3 3841,1
11 4110,7 4123,9 3652,2 5743,2 5442,7 5089,2
12 6208,3 5816,1 4862,3 7682,5 7008,3 6715,1
13 8358,9 7572,7 6542,3 10081,1 9054,6 8872,4
14 10986,2 9658,1 8946,9 13176,0 12156,2 11812,9
15 16166,7 13132,3 12478,7 17337,5 18436,9 15978,0
16 24229,6 19243,0 17765,1 23137,4 - -
17 33714,8 27884,1 25767,3 31446,2 - -
18 43005,5 38150,7 37932,1 43568,0 - -
19 51229,7 48545,5 56398,4 61423,6 - -
20 58264,0 57035,5 84273,2 87796,1 - -
Tabla 4.9. Comparación de los valores de las resistencias obtenidos mediante los distintos métodos de
predicción estadística.
127
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
90000
80000
70000
60000
Resistencia (kg)
50000
40000
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 4.7. Comparación de las curvas de resistencia obtenidas mediante los distintos métodos de predicción estadística.
128
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
Como ya se comentó en el Capítulo II cuando se habló de la elección de los
métodos numéricos, no existe un único método que sea preciso para un amplio rango de
velocidades. Esto se hace visible en el hecho de que no aparezcan los valores de la
resistencia para algunas velocidades en algunos métodos debido a que en la realización
de los cálculos de dichos métodos se obtenían unos valores que no eran coherentes. Esto
se debe a que cada método tiene un rango de aplicación determinado, lo que hace que
para un número de Froude concreto unos se ajusten mejor que otros.
A esto hay que añadir que los métodos de predicción estadística están enfocados a
buques pesqueros y, como ya se comentó en el Capítulo I, un ‘buque pesquero’ abarca
una gama amplísima de tipos y tamaños de embarcaciones, por lo que puede ocurrir
que, aunque el buque en cuestión cumpla con los rangos especificados por cada método,
no se obtengan unos valores razonables.
Una vez hecha esta aclaración, como se puede ver en la gráfica anterior los
valores de la resistencia al avance de todos los métodos son en general bastante
similares en el rango de velocidades estudiado.
El método de Van Oortmersen tiene un rango del número de Froude de 0 – 0,50,
lo que equivale a una velocidad máxima de 20 nudos. Aunque para números de Froude
inferiores a 0,30 (12 nudos) posee el inconveniente de mostrar inestabilidades, viendo
los resultados obtenidos no parece que posea ninguna.
Por otro lado el método de Holtrop tiene un límite máximo del número de Froude
de 0,38, lo que equivale a una velocidad ‘tope’ de unos 16 nudos. Sin embargo se
obtienen unos valores bastante similares que los de Van Oortmersen, por lo que parece
que se ajusta bien para velocidades superiores a 0,38.
Ambos métodos a partir de los 18 nudos presentan un cambio en la curvatura de la
curva resistencia-velocidad.
El método de Amadeo García tiene un rango de aplicación del número de Froude
de 0,25 – 0,40 (hasta 16 nudos). Sin embargo sólo se ha alcanzado en los cálculos una
velocidad de 15 nudos debido a que, a partir de este punto, se obtenían unos resultados
incoherentes. Cabe añadir que tanto para el caso del buque con bulbo y sin él se han
conseguido unos resultados bastante similares.
Finalmente el método de Casisal, como ya se comentó, es un método actual en el
que las formas base son similares a los buques actuales y que ajusta bien para números
de Froude por debajo de 0,30. Prácticamente se obtienen unos resultados muy similares
respecto a los otros métodos empleados. Sin embargo, a partir de aproximadamente 18
nudos se puede apreciar cómo aparecen claras diferencias en comparación con el
método de Calisal y McGreer, llegando a dar valores de unos 30000 kg de diferencia a
una velocidad de 20 nudos.
129
CAPÍTULO IV – APLICACIÓN DE LOS MÉTODOS DE PREDICCIÓN ESTADÍSTICA SELECCIONADOS
7. CONCLUSIONES
En base a los resultados obtenidos, todos ellos alcanzan valores bastante parecidos
entre ellos. Sin embargo, como se ha comentado anteriormente, a partir de unos 18
nudos los resultados discrepan entre ellos.
Por un lado están los obtenidos por Van Oortmersen y Holtrop, los cuales son
bastante similares entre sí; y, por otro, los obtenidos mediante Calisal.
Supuestamente el método de Van Oortmersen tiene un rango mayor del número
de Froude (hasta 0,50, equivalente a 20 nudos) y sus resultados han coincidido
prácticamente con los de Holtrop (cuyo 0,38, equivalente a 16 nudos). Por el
contrario, la curva de resistencia-velocidad a partir de 18 nudos adquiere una pendiente
menor.
Sin embargo, en el método de Calisal, la pendiente de la curva resistencia-
velocidad a partir de los 18 nudos continúa creciendo. Además es un método actual en
el que las formas base son similares a los buques actuales.
Es por esto último que, a falta de resultados experimentales, lo cual sería lo más
fiable de cara a la elección del método más aproximado, se tomará como más correcto el
método de Sander M. Calisal, aunque posteriormente se comparará con los obtenidos
mediante modelos matemáticos (CFD).
130
CAPÍTULO V
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1. INTRODUCCIÓN
132
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
popularidad son los métodos FVM (métodos de los volúmenes finitos) que permiten la
integración de las ecuaciones de masa y momentum para cada celda individual y
aproximar su valor en el nodo central.
El problema de la resistencia, entendida como el efecto que una corriente fluida
ejercía sobre un cuerpo inmerso en su seno fue uno de los dos pilares en donde se
sustentó la génesis y evolución de la mecánica de fluidos en los siglos XVII y XVIII. El
interés que despertó este problema se explica, aunque no exclusivamente, por su
evidente utilidad práctica, ya que barcos, molinos y otras máquinas se encontraban
sometidos a estas fuerzas y cumplían sus funciones por la existencia de las mismas.
Los primeros enfoques para explicar la resistencia, allá por el último tercio del
siglo XVII, siguieron unas pautas muy diferentes y fue a mediados del siglo XVIII
cuando se llegó a las ecuaciones básicas del movimiento, aunque la falta de capacidad
de cálculo matemático y de los conocimientos adecuados de la viscosidad
imposibilitaron el avance en esa dirección.
Entre finales del siglo XVII y principios del siglo XVIII se constituyó la ‘teoría
del impacto’; teoría basada en dos hipótesis: la primera que suponía que el fluido estaba
formado por un conjunto de partículas individualizadas e independientes; la segunda
que asimilaba la resistencia a la suma de los efectos mecánicos de todos los choques.
Esta teoría no daba explicación a las fuerzas en las popas y así Euler a mediados
del siglo XVIII conjeturó que el fluido contorneaba el cuerpo introduciendo el modelo
currentilíneo con lo cual se aproximó a las concepciones actuales.
Por esa misma época Jean le Rond D’Alambert expone su idea sobre el
tratamiento de la estructura del movimiento de una corriente fluida alrededor de un
cuerpo. El fluido, al que supone animado con un movimiento uniforme aguas arriba, al
acercarse al cuerpo empieza a desviarse formando unas líneas de corriente que lo
contornearán para luego volver aguas abajo a la condición de uniforme. El cambio es
continuo y gradual, y demuestra que la estructura del campo fluido, esto es, la forma
geométrica de las líneas de corriente es independiente de la velocidad, y sólo depende
de la forma del cuerpo.
Como consecuencia la velocidad en un punto concreto será proporcional a la
existente aguas arriba, hecho que permite adimensionalizar el movimiento. Por otro
lado, el efecto del fluido sobre el cuerpo será el producido por la presión y no el choque.
Para obtener el campo fluido parte de dos premisas: la continuidad y la dinámica. Así
hasta la llegada de los ordenadores su resolución quedó paralizada debido a su
dificultad.
En la actualidad son diversas las parcelas hidrodinámicas que dan cabida a los
cálculos mediante CFD, si bien es cierto que el campo de la aeronáutica fue el pionero
en este tipo de cálculos.
La razón de lo anterior la debemos encontrar en el tipo de fluido, aire en lugar de
agua, y en la no existencia de la interface aire-agua con la formación de olas que ello
conlleva.
133
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
La diferencia de fluido hace resaltar la presencia de la viscosidad y de ahí que los
aspectos hidrodinámicos que mejor pueden ser abordados por la CFD sean aquellos
problemas en donde los efectos viscosos sean pequeños, como pueden ser la predicción
de los movimientos del buque, la componente de la resistencia debida a la formación de
olas y la acción del propulsor.
En la hidrodinámica aparece en los siguientes campos:
134
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.3 Ventajas e inconvenientes de la CFD
El análisis con técnicas CFD reduce de forma sustancial en tiempo y costes en los
nuevos diseños, posibilita analizar sistemas y condiciones muy difíciles de simular
experimentalmente y es capaz de estudiar sistemas bajo condiciones peligrosas o más
allá de sus condiciones límites de funcionamiento (accidentes con sustancias tóxicas).
Por otro lado, el nivel de detalle prácticamente ilimitado. Los métodos
experimentales son tanto más caros cuanto mayor es el número de puntos de medida,
mientras que los programas CFD pueden generar gran cantidad de información sin coste
añadido y con posibilidad de hacer estudios paramétricos.
Además añade valor al producto, es decir, existe la capacidad de generación de
gráficos que permiten una mejor comprensión del resultado y así estimular la compra
del producto.
Sin embargo, no todo pueden ser ventajas en este tipo de herramientas, uno de los
principales inconvenientes de las técnicas CFD, es el gran conocimiento de las
ecuaciones que modelan ciertos fenómenos físicos, necesitando personal con grandes
conocimientos en la materia.
Otra desventaja es que no siempre es posible llegar a resultados lo suficientemente
precisos, dando lugar a grandes errores en cuestiones básicas. También es necesaria la
simplificación del fenómeno a estudiar para que el hardware y el software puedan
abordarlo. El resultado será tanto más preciso cuanto más adecuadas hayan sido las
hipótesis y simplificaciones realizadas.
Además está la existencia de insuficientes e incompletos modelos para simular el
efecto de la turbulencia, fenómenos multifásicos o la combustión, entre otros.
Por último y no menos importante, la tendencia humana a creerse todo lo
procedente del ordenador, sobre todo cuando se presentan resultados de forma atractiva.
Con respecto a la necesidad de la simplificación de los fenómenos también se
debe tener en cuenta la aproximación a la realidad. Se trata de las distintas
aproximaciones que se realizan y se llevan a cabo en ingeniería para la resolución
numérica de las ecuaciones de gobierno que modelan un determinado fenómeno.
Conocidas las ecuaciones de gobierno del fenómeno físico a estudiar será
necesario introducir simplificaciones al problema que ayuden a obtener un resultado
lo más preciso posible pero con un coste y dificultad tan reducido como sea posible.
Así pues, como de antemano no se conoce la solución analítica de las ecuaciones
resulta muy difícil conocer cuáles de las hipótesis realizadas son correctas y cuales
descartables.
También será necesario disponer las ecuaciones de los fenómenos físicos de los
modelos, de manera que el ordenador sea capaz de resolverlas, ya que muchas de ellas
son ecuaciones diferenciales en derivadas parciales, muchas con términos no lineales,
que dificultarán su resolución. Otro aspecto a discretizar será la geometría de estudio,
cuya simplificación ayudará a la resolución del fenómeno. Aunque el campo fluido sea
continuo, las ecuaciones no pueden ser resueltas en tantos puntos de fluido como
queramos, así que será, necesario discretizar el dominio en determinados puntos, donde
serán resueltas las ecuaciones, a este paso se le llama discretización espacial.
135
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Por otro lado, cuando el fenómeno está regido por ecuaciones que dependen del
tiempo, tampoco es posible el estudio del fenómeno de forma continua, así que será
necesaria una discretización temporal.
Por último, de las ecuaciones de gobierno se eliminarán aquellos términos que no
sean necesarios o despreciables, con lo que tendremos un nivel más de aproximación, al
que llamaremos aproximación dinámica.
Con todo ello, ya se estará en disposición de resolver las ecuaciones que
modelan determinados fenómenos físicos, empleando distintas técnicas matemáticas.
Soluciones
Soluciones
matemáticas:
aproximadas
modelos
Figura 5.1. Representación de la zona objetivo de las técnicas numéricas.
136
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
errores numéricos podemos señalar: la transformación del sistema de coordenadas, la
discretización, el acoplamiento velocidad-presión, la disipación artificial, los procesos
iterativos y de mallado no convergentes, aproximaciones geométricas y los redondeos
del ordenador.
Según G. W. Dolphin [13], los procedimientos para cuantificar estos errores
pueden agruparse convenientemente en: documentación, verificación y validación. La
documentación recoge una detallada presentación de las ecuaciones matemáticas y los
métodos numéricos utilizados. La minimización de los errores numéricos a través de
estudios de convergencia y orden de exactitud constituyen la base de la verificación,
mientras que la validación es la determinación de la verdadera magnitud del error
cometido a través de la comparación de los resultados con datos experimentales
suficientemente contrastados.
Se requieren estas tres facetas (documentación, verificación y validación) para
identificar errores en el modelo y reducirlos a través de una mejora en la representación
física del problema o en la representación matemática. No conviene olvidar que la
aplicación rigurosa de estos procedimientos permitirá aumentar la confianza en las
soluciones proporcionadas por la CFD, pasando de las aulas universitarias y Centros de
Investigación a las Oficinas Técnicas.
Tendremos que distinguir entre verificación y validación:
Las características según las normas ISO que deben cumplir un software CFD son
las siguientes:
137
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.5 Aspectos teóricos
138
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.5.1.2 Tipos de fluidos
Existen diversas formas de clasificar los fluidos, de las cuales se mencionan las
más importantes.
Gases y líquidos
Como ya se ha comentado anteriormente los líquidos se caracterizan por poseer
fuerzas de cohesión entre sus moléculas, que provocan una fuerte tendencia a que el
fluido mantenga su volumen constante.
Por el contrario, los gases apenas poseen fuerzas de cohesión, por lo que se
pueden expandir tratando de ocupar el máximo volumen.
Compresibles y no compresibles
Todos los fluidos son compresibles en mayor o menor medida, por lo que en
realidad cuando se habla de fluidos no compresibles se refiere a fluidos en los que esta
propiedad es despreciable.
Un fluido es más compresible cuanto mayor es el cambio que sufre en su densidad
al aplicarle una presión externa. Por tanto, los gases son fluidos muy compresibles,
mientras que los líquidos en general se suelen considerar no compresibles.
Newtonianos y no newtonianos
Un fluido newtoniano es aquel en el que su viscosidad no varía en función de la
tensión que se le aplica. Por tanto, su viscosidad es un valor constante.
En los fluidos no newtonianos, la viscosidad no es constante, por lo que se puede
decir que en cierto modo se comportan como sólidos cuando se les aplica una presión, y
como líquidos cuanto menor es la presión
Temperatura
Es un escalar que representa la actividad interna de una sustancia. Este concepto
está ligado al transporte de energía en forma de calor.
139
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Velocidad
Es un vector que representa la dirección, sentido y magnitud de la rapidez de
movimiento del fluido.
Densidad
Suele denotarse con la letra , y es una magnitud escalar referida a la cantidad de
masa contenida en un determinado volumen de una sustancia.
, (5.1)
Viscosidad
La viscosidad de un fluido es aquella propiedad que determina la cantidad de
resistencia opuesta a las fuerzas cortantes. También se puede definir como la facilidad
que tiene un fluido para fluir cuando se le aplica una fuerza externa.
Figura 5.2. Fluido entre dos placas planas (Fuente:
http://webdelprofesor.ula.ve/ingenieria/djean/index_archivos/INST_Flujo/objetos/visco.JPG).
∝ (5.2)
140
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Por semejanza de triángulos en la figura anterior / se puede reemplazar por
el gradiente de velocidades / . Si se introduce ahora una constante de
proporcionalidad, la viscosidad dinámica o absoluta , el esfuerzo cortante entre dos
capas finas de fluido cualesquiera se puede expresar como:
(5.3)
, (5.4)
/
donde es la viscosidad ‘absoluta’ y la densidad.
(5.5)
, (5.6)
donde .
Las fuerzas superficiales son efectuadas sobre el contorno del fluido por el resto
del fluido, por el recipiente o cualquier otra sustancia, mediante contacto directo.
Considerando una fuerza que actúa sobre una superficie infinitesimal , esta
fuerza puede descomponerse en dos componentes:
, (5.7)
donde es la componente normal a la superficie y la componente tangente.
Figura 5.3. Fuerzas de superficie (Fuente: Sámano Tirado, D. A. & Sen, M., 15 de julio de 2009,
Mecánica de fluidos).
141
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
El esfuerzo se define como la fuerza que actúa en el área unitaria. En este caso se
pueden definir dos tipos de esfuerzos:
(5.8)
(5.9)
, (5.10)
Figura 5.4. Perfil de velocidades en un flujo ideal y real (Fuente: Sánchez Vázquez, J., Simulación
numérica directa en paralelo de las ecuaciones de Navier-Stokes en flujos con capa límite desprendida.
Aplicaciones en instalaciones deportivas con gradas cubiertas, Capítulo 2, Ecuaciones de Navier Stokes.
Fénomenos Turbulentos, pp.31-45).
, (5.11)
142
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
donde es la presión absoluta, la densidad, el volumen específico (1/ ), la
constante del gas y la temperatura absoluta.
143
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Figura 5.5. Flujo laminar y turbulento (Fuente: http://fisica.laguia2000.com/wp-
content/uploads/2011/07/Flujo-Laminar-2.gif).
, (5.12)
144
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.5.2.2 Ecuaciones de Navier-Stokes
0 (5.14)
(5.15)
Componente y:
(5.16)
Componente z:
(5.17)
2 2
(5.18)
145
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
de volumen, es la conductividad térmica, es la deposición de calor por unidad de
masa y son las fuerzas viscosas.
Estas expresiones representan cinco ecuaciones de transporte con siete incógnitas
, , , , , , , por lo que son completadas añadiendo dos ecuaciones algebraicas;
una que relaciona la densidad con la temperatura y la presión:
, (5.19)
y otra que relaciona la entalpía con la temperatura y la presión:
, (5.20)
(5.22)
Componente y:
(5.23)
Componente z:
(5.24)
(5.25)
La capa límite es la zona del campo fluido próxima a un contorno sólido en la que
se manifiestan especialmente los efectos viscosos. Debido a la viscosidad y a la
condición de no deslizamiento, cerca de cualquier contorno sólido aparece un gradiente
de velocidades en la dirección normal a dicho contorno. Fuera de esta capa, los
gradientes de velocidad disminuyen rápidamente y el esfuerzo tangencial viscoso es
pequeño. Este gradiente de velocidades condiciona el intercambio energético entre las
146
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
distintas partículas de fluido con velocidades diferentes, originando vorticidad y
turbulencia.
Figura 5.6. Evolución de la capa límite (Fuente: http://www.wfis.uni.lodz.pl/edu/Proposal/image085.gif).
Figura 5.7. Flujo laminar y turbulento sobre una placa plana lisa (Fuente:
http://3.bp.blogspot.com/_ZZuvUtGsTew/Sj4Hx1hPe2I/AAAAAAAAAGE/Xve84bp2HiA/s320/figura+
1.jpg).
147
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.5.3.1 Desprendimiento de la capa límite
Figura 5.8. Evolución de la capa límite hacia un punto de separación (Fuente: Sánchez Nieto, M. M.,
2013, Apuntes Mecánica de Fluidos. ETSINO, Cartagena).
148
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.6 Modelos físico-matemáticos empleados en CFD’s
Los métodos que simulan el perfil de la ola alrededor del casco para la
determinación de la resistencia por formación de olas y para el estudio de la superficie
libre pertenecen al primer grupo. El segundo grupo incluye los métodos que estiman la
resistencia de fricción, en algunos casos la resistencia total y las estelas.
Las ecuaciones de conservación de masa y momentum resultan más complejas de
lo que parecen ya que existen términos no lineales, acoplados y difíciles de resolver. La
experiencia ha demostrado que las ecuaciones de Navier-Stokes describen con exactitud
el movimiento del fluido, sin embargo, sólo han sido capaces de resolverse unos pocos
casos en geometrías simples y cuya aplicación a la industria resulta limitada.
0 (5.28)
149
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1 0 (5.29)
Continuidad:
0 (5.31)
Momentum:
(5.32)
Es un tipo de modelos muy similar a los del fluido ideal. Consiste en reducir el
número de variables que intervienen en los cálculos introduciendo la vorticidad en las
ecuaciones de cantidad de movimiento y de la energía. No se consideran aquí ni las
pérdidas por viscosidad en la capa límite ni los efectos de la turbulencia.
150
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
- El flujo es estacionario.
- Los gradientes de esfuerzos viscosos son despreciables en la dirección de las
líneas de corriente.
- Los gradientes de presión en la dirección de las líneas de corriente se aproximan
por su valor en zonas de capa límite cercanas. La aplicación de estas técnicas
queda restringida al campo aeronáutico.
151
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.6.7 Modelos para la resolución de las turbulencias
152
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1.6.7.3 Modelos URANS (Unsteady Reynolds Averaged Navier-Stokes)
Son modelos matemáticos muy similares a los métodos RANS, sin embargo, en
este caso trabajan en régimen transitorio en lugar de permanente como lo hacen los
métodos RANS. Así la solución del problema depende del tiempo. Este tipo de métodos
está sufriendo gran desarrollo en la actualidad.
1
, lim , (5.37)
→
153
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1
(5.39)
Figura 5.9. Promediado de una ecuación (Fuente: Gutiérrez Romero, J. E., 17 de enero de 2011,
Introduction to CFD in Naval Hydrodinamic: A little review).
154
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1
La turbulencia causa la formación de torbellinos de muchos tamaños diferentes.
A bajas velocidades el flujo es laminar, el fluido es suave, puede sin embargo producir
vórtices de gran tamaño, pero conforme la velocidad se incrementa, en un punto se
pasará a tener flujo turbulento. Una vez este flujo turbulento está totalmente
desarrollado, estos grandes vórtices se contraen a lo largo de la dirección perpendicular
a la vorticidad, y se alargan en la dirección de la vorticidad para ser capaces de
conservar el momento angular, como se puede observar en la siguiente figura.
La complejidad de este fenómeno hace que entender la turbulencia, su magnitud,
predicción, simulación y control se haya convertido en uno de los problemas más
complejos e importantes de la ciencia e ingeniería.
1
Los torbellinos son vórtices de fluido y la corriente inversa que generan cuando el flujo pasa un obstáculo.
155
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Este tipo de modelos resultan demasiado simples para su utilización en
simulaciones en general, por ello es necesario tener un especial cuidado de utilizarlo en
geometrías complicadas o procesos complejos, ya que los resultados no serían válidos.
Entre otros modelos algebraicos se pueden encontrar:
1. Modelo de Cebeci-Smith
2. Modelo de Baldwin-Lomax
3. Modelo de Johnson-King
4. Modelos dependientes de la rugosidad
5. Modelos de longitud de mezcla
Dentro de este punto se incluyen tanto los modelos de una ecuación de cierre y de
dos ecuaciones. Ambos modelos retienen las aproximaciones de Boussinesq para la
viscosidad turbulenta pero difieren en un importante aspecto, mientras que los modelos
de una ecuación resultan incompletos ya que relacionan la longitud de escala turbulenta
con alguna dimensión típica del flujo. Por el contrario, los modelos de dos ecuaciones
provienen de una ecuación más para la longitud de escala turbulenta, lo que equivale a
decir que presentan un modelo completo.
156
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Si se estudia que representa cada uno de los términos se tendría que:
- es la variación de la energía cinética turbulenta,
Más recientemente Baldwin y Barth y por otro lado Spalart y Allmaras han
desarrollado métodos más sofisticados, así, por ejemplo, el modelo de Baldwin-Barth
incluye siete coeficientes de cierre, dos funciones empíricas de amortiguamiento y
alguna otra función que describe la longitud de escala turbulenta.
Originalmente resultaban efectivos con bajo número de Reynolds; en los cuales se
requería una buena resolución de la región viscosa afectada por la capa límite.
157
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Estos modelos están implementados usando funciones de contorno en donde la
resolución de la malla no es suficientemente buena, lo que los hace ideales para obtener
simulaciones en donde la malla no es suficientemente fina y en donde el cálculo
turbulento no es demasiado crítico.
Son relativamente modernos y no se sabe con certeza si son aplicables a cualquier
tipo de fluido complejo en ingeniería; así por ejemplo: no se puede predecir la caída de
los flujos turbulentos, isotrópicos y homogéneos. Además estos modelos son criticados
por la imposibilidad de adaptarlos rápidamente a los cambios de escala ya que en
configuraciones geométricas complejas resulta difícil describir , porque es dependiente
de cantidades no locales como del espesor de capa límite, desplazamiento del espesor...
e introduce una incertidumbre algo similar al modelo de ecuación algebraica.
Este tipo de modelos de turbulencia han sido los más estudiados durante las
décadas de los 80 y 90 de siglo XX, gran parte de los estudios realizados en estos años
son con este tipo de modelos, los cuales permiten una mayor precisión en flujos
complicados, pero que conllevan un gasto computacional elevado, aunque actualmente
esto ya no resulta un inconveniente por el desarrollo informático.
Este tipo de modelos provienen no sólo de una ecuación para el cálculo de sino
que también nos ofrecen ecuación otra para el cálculo de la longitud de escala turbulenta
. Como consecuencia tenemos modelos que permiten predecir los flujos turbulentos
con más precisión. Como punto de partida de estos métodos, tenemos la hipótesis de
Boussinesq junto con la ecuación de energía cinética vista anteriormente. Kolmogorov
en los años 40, definió entonces la segunda variable, a la que se llamó rango específico
de disipación . Esta cantidad tiene dimensiones de , así de esta manera los valores
de la viscosidad turbulenta, longitud de escala turbulenta y disipación turbulenta pueden
ser determinadas como:
/
~ ~ ~ (5.45)
Hemos de notar que este tipo de modelos de turbulencia no son una herramienta
universal para la obtención de las características turbulentas del fluido, dando resultados
más o menos aceptables dependiendo de dónde se utilicen.
A continuación nos centraremos en un modelo de dos ecuaciones que resulta
ampliamente utilizado en diversas simulaciones, este modelo es el .
Como modelo de cierre o estrategia numérica para resolver de forma aproximada
las ecuaciones de Navier-Stokes, se desarrollan dos ecuaciones de transporte
adicionales, una para la energía cinética turbulenta ( ) y otra para la tasa de disipación
de energía cinética turbulenta ( ). Estas variables se definen según las expresiones:
158
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1 1
(5.46)
2 2
2 , (5.47)
donde es la parte fluctuante del tensor de velocidad de deformación.
Las ecuaciones de transporte para y se basan en el conocimiento de los
procesos que producen los cambios en esas variables y son:
2 (5.48)
2 , (5.49)
ó ó (5.50)
, (5.52)
159
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
produce pobres resultados en casos como: flujos no confinados, flujos con grandes
gradientes longitudinales y flujos turbulentos completamente desarrollados en
conductos no circulares.
Una variación del método , denominada , fue publicada por Wilkox en
1993. Se trata de un método RANS en el cual se elimina una variable dividiendo entre
ella todas la demás del problema, simplificando la resolución del problema particular.
Este método fue diseñado para trabajar con modelos a bajo o moderado número de
Reynolds.
Pero da muy buenos resultados en las capas más cercanas a la pared en la capa
límite y da mejores resultados que los modelos en cuanto a la predicción de la
separación de la capa límite. Sin embargo este método tiene el inconveniente de
presentar singularidades, cuando la variable por la que se divide es cero y es muy
sensible a la predicción que se haga de la variable . El modelo es menos
sensible a este valor pero predice peor los resultados cuando hay grandes gradientes de
presión.
Por ello este modelo de turbulencia ( ) resulta más acertado en el campo
aeronáutico que en el naval.
Por tanto, entre modelos de dos ecuaciones de cierre se pueden encontrar:
1. Modelos de k-epsilon
2. Modelos de k-omega
Además de todos estos modelos existen otros que podemos encontrar dentro de la
variedad de códigos existentes en el mercado actual de la CFD:
Modelos v2-f
Modelos de tensiones de Reynolds (RSM)
Modelos de simulación de grandes vórtices (Large Eddy Simulation, LES)
1. Modelo de Smagorinsky
2. Submalla dinámica.
3. Modelo RNG-LES
4. Modelo de adaptación a la viscosidad turbulenta cercana a la pared (WALE)
Modelos de simulación directa
160
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
2. METODOLOGÍA DEL ENSAYO DE REMOLQUE DE UNA
CARENA
Los diferentes módulos están integrados dentro del entorno del programa, siendo
usado éste como uno solo.
En nuestro caso se utilizará el módulo RANSOL para tratar problemas de
dinámica de fluidos y el módulo NAVAL para resolver las ecuaciones de superficie
libre.
1. Pre-proceso
- Creación (o importación) de la geometría a analizar.
- Asignación de condiciones: condiciones de límites, materiales y definición
de datos (velocidad del fluido, rugosidad del casco, tipo de capa límite,
campos de presiones y velocidades, aceleración de la gravedad, etc.).
- Generación del mallado.
2. Cálculo
3. Post-proceso
161
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
• Creación o importación de la
geometría a analizar
• Aplicación de las condiciones
Pre- (condiciones de contorno,
Proceso materiales y definición de
datos)
• Discretización espacial
(mallado)
Figura 5.11. Pasos a seguir en el análisis de un problema de elementos finitos.
2.1.1 Pre-proceso
162
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
2.1.1.1.1 Creación del volumen de control
El tamaño del volumen de control dependerá de la eslora del buque: así para la
manga y el calado el tamaño será de aproximadamente de un 80%·L, aguas arriba (a
proa) un 70%·L y aguas abajo (a popa) un 140%·L.
A continuación se procederá a la identificación de las entidades geométricas, tales
como la obra viva (superficie mojada), obra muerta, superficie libre, superficies de
control y el volumen de control (entre otras).
En los casos en los que sea necesario se hará una unificación de superficies,
obteniendo formas simples de las superficies NURBS.
A continuación una vez que ya se tiene definida la geometría de estudio, junto con
su volumen de control, se pasará a la asignación de las condiciones de contorno
necesarias para la resolución del problema. Las condiciones son todas las propiedades
de un problema, (excepto materiales y límites) que pueden ser asignadas a una entidad
(punto, línea, superficie o sólido). Definiendo así las condiciones básicas de los límites
de un problema.
Tdyn dispone de múltiples condiciones de contorno referentes a cada uno de los
módulos disponibles. Sólo se ha hecho referencia a las opciones propias del RANSOL y
del módulo NAVAL. Dentro de estas condiciones no se engloban las condiciones de
contorno incluidas en las propiedades del fluido.
Las condiciones de contorno necesarias para la resolución de los problemas son
las que se muestran a continuación.
163
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
cierre para evaluar la viscosidad turbulenta. Si es posible estas variables deben definirse
mediante experimentos iniciales. Si no hay datos experimentales disponibles (como es
este caso) los valores se estimarán siguiendo las suposiciones que se exponen a
continuación. El hecho de que la elección de estas suposiciones sea acertada se deberá
comprobar haciendo diferentes simulaciones aunque en la mayoría de los casos un
campo uniforme de las variables de turbulencia compatible con las condiciones de
contorno es suficiente.
Normalmente un valor apropiado de se suele calcular usando la llamada nivel de
intensidad de turbulencia , el cual se define como el ratio de fluctuación entre la
componente de la velocidad local ( ) con respecto a la velocidad media del fluido ( ):
(5.53)
1 2
(5.54)
3 3
Para flujo sobre aviones, coches o submarinos este nivel de turbulencia debe ser
bastante menor del 1%, normalmente se estimará un valor del 0,3% ( 0,003).
En fenómenos de estudio de la capa límite, como se puede observar en la
siguiente tabla, se puede llegar a alcanzar niveles del 25% ( 0,25).
Utilizando la tabla anterior los valores del nivel de turbulencia cambian con la
velocidad y serán mayores a medida que ésta aumente. Esta tabla se ha aproximado a
una curva polinómica de tercer grado con un coeficiente de determinación de
0,9997, lo que corresponde a un ajuste muy bueno. Dicha curva polinómica tiene
la siguiente expresión
0,1051 0,0018 0,0107 (5.55)
Una vez se tiene se calcula la energía cinética turbulenta de la siguiente
forma:
3
(5.56)
2
La longitud de escala turbulenta la estimamos mediante el ratio de viscosidad
turbulenta/viscosidad, para este valor se suele tomar los siguientes valores:
1 10 (5.57)
Una vez estimado este valor obtendremos el valor de longitud de escala turbulenta
con la siguiente fórmula:
164
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
, (5.58)
Campo de velocidades
Cuanto mayor sea la velocidad, mayor será el número de Reynolds ( ) y el
número de Froude ( ), lo que significa que se necesitan modelos más complejos para
resolver un problema dado.
(5.59)
, (5.60)
Campo de presiones
Abarca la presión inicial y de referencia en el campo de presión dinámico,
pudiendo ser una constante o una función. En este caso se trabaja con valores de presión
relativa, es decir, si se fija la presión en un punto en 0 Pa, en realidad la se está fijando
al valor de la presión atmosférica 101325 Pa.
El campo de presiones debe definirse al menos en las líneas de salida, es decir, las
últimas líneas aguas abajo de la superficie de control. El campo de presiones puede
tener un valor fijo o ser una función del tiempo. En este caso la presión en la línea de
salida tendrá un valor fijo de 0 Pa.
165
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
de forma explícita con el modelo numérico, lo que hacen es que abordan esta zona
utilizando lo que se conocen como funciones de ley de la pared.
Para construir estas funciones la región cercana a la pared es caracterizada en
términos de variables adimensionalizadas con respecto a condiciones en la pared, las
cuales son definidas en términos de la velocidad de fricción en la pared , la cual a su
vez depende del esfuerzo cortante en la misma:
, (5.61)
(5.64)
(5.65)
Con esto se logra que si el flujo cercano a la pared depende de las condiciones en
la pared se puede obtener de e funciones universales, dentro de unos límites de
que permitan obtener algunas de las características de la capa límite; de hecho en la
práctica se comprueba las relaciones entre estas dos magnitudes.
Las funciones estándar de la ley de la pared son válidas para paredes lisas, pero
pueden ser modificadas de acuerdo a la rugosidad de la pared mediante el ajuste de las
constantes de la ‘Ley de la pared’. Si la rugosidad de la pared está siendo modelada, la
distancia a la pared en ésta es no dimensionalizada con una equivalente altura de
rugosidad. Así algunos códigos incluyen una ley específica de la ‘Ley de la Pared’ para
modelar de acuerdo a este efecto (Rough Wall field).
Estas funciones universales pueden ser usadas para relacionar las variables fluidas
en una primera malla de puntos computacional, desplazando alguna distancia desde la
pared, el esfuerzo cortante de la pared es obtenido directamente sin resolver la
estructura situada en medio.
Las funciones estándar de la ‘Ley de la pared’ son una gran fuente de falsas ideas
en la computación de flujos turbulentos, incluso para usuarios con experiencia. El
propósito es obtener una capa viscosa extremadamente delgada cerca de la superficie.
Tdyn utiliza la ley de Reychardt para modelizar la capa límite:
166
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
donde es la velocidad adimensionalizada e la distancia adimensionalizada
normal a la pared, la cual viene dada por:
(5.67)
Los valores de suelen ser un 10% del espesor de capa límite ( 0,1 ), que
puede determinarse, por ejemplo, a través de la Ley de Blasius:
5
, (5.68)
√
donde es el número de Reynolds del problema y es la longitud característica del
problema. Esta longitud se refiere a la longitud de flujo laminar que hay en la proa del
buque. En el caso en que se trate de flujo laminar o muy cercano a ser laminar esta
casi coincidirá con la eslora del vehículo. Sin embargo, si se trata de flujo turbulento la
tendrá un valor aproximado de entre un 10 y un 30% de la eslora de la embarcación.
Figura 5.12. Esquema del comportamiento del flujo hidrodinámico alrededor de una carena (Fuente:
http://www.magisternavis.com/KM7.jpg).
Estas funciones universales pueden utilizarse para relacionar las variables fluidas
en el primer nodo de la malla desplazado una distancia de la pared (punto C en la
figura inferior de la izquierda) directamente a los esfuerzos cortantes en la pared sin la
necesidad de resolver la estructura que se interpone. De hecho, estas funciones pueden
ser usadas para calcular los esfuerzos cortantes en los puntos más cercanos a la pared
(puntos A, B y C en la figura de la izquierda). Este valor puede ser tomado como una
condición de contorno para la solución del flujo en el resto de nodos (puntos C, D y E
en la figura inferior de la izquierda) empleando el método numérico del modelo de
turbulencia correspondiente.
En conclusión, la implementación de la ley de la pared consiste en resolver los
puntos más cercanos a la pared (A, B y C en la figura inferior de la izquierda) usando la
formulación analítica (universal pero aproximada) expresada arriba (la fórmula de
Reychardt), mientras que el método numérico estándar es usado para el resto de puntos
(C, D y E).
167
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Figura 5.11. Dominio real cerca de la pared (izquierda) y aproximación de la ley de la pared (Fuente:
Compass Ingeniería y Sistemas, Turbulence Handbook. Environment for Multi- Physics simulation,
including Fluid Dynamics, Turbulence, Heat Transfer, Advection of Species, Structural mechanics, Free
surface and user defined PDE solvers).
Sin embargo cuando se utiliza la ley de la pared el flujo calculado debe ser
consecuente con las asunciones hechas para llegar a las ecuaciones de pared, ya que
fuera del rango de aplicación de estas funciones; estas funciones pueden ser
inadecuadas. Además se debe señalar el valor de dentro de los primeros puntos de la
malla para validar las funciones de pared, las funciones de la ley de la pared no liberan
al usuario de la necesidad de una adecuada resolución de una porción turbulenta de la
capa límite y para generar una malla con un distribución de pre-especificada resulta
difícil, ya que de depende la solución.
Así será necesario en el modelado de la capa límite turbulenta chequear el límite
inferior (el cual se sitúa entre = 20 – 30) y superior de . Nótese que si el número
de Reynolds es moderado y la capa límite sólo se extiende en entre 300 y 500, no
hay oportunidad de resolver con exactitud la capa límite; si el primer punto de
integración está emplazado con un valor de más pequeño de 100.
Además habrá que comprobar la resolución de la capa límite (si los efectos de la
capa límite son importantes, es recomendable chequear la resolución de la capa límite
después de la simulación) así como que las funciones de pared han sido usadas
correctamente de acuerdo con la rugosidad de la misma.
Superficie libre
El problema de la superficie libre, es uno de los principales obstáculos en la
computación con CFD en hidrodinámica marina. Este problema puede estudiarse desde
diversos puntos de vista: flujo potencial o flujo viscoso.
El problema de la superficie libre puede ser considerado como un caso particular
de un problema más general, la predicción de la interfase entre dos fluidos inmiscibles:
un líquido que fluye (por lo general agua) y aire. El cálculo de la interfase entre dos
168
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
fluidos inmiscibles resulta dificultoso porque no son conocidos a priori ni la forma ni la
posición de la interfase. Básicamente, existen dos enfoques para el cálculo de
superficies libres: uno con flujo potencial y otro con flujo viscoso y dentro de éste
último grupo encontramos los métodos de seguimiento de la interfaz y los de captura de
la interfaz. El primero calcula el movimiento de las partículas de flujo basado en el
enfoque de la función de Lagrange, dónde el dominio numérico se adapta a la forma y
posición de la superficie libre. Aquí aparecen diferentes técnicas numéricas, como por
ejemplo los métodos SPH (Smooth Particle Hydrodinamics) y los métodos PFEM
(Particle Finite Element Methods). En los métodos de seguimiento de la superficie
libre, ésta se considera como un límite del dominio de cálculo dónde las condiciones de
contorno cinemáticas y dinámicas son aplicadas. El principal problema de este enfoque
son los grandes esfuerzos computacionales que son requeridos debido a la necesidad de
actualizar el dominio de análisis de cada paso de tiempo y la dificultad en la imposición
de continuidad de la masa de una manera precisa.
Los métodos de captura de la interfaz consideran a ambos fluidos como un único
fluido con propiedades variables. La interfaz es considerada como una región de cambio
repentino en las propiedades del fluido. Este enfoque requiere un modelado preciso del
salto en las propiedades de los dos fluidos teniendo en cuenta que la superficie libre se
puede mover, doblar y volver a conectar de manera arbitraria, siendo necesario además
que la imposición de las condiciones de contorno en la interfaz sea simplificada.
En este grupo podemos encontrar técnicas como: métodos de marcadores y celdas
(Marked and cell); métodos VOF (Volumen Of Fluid) dónde los fluidos son
considerados como uno y la posición de cada fase del fluido se describe mediante la
asignación de una fracción del volumen total con un número entre 0 y 1; métodos level
set, aquí se define un función de posición de la superficie libre desde el nivel de
referencia para cada una de las celdas, pudiendo tomar valores positivos o negativos.
Por último podemos encontrar técnicas híbridas que tratan de sobrellevar los
aspectos negativos de cada una de las anteriores y mejorar en los positivos, así
encontramos la técnica ODDLS (Overlapping Domain Decomposition Level Set) dónde
se introduce el concepto de solapamiento de dominios, que permite mejorar la precisión
en la captura de la superficie libre ahorrando coste computacional.
Modelización de la popa
Como se observa en “Un método de elementos finitos para análisis
hidrodinámico de estructuras navales”, de Julio García Espinosa [14], hay un aspecto
geométrico que presenta dificultades en la simulación de los fenómenos que se
desarrollan en presencia de una popa con espejo. En la figura siguiente se pueden
observar los tres tipos de flujos que se pueden encontrar en presencia de una popa de
espejo. La aparición de uno u otro dependerá de la configuración del buque en cada
momento.
El flujo que se ha denominado de popa tradicional es similar al que se produce en
popas de crucero y no presenta dificultades especiales, debido a que el espejo queda
sobre la línea de flotación inicial y no afecta en forma alguna. Sin embargo, las popas
habituales en la mayoría de los buques modernos presentan los dos tipos de flujos
adicionales que se presentan en la figura.
El flujo con aguas muertas aparece para velocidades bajas, desarrollándose una
zona de flujo turbulento localizado a popa del buque. Sin embargo, a medida que se
169
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
aumenta la velocidad, la popa de espejo comienza a secarse hasta llegarse a una
situación como la indicada bajo el nombre de popa de espejo, siendo en esta
configuración en la que suele situarse la condición de diseño.
Figura 5.12. Esquemas de los posibles regímenes de flujo que pueden aparecer en una popa de espejo
(Fuente: García Espinosa, J., 2000, Un método de elementos finitos para análisis hidrodinámico de
estructuras navales, Universitat Politècnica de Catalunya).
Todas las situaciones referidas pueden ser perfectamente simuladas. Sin embargo,
existe una problemática asociada a la resolución de la situación típica de la popa de
espejo.
En este caso la superficie libre que se observa experimentalmente presenta dos
zonas claramente diferenciadas: una que se extiende desde la popa de espejo hacia
aguas abajo y otra está formada por el resto de la superficie libre.
En esta situación, el problema que se presenta en la resolución del problema de
superficie libre, por su carácter de ecuación de convección pura, es la necesidad de
imponer una condición de contorno para que la solución esté definida. En este caso es
pues necesario prescribir la elevación de la superficie libre de manera que coincida con
la salida del flujo fluido, como se observa a continuación.
170
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Figura 5.13. Situación de la zona donde es necesario prescribir en los flujos con popa de espejo
(Fuente: García Espinosa, J., 2000, Un método de elementos finitos para análisis hidrodinámico de
estructuras navales, Universitat Politècnica de Catalunya).
Tal prescripción es suficiente para poder simular los fenómenos en este tipo de
flujos, sin necesidad de añadir ningún artificio especial. Sin embargo, la necesidad de
prescribir el valor de la elevación de la superficie libre, complica la simulación del
proceso de transición de flujo con aguas muertas a flujo en popa de espejo. La manera
de proceder en estos casos es partir de velocidades elevadas, donde es seguro que no
aparezca el flujo con aguas muertas, para ir descendiendo la velocidad gradualmente.
De esta manera se llegará a un punto donde el algoritmo no converja, que corresponderá
a la situación en que el espejo comience a mojarse, por lo que la imposición hecha sobre
la superficie libre no es ya física. A partir de este punto el problema se resuelve de la
manera habitual, sin ningún tipo de imposiciones excepcionales.
Creación de la cuchara
Se denomina comúnmente ‘cuchara’ al efecto que acompaña a la popa del buque
cuya finalidad es simular la forma de la ola que se produce a popa del mismo para
aquellos con popas de espejo o formas que se cortan abruptamente en popa (a
determinada velocidad la ola generada deja el espejo seco). Esta estructura situada en
popa no varía en nada las formas del barco.
La forma normal de la cuchara suele ser la de una superficie tangente a la parte
inferior de espejo de popa por un lado y a la superficie libre en el otro, con una longitud
proporcional a la velocidad del buque dada por las siguientes ecuaciones:
donde es la longitud de onda característica del problema dada por la teoría de olas
2 (5.69)
Los valores dados por estas fórmulas son orientativos y están basados en la
experiencia.
171
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
La modelización de este efecto no será necesaria para la correcta simulación del
problema a estudiar debido a que, en este caso, al calado de proyecto el agua no ‘moja’
al espejo del casco a ensayar.
2.1.1.3 Materiales
Límites
Los límites fluidos del cuerpo son grupos de condiciones límite, de características
geométricas y de otros datos y propiedades, que identifican un cuerpo como límite de un
líquido para el análisis. Estas propiedades pueden ser asignadas a superficies (2D) y
volúmenes (3D).
Existe también una base de datos que puede usarse para la definición de los
límites del análisis. El usuario también puede crear nuevos límites a partir de los
existentes y asignarlos también.
Tdyn ofrece una gran variedad de condiciones de capa límite. Las más relevantes
son las siguientes:
172
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
bajos valores del número de Reynolds. Aplicar esta condición fuera de su rango
de aplicación puede generar en grandes imprecisiones en la resolución de la capa
límite y la evaluación de fuerzas de fricción.
- DeltaWall: La ley de la pared usada se extiende en el contorno a una distancia y
de la pared. La distancia a la pared debe ser introducida in el campo Delta .
Este se calcula con las fórmulas expresadas anteriormente. La implementación
de esta condición, basada en la ley de la pared de Reychardt, se considera válida
para valores de 0 300, donde 0,1 , donde es el espesor de la
capa límite.
- RoughWall: Esta condición de ley de la pared tiene en cuenta la rugosidad de la
pared, es aplicada a una distancia y de la pared. La distancia a la pared debe ser
introducida en el campo Delta . La implementación de esta condición, basada
en la ley de la pared con rugosidad, se considera válida para valores de 30
100.
- YplusWall: Esta condición de ley de la pared se extiende en el contorno a una
distancia adimensional de la pared . Esta distancia adimensional debe ser
introducida en el campo Yplus. Esta es una implementación lineal simplificada
de la condición de contorno, basada en la ley de la pared de Reychardt. Esta
formulación se considera válida para valores de 30 100. Esta condición
es la más ampliamente usada para modelos sobre todo para elevados números de
Reynolds, debido a que ha sido corregida mediante diversas pruebas
experimentales.
- ITTC Wall: Esta condición de ley de la pared se extiende en el contorno a una
distancia adimensional de la pared . Esta distancia adimensional debe ser
introducida en el campo Yplus. Esta condición es similar a YplusWall, pero está
corregida basada en experimentos numéricos, para coincidir con la fuerza de
fricción predicha por la fórmula de la ITTC-57. Esta implementación la hace
muy útil para análisis navales en tanques de remolque. Esta es una
implementación lineal simplificada de la condición de contorno, basada en la ley
de la pared de Reychardt. Esta formulación se considera válida para valores de
30 100.
Los datos del problema es toda la información requerida para llevar a cabo el
análisis que no trata de ninguna entidad geométrica en especial. Esto es diferente de las
definiciones previas de condiciones y las propiedades de materiales, que son asignadas a
entidades diferentes. Algunos ejemplos de los datos del problema generales son: el tipo
de solución algorítmica usada junto al programa de resolución, el valor de los
incrementos de tiempos entre los pasos, las condiciones de convergencia, etc.
Los datos del problema están organizados en tres grupos diferentes:
173
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
3. Unidades: Este grupo de datos hace referencia a todas las unidades generales
requeridas para identificar los datos del problema.
En este apartado se tratará de especificar cómo debe ser el mallado (tipo y tamaño
de los elementos de malla) así como unos patrones para fijar un tamaño de malla inicial
con el que empezar a realizar los cálculos.
Figura 5.14. Método de los elementos de contorno (Fuente: Gutiérrez Romero, J. E., 2011, Introducción
a los CFD, Presentación Prezi).
174
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
puede llegar a ser importante. Así, cada vez, son menos utilizados, siendo sustituidos
por métodos como los que se presentan en los siguientes puntos.
Figura 5.15. Método de las diferencias finitas (Fuente: Gutiérrez Romero, J. E., 2011, Introducción a los
CFD, Presentación Prezi).
175
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
error en las integrales de la misma manera que en el método de Galerkin, con lo que
método pierde elegancia.
Figura 5.16. Método de los elementos finitos (Fuente: Gutiérrez Romero, J. E., 2011, Introducción a los
CFD, Presentación Prezi).
Figura 5.17. Método de los volúmenes finitos (Fuente: Gutiérrez Romero, J. E., 2011, Introducción a los
CFD, Presentación Prezi).
Métodos espectrales
Los métodos espectrales utilizar el mismo enfoque general que el de diferencias
finitas y métodos de elementos finitos sustituyendo las incógnitas de la ecuación que
rige por series truncadas.
La diferencia es que, cuando los dos métodos utilizan aproximaciones locales, la
aproximación método espectral es válida en todo el dominio. La aproximación es, ya
sea por medio de series de Fourier truncada o por una serie de polinomios de
Chebyshev. La discrepancia entre la solución exacta y la aproximación se aborda con un
concepto similar al de residuos ponderados de los elementos finitos.
176
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Tipos de mallado
Dependiendo del tipo de problema se debe optar por un tipo de mallado u otro.
Los tipos de mallado se pueden clasificar en varios grupos dependiendo de su
distribución (estructurados, no estructurados, o multi-bloque), o según el tipo de
elementos (triangulares, rectangulares, hexagonales, tetraédricos). A continuación se
analizarán los diferentes tipos de mallado atendiendo a su distribución.
Mallado estructurado
El mallado estructurado conduce a resoluciones más simples y rápidas de las
ecuaciones de Navier Stokes. Además se requiere de menos memoria para el
almacenamiento de las variables y la suavización y la ortogonalidad pueden ser
controlados fácilmente.
Los métodos implícitos que utilizan la estructura del mallado estructurado son
métodos de relajación lineal, relajación en el plano.
Sin embargo, el manejo de geometrías complicadas no es flexible, la
adaptabilidad sólo es posible añadiendo o moviendo líneas de mallado (2D) o
superficies de mallado (3D), lo cual no lo hace flexible. Además, el movimiento de
contornos es difícil de manejar.
Debido a estas características, el mallado estructurado no se puede utilizar sólo,
hay que utilizarlo en combinación con el mallado no estructurado. La opción por la que
se suele optar es por un mallado multibloque.
Figura 5.20. Ejemplo de mallado estructurado (Fuente: Gutiérrez Romero, J. E., 2011, Introducción a los
CFD, Presentación Prezi).
Mallado no estructurado
El mallado no estructurado ofrece flexibilidad para manejar geometrías
complicadas, adaptabilidad y movimiento de contornos. La generación del mallado
automática es más fácil, incluso para geometrías de tres dimensiones muy complicadas.
Sin embargo requiere mayor almacenamiento de variables y precisa un
direccionamiento indirecto, las resoluciones son más complicadas en general y la
exactitud es más baja, en general, debido a la falta de suavidad del mallado.
177
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Debido a estas características, este mallado se tiene que usar necesariamente en
zonas cercanas a la superficie del vehículo, ya que la geometría allí suele ser compleja y
el programa no permite usar otro tipo de mallado.
Figura 5.21. Ejemplo de mallado no estructurado (Fuente: Gutiérrez Romero, J. E., 2011, Introducción a
los CFD, Presentación Prezi).
Mallado Multibloque
En este tipo de mallado se combinan los dos mallados anteriores, aprovechando
así los beneficios de ambos tipos de mallados, estructurados y no estructurados, hay
varias formas de construir este mallado, dependiendo de que se utilicen elementos
triangulares o cuadrangulares para las zonas de mallado no estructurado, la zona del
mallado no estructurado utilizará siempre elementos triangulares o tetraedros
dependiendo de que estemos en una superficie o un volumen. Para construir este tipo de
malla hay que dividir las superficies y volúmenes previamente en zonas más pequeñas
para así en cada una asignar la condición que se desee, es importante que la transición
entre las zonas estructuras y no estructuradas sea suave, es decir, que no haya una gran
variación del tamaño de los elementos entre zonas en contacto.
Con el mallado multibloque se pueden manejar geometrías más complicadas, lo
que aumenta la flexibilidad en general, permite maneras fáciles de paralelizar en un
multiprocesador y ahorra memoria en máquinas de secuencias.
Por el contrario, la generación de mallados es difícil (la forma de especificar los
contornos entre bloques de forma adecuada exige alto grado de experiencia) y es
todavía menos flexible que los mallados sin estructurar cuando existe movimiento de
contornos.
178
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Figura 5.22. Ejemplo de mallado multibloque (Fuente: Gutiérrez Romero, J. E., 2011, Introducción a los
CFD, Presentación Prezi).
Generación de la malla
Los programas CFD que emplean métodos FEM o FVM suelen estar formados
por tres módulos: un preprocesador, un módulo de análisis o solución y un post-
procesador. El preprocesador es el encargado de llevar a cabo el mallado del problema.
Así, los mallados se pueden clasificar atendiendo a dos criterios: la forma de definir las
fronteras del dominio y la conectividad entre los distintos puntos del mallado.
En cuanto a la forma de definir las fronteras, están los mallados conformes (no
interpolan en la definición de la frontera) y los no conformes (que no definen la frontera
real sino una aproximación numérica de la misma).
Atendiendo al criterio de estructura de los datos del mallado, existen dos métodos
de generación de mallados que darían lugar a los dos tipos de mallados básicos: los
estructurados y los no estructurados. En los mallados estructurados, los puntos se
colocan siguiendo una red de familias de líneas coordenadas que permiten visualizar la
relación entre unos y otros de forma directa. Este hecho simplifica mucho los
algoritmos respecto a lo que sería un mallado no estructurado.
En cualquiera de los casos (estructurado o no estructurado), el mallado debe
ajustarse a las fronteras de la región a estudiar de forma que las condiciones de contorno
queden representadas con la mayor exactitud posible.
El mallado debe distribuirse localmente de la forma más regular posible, con
variaciones suaves de su densidad. La densidad de un mallado se define como el
número de puntos por unidad de superficie o volumen. La asignación del tamaño de los
elementos generados en el paso del mallado es de importancia crítica y se realizará
dependiendo de los resultados que quiera obtener. En la mayoría de los casos se
asignarán valores de tamaño menor a la carena y determinadas partes de interés e irá
aumentado el valor del tamaño de malla al alejarse de la carena. Los elementos
demasiado grandes pueden conducir a los malos resultados, mientras que los elementos
179
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
demasiado pequeños pueden aumentar dramáticamente el tiempo de cómputo sin
mejorar la calidad del resultado.
La mayor densidad del mallado se debe localizar donde se espere que las
variaciones espaciales de la solución sean mayores y el mallado debería ajustarse
dinámicamente a las variaciones de las variables en la solución del flujo.
Se pueden presentar una serie de dificultades en la generación de la malla como
una excesiva distorsión de elementos. En el MEF se usan elementos básicos (uni, bi y
tridimensionales) de formas sencillas (rectas, triángulos, rectángulos, cubos, etc.) que
hay que transformar en otras formas de configuración más arbitraria.
Se pueden también mezclar de forma incorrecta distintos tipos de elementos (la
mayoría de los programas de elementos finitos dan libertad completa al usuario en la
elección de elementos en la malla, pero no todos los elementos son compatibles entre sí)
además de realizar conexiones incorrectas entre elementos (para que haya convergencia
en el MEF debe haber continuidad.
En el interior de los elementos la continuidad está asegurada, pero en la frontera
entre elementos es el usuario el que debe conseguirla mediante una elección correcta de
la malla, además de conocer la respuesta del elemento y las aproximaciones utilizadas
en su formulación, teniendo en cuenta que los elementos lineales requieren mallas más
finas que los cuadráticos, y éstos más que los cúbicos.
180
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Función geométrica 1 ⋯
Función exponencial
Función especial 1 1 1
(5.71)
20
2 (5.72)
donde es la longitud de la ola asociada a una determinada velocidad y es la eslora
del buque.
(5.76)
181
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Pasos iniciales (Initial steps)
Durante los primeros pasos iniciales se llevan algunos controles en los algoritmos
para estabilizar el problema durante el inicio del problema. Se recomienda definir que
los paso iniciales sean de aproximadamente el 10% del en el caso de no haber
superficie libre, aunque en los ensayos se llegarán hasta valores del 30% por tener
superficie libre.
Esta estabilización en los pasos iniciales se puede realizar mediante la velocidad
(incrementado poco a poco la velocidad del flujo antes de empezar con el cálculo), el
tiempo (va incrementando poco a poco el incremento de tiempo hasta llegar al
incremento con el que se realizarán los cálculos) o una combinación de ambos.
Para la elección del modelo de turbulencia más adecuado hay que tener cuidado
con las simulaciones simétricas, muchas veces se hace esto porque se reducen los
tiempos de computación, sobre todo cuando la geometría es simétrica. Pero hay muchas
aplicaciones para las que un modelo simétrico da soluciones asimétricas, por ejemplo la
estela en un cilindro moviéndose a bajos números de Reynolds, conocido como vórtices
de Von Karman. Esto puede ser un factor importante a la hora de predecir la forma de
flujo más realista posible.
Además, para poder estudiar qué modelo se adapta mejor a cada situación, se
necesita hacer un estudio y ver como varían los resultados con respecto al mallado.
Debería simularse cada modelo de turbulencia con varios tamaños de malla y
comprobar la sensibilidad de los resultados a este cambio.
Por otra parte se deberían estudiar varios modelos de turbulencia y comprobar la
sensibilidad de los resultados a este cambio. Aún no hay modelo de turbulencia que
pueda expandirse a todos los casos, el modelo de turbulencia a usar será diferente
dependiendo de las características del flujo.
El modelo de turbulencia debe ser escogido para cada aplicación en concreto y se
escogerá aquel que genere las predicciones más precisas, comparándolas con datos
experimentales fiables de simulaciones parecidas.
Cuanto más complejo (y posiblemente más realista) sea el modelo, más fina
deberá de ser la malla y esto acarreará mayores gastos. Por lo que normalmente en
ingeniería se debe optar por el modelo que resulte más práctico.
2.1.2 Cálculo
182
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
2.1.3 Post-proceso
Figura 5.23. Ejemplo de convergencia de las fuerzas de presión para la velocidad de 20 nudos (Fuente:
Tdyn).
Parámetros
iniciales
3
2 2
1000 0,1 1 10
20 0,1051
0,0018
0,0107
5
2 0,5
0,075
0,01 0,1 2
Figura 5.24. Parámetros iniciales a introducir en Tdyn.
183
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
3. APLICACIÓN EN UN CASO PRÁCTICO
En este apartado se realizará una descripción de todos los pasos a seguir para
llevar a cabo el estudio mediante la utilización de Tdyn aplicados al casco del buque
arrastrero. Dichos ensayos se realizarán a distintos números de Froude correspondientes
a las velocidades a las que faenará en las distintas condiciones de trabajo.
Figura 5.25. Casco del buque dividido a una altura de 4 metros (Fuente: Tdyn).
184
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
3.1.1 Generación del volumen de control
Coordenada
Punto X (m) Y (m) Z (m)
1 -70 0 4
2 -70 40 4
3 -70 40 -40
4 -70 0 -40
5 80 0 4
6 80 40 4
7 80 40 -40
8 80 0 -40
Tabla 5.2. Coordenadas de los puntos que forman el volumen de control.
Lo siguiente que se hará es unir dichos puntos mediante líneas, sin olvidar los
puntos ‘extremos’. Una vez que se han unido todos los puntos, el siguiente paso es la
creación de superficies necesarias para luego crear el volumen de control.
Figura 5.26. Volumen de control (Fuente: Tdyn).
185
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
A continuación la geometría y el volumen de control deben ser divididas en
diferentes capas. Estas capas reducen la complejidad visual y facilita la visualización de
los resultados obtenidos. Dichas capas son:
Figura 5.27. Identificación de las entidades geométricas en capas (Fuente: Tdyn).
186
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Figura 5.28. Campo de velocidades (Fuente: Tdyn).
Figura 5.29. Campo de presiones (Fuente: Tdyn).
187
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
ayuda de una página web dedicada a la mecánica de fluidos computacional (www.cfd-
online.com), para cada una de las velocidades se ha obtenido el valor de dicho
porcentaje, el cual se muestra a continuación.
Velocidad ‘x’
(nudos) (%)
3 16,5
5 16
10 15
15 13
17,5 12,8
20 12,5
Tabla 5.3. Valores de la longitud del flujo laminar ‘x’ para las distintas velocidades.
188
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
3.2.3 Mallado
a)
189
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
b)
c)
d)
e)
Figura 5.30. Comparación entre los diferentes mallados (Fuente: Tdyn).
Para cada modelo con cada una de las mallas se han obtenido unos valores de la
resistencia en función del número de elementos.
190
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
K‐Epsilon Low Reynolds K‐Epsilon High Reynolds K‐Omega
K‐Omega SST Smagorinsky K‐Energy Two Layers
K‐Epsilon Two Layers
110000
105000
100000
Resistencia (kg)
95000
90000
85000
80000
75000
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.1. Comparación de los distintos modelos de turbulencia empleados en el análisis de
sensibilidad.
Se puede apreciar que los resultados varían mucho conforme se añaden elementos,
pero que llegado un punto estos resultados se estabilizan (a partir de unos 150000
elementos). Se observa además que para unos modelos de turbulencia los resultados
convergen más rápidamente que para el modelo de Smagorinsky y el de K-Energy Two
Layers.
Para el mallado número 4 (aproximadamente 100000 elementos) todavía no se ha
alcanzado dicha estabilidad en los resultados y para mallados más densos que el del
mallado 5 se obtendrían unos resultados similares al del mismo pero con un gasto
computacional mucho mayor sin necesidad. Es por esto que se escogerá el mallado 5
para los ensayos.
191
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
origen en la viscosidad lo hace en menor proporción, por lo que la resistencia de origen
viscoso debe ser mucho menor que la resistencia de presión. Además se tendrá en
cuenta la estabilidad de los resultados según el aumento del número de elementos de
malla.
A continuación se van a analizar los resultados obtenidos con 7 modelos de
turbulencia distintos.
120000
100000
Resistencia (kg)
80000
60000
40000
20000
0
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.2. Curvas de resistencia del modelo de turbulencia K-Epsilon High Reynolds.
192
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
también la distribución de resistencia de presión y viscosa es la que cabría esperar y
tiene una buena convergencia.
120000
100000
Resistencia (kg)
80000
60000
40000
20000
0
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.3. Curvas de resistencia del modelo de turbulencia K-Epsilon Low Reynolds.
120000
100000
Resistencia (kg)
80000
60000
40000
20000
0
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.4. Curvas de resistencia del modelo de turbulencia K-Epsilon Two Layers.
193
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
IV. K-Omega
El modelo K-Omega es un modelo de dos ecuaciones para flujos
turbulentos con integración del contorno de fluido a la pared. El modelo implementado
está basado en la descripción hecha por Wilcox [5] con algunas modificaciones para ser
usado con diferentes condiciones en los límites de la pared.
Para el caso ensayado este modelo de turbulencia obtiene muy buenos resultados
al obtener la resistencia más próxima a los valores obtenidos mediante los métodos de
predicción estadística de entre los otros modelos de turbulencia empleados.
En cuanto a la distribución de resistencias la resistencia de presión es mayor que
la resistencia viscosa, lo que se ajusta a lo predecible, presentando además una buena
convergencia.
120000
100000
Resistencia (kg)
80000
60000
40000
20000
0
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.5. Curvas de resistencia del modelo de turbulencia K-Omega.
V. K-Omega SST
El modelo K-Omega SST es un modelo con dos ecuaciones combinadas en una
expresión que se obtiene a partir de las ecuaciones de energía y de disipación ,
mientras que la otra ecuación se obtiene del modelo , es decir, utiliza las
propiedades de los modelos cerca de la pared y gradualmente mezclándolo con el
modelo conforme se aleja de ella, solucionando el problema de la sensibilidad al
valor de sin perder los buenos resultados que ofrece cerca de la pared,
pudiendo ser usado para todo tipo de contorno fluido.
Para el caso del casco ensayado este modelo de turbulencia obtiene unos
resultados similares a los obtenidos con el modelo K-Epsilon High Reynolds. Además
la distribución de resistencia de presión y viscosa es aceptable, presentando también una
buena estabilidad.
194
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Resistencia total Resistencia de presión Resistencia viscosa
120000
100000
Resistencia (kg)
80000
60000
40000
20000
0
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.6. Curvas de resistencia del modelo de turbulencia K-Omega SST.
VI. Smagorinsky
Este modelo ofrece una resolución de la turbulencia mediante el modelo Large
Eddy Simulation. La implementación incluye una amortiguación de la viscosidad
turbulenta en la capa límite.
Así en este modelo se calculan los remolinos de manera independiente mediante
simulaciones dependientes del tiempo que ‘filtran’ esos remolinos. Esta filtración
consiste simplemente en la manipulación de las ecuaciones de Navier-Stokes y borrar
los remolinos más pequeños de una manera prefijada por la malla.
Este modelo de turbulencia no presenta convergencia, necesitando una malla muy
densa, posiblemente si se redujese aún más el tamaño de malla se obtendrían buenos
resultados pero con un coste computacional mucho más elevado.
Además, aunque la resistencia de presión sólo es del orden del doble de la
resistencia viscosa cuando quizá cabría esperar que la resistencia de presión fuese algo
mayor en proporción.
120000
100000
Resistencia (kg)
80000
60000
40000
20000
0
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.7. Curvas de resistencia del modelo de turbulencia Smagorinsky.
195
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
VII. K-Energy Two Layers
El modelo K Energy Two Layers es un modelo de una ecuación de Prandtl ( )
para flujos turbulentos con integración en la pared, donde la longitud de escala
turbulenta ( ) viene dada por el campo Elen Field. La implementación de este modelo
incluye una amortiguación de la viscosidad turbulenta en la zona de la capa límite.
Al igual que con el modelo anterior no se presenta una estabilidad en los
resultados, los cuales son demasiados elevados respecto a los estimados, y la
distribución de resistencias tampoco es la esperada.
120000
100000
Resistencia (kg)
80000
60000
40000
20000
0
0 50000 100000 150000 200000 250000
Nº de Elementos
Gráfica 5.8. Curvas de resistencia del modelo de turbulencia K-Energy Two Layers.
196
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
197
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
En lo referente a la discretización espacial, el tamaño de los elementos de malla en
el casco es de 0,2 metros. Este valor se ha obtenido probando de forma que el número
de Courant se encontrara entre 0,01 0,1, como se puede comprobar a
continuación.
En cuanto a la discretización temporal hay que recordar que los valores del
incremento de tiempo y del número mínimo de pasos obtenidos serán modificados de
forma que el tiempo de simulación ronde los 120 segundos, por lo que los valores de los
parámetros son los siguientes:
Tabla 5.8. Valores de los parámetros referentes a la discretización temporal para las distintas velocidades.
198
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
3.3 Visión, análisis y comparación de resultados
199
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
100000
90000
80000
70000
60000
Resistencia (kg)
50000
40000
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 5.9. Comparación de las resistencias obtenidas por los distintos modelos de turbulencia seleccionados.
200
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Como se puede observar los valores obtenidos por los modelos K-Epsilon son
superiores a los obtenidos por los K-Omega. Entre estos últimos los resultados son
prácticamente similares, sin embargo, para 20 nudos la diferencia entre las resistencias
totales (unos 5000 kg aproximadamente) es superior a las del resto de velocidades.
Como ya se comentó anteriormente los resultados obtenidos se deberían comparar
con datos experimentales pero, al no disponer de ellos, se escogerá aquel modelo que
más se acerque a los valores hallados por medio de los métodos de predicción
estadísticos. Atendiendo a este criterio el modelo escogido sería el K-Omega. Aun así
los valores obtenidos por este modelo son algo elevados ya que también es posible que
haya cierto trimado dinámico y, como resultado, genere una disminución de la
resistencia viscosa o de presión, que haga que los resultados sean superiores a los
esperados.
Sin embargo, posteriormente se han comprobado la distribución de presiones y la
de velocidades, así como los perfiles de ola, y se ha observado que el modelo que ofrece
unas distribuciones más coherentes es el K-Epsilon High Reynolds.
Es por esta última razón que, aunque las resistencias obtenidas con este último
modelo sean bastante superiores a las de los métodos de predicción estadística, se
utilizará este modelo de turbulencia.
Como se comentó en el Capítulo I, la magnitud relativa de cada una de las
componentes de la resistencia al avance de un barco depende de la velocidad a la que se
navega y de las formas de la carena.
Para valores pequeños de la velocidad la resistencia ocasionada por la viscosidad
predomina sobre la resistencia por formación de olas, llegando la resistencia viscosa a
ser aproximadamente un 61,5% de la resistencia total. A medida que la velocidad crece
la resistencia por formación de olas también lo hace, mientras que la que tiene su origen
en la viscosidad lo hace en menor proporción. Cuando la velocidad llega a ser
suficientemente alta, el valor que toma la resistencia debida a la formación de olas es
tan alto que llega a ser un 77,5% respecto de la resistencia total.
201
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Rw Rv
90
80
70
Resistencia (%RT)
60
50
40
30
20
10
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 5.10. Curvas de la resistencia viscosa y por formación de olas respecto a la resistencia total.
202
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de presiones a 3 nudos
Figura 5.31. Distribución de presiones en el casco a 3 nudos (Fuente: Tdyn).
Figura 5.32. Distribución de presiones en la superficie libre a 3 nudos (Fuente: Tdyn).
Figura 5.33. Distribución de presiones en el volumen de control a 3 nudos (Fuente: Tdyn).
203
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de velocidades a 3 nudos
Figura 5.34. Distribución de velocidades en el casco a 3 nudos (Fuente: Tdyn).
Figura 5.35. Distribución de velocidades en la superficie libre a 3 nudos (Fuente: Tdyn).
Figura 5.36. Distribución de velocidades en el volumen de control a 3 nudos (Fuente: Tdyn).
204
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Mapa de olas a 3 nudos
Figura 5.37. Mapa de olas en la superficie libre a 3 nudos (Fuente: Tdyn).
Figura 5.38. Mapa de olas en la superficie libre a 3 nudos (Fuente: Tdyn).
Figura 5.39. Mapa de olas en el volumen de control a 3 nudos (Fuente: Tdyn).
205
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
0,06
0,04
0,02
Altura (m)
0
‐70 ‐60 ‐50 ‐40 ‐30 ‐20 ‐10 0 10 20 30 40 50 60 70 80
‐0,02
‐0,04
‐0,06
Longitud (m)
Gráfica 5.11. Perfil de olas a una velocidad de 3 nudos.
206
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de presiones a 5 nudos
Figura 5.40. Distribución de presiones en el casco a 5 nudos (Fuente: Tdyn).
Figura 5.41. Distribución de presiones en la superficie libre a 5 nudos (Fuente: Tdyn).
Figura 5.42. Distribución de presiones en el volumen de control a 5 nudos (Fuente: Tdyn).
207
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de velocidades a 5 nudos
Figura 5.43. Distribución de velocidades en el casco a 5 nudos (Fuente: Tdyn).
Figura 5.44. Distribución de velocidades en la superficie libre a 5 nudos (Fuente: Tdyn).
Figura 5.45. Distribución de velocidades en el volumen de control a 5 nudos (Fuente: Tdyn).
208
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Mapa de olas a 5 nudos
Figura 5.46. Mapa de olas en la superficie libre a 5 nudos (Fuente: Tdyn).
Figura 5.47. Mapa de olas en la superficie libre a 5 nudos (Fuente: Tdyn).
Figura 5.48. Mapa de olas en el volumen de control a 5 nudos (Fuente: Tdyn).
209
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
0,10
0,08
0,06
0,04
0,02
Altura (m)
0,00
‐70 ‐60 ‐50 ‐40 ‐30 ‐20 ‐10 0 10 20 30 40 50 60 70 80
‐0,02
‐0,04
‐0,06
‐0,08
‐0,10
Longitud (m)
Gráfica 5.12. Perfil de olas a una velocidad de 5 nudos.
210
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de presiones a 10 nudos
Figura 5.49. Distribución de presiones en el casco a 10 nudos (Fuente: Tdyn).
Figura 5.50. Distribución de presiones en la superficie libre a 10 nudos (Fuente: Tdyn).
Figura 5.51. Distribución de presiones en el volumen de control a 10 nudos (Fuente: Tdyn).
211
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de velocidades a 10 nudos
Figura 5.52. Distribución de velocidades en el casco a 10 nudos (Fuente: Tdyn).
Figura 5.53. Distribución de velocidades en la superficie libre a 10 nudos (Fuente: Tdyn).
Figura 5.54. Distribución de velocidades en el volumen de control a 10 nudos (Fuente: Tdyn).
212
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Mapa de olas a 10 nudos
Figura 5.55. Mapa de olas en la superficie libre a 10 nudos (Fuente: Tdyn).
Figura 5.56. Mapa de olas en la superficie libre a 10 nudos (Fuente: Tdyn).
Figura 5.57. Mapa de olas en el volumen de control a 10 nudos (Fuente: Tdyn).
213
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
0,6
0,4
0,2
0
Altura (m)
‐0,2
‐0,4
‐0,6
‐0,8
Longitud (m)
Gráfica 5.13. Perfil de olas a una velocidad de 10 nudos.
214
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de presiones a 15 nudos
Figura 5.58. Distribución de presiones en el casco a 15 nudos (Fuente: Tdyn).
Figura 5.59. Distribución de presiones en la superficie libre a 15 nudos (Fuente: Tdyn).
Figura 5.60. Distribución de presiones en el volumen de control a 15 nudos (Fuente: Tdyn).
215
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de velocidades a 15 nudos
Figura 5.61. Distribución de velocidades en el casco a 15 nudos (Fuente: Tdyn).
Figura 5.62. Distribución de velocidades en la superficie libre a 15 nudos (Fuente: Tdyn).
Figura 5.63. Distribución de velocidades en el volumen de control a 15 nudos (Fuente: Tdyn).
216
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Mapa de olas a 15 nudos
Figura 5.64. Mapa de olas en la superficie libre a 15 nudos (Fuente: Tdyn).
Figura 5.65. Mapa de olas en la superficie libre a 15 nudos (Fuente: Tdyn).
Figura 5.66. Mapa de olas en el volumen de control a 15 nudos (Fuente: Tdyn).
217
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
1,2
1
0,8
0,6
0,4
0,2
Altura (m)
0
‐70 ‐60 ‐50 ‐40 ‐30 ‐20 ‐10 0 10 20 30 40 50 60 70 80
‐0,2
‐0,4
‐0,6
‐0,8
‐1
‐1,2
‐1,4
Longitud (m)
Gráfica 5.14. Perfil de olas a una velocidad de 15 nudos.
218
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de presiones a 20 nudos
Figura 5.67. Distribución de presiones en el casco a 20 nudos (Fuente: Tdyn).
Figura 5.68. Distribución de presiones en la superficie libre a 20 nudos (Fuente: Tdyn).
Figura 5.69. Distribución de presiones en el volumen de control a 20 nudos (Fuente: Tdyn).
219
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Distribución de velocidades a 20 nudos
Figura 5.70. Distribución de velocidades en el casco a 20 nudos (Fuente: Tdyn).
Figura 5.71. Distribución de velocidades en la superficie libre a 20 nudos (Fuente: Tdyn).
Figura 5.72. Distribución de velocidades en el volumen de control a 20 nudos (Fuente: Tdyn).
220
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Mapa de olas a 20 nudos
Figura 5.73. Mapa de olas en la superficie libre a 20 nudos (Fuente: Tdyn).
Figura 5.74. Mapa de olas en la superficie libre a 20 nudos (Fuente: Tdyn).
Figura 5.75. Mapa de olas en el volumen de control a 20 nudos (Fuente: Tdyn).
221
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
2
1,7
1,4
1,1
0,8
0,5
0,2
Altura (m)
‐0,1
‐70 ‐60 ‐50 ‐40 ‐30 ‐20 ‐10 0 10 20 30 40 50 60 70 80
‐0,4
‐0,7
‐1
‐1,3
‐1,6
‐1,9
‐2,2
Longitud (m)
Gráfica 5.15. Perfil de olas a una velocidad de 20 nudos.
222
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
Comparando los resultados obtenidos en lo que se refiere a la velocidad y la
presión, se puede observar que se cumple el Teorema de Bernoulli, es decir, que las
zonas que experimentan altas presiones registran bajas velocidades y viceversa.
Donde se observa bien este fenómeno es en la zona de proa, donde en la
distribución de presiones posee un color rojo (lo cual representa una zona de alta
presión) y en la de velocidades presenta un color azul (correspondiente a una zona de
baja velocidad).
A lo largo de todo el casco se observa este mismo fenómeno, el cual se refleja
mediante los diferentes colores de las leyendas.
Figura 5.76. Distribución de presiones (izquierda) y velocidades (derecha) en la zona de proa (Fuente:
Tdyn).
20 kn 10 kn
2
1,5
1
0,5
Altura (m)
0
‐70 ‐60 ‐50 ‐40 ‐30 ‐20 ‐10
‐0,5 0 10 20 30 40 50 60 70 80
‐1
‐1,5
‐2
‐2,5
Longitud (m)
Gráfica 5.16. Perfil de olas a una velocidad de 20 y 10 nudos.
223
CAPÍTULO V - MÉTODOS DE ESTIMACIÓN DE POTENCIA MEDIANTE CFD
4. CONCLUSIONES
224
CAPÍTULO VI
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
1. COMPARACIÓN DE MÉTODOS
A continuación se van a mostrar los resultados obtenidos tanto por los métodos de
predicción estadística como mediante CFD.
Comparando estos últimos con el método de Calisal, el cual ya se comentó que se
iba a tomar como el método numérico de referencia en el Capítulo IV debido a que no
se dispone de resultados experimentales y a la forma de la curva para velocidades a
partir de 18 nudos (aparte de ser un método actual en el que las formas base son
similares a los buques actuales), los resultados ofrecidos por Tdyn para el modelo K-
Epsilon High Reynolds son elevados en comparación con el resto de métodos
numéricos en general. Hasta los 5 nudos y para los 20 nudos, los valores de la
resistencia al avance son más parecidos. Sin embargo, para el resto de velocidades
intermedias los resultados distan bastante entre ellos, llegando incluso a los 20000kg de
diferencia (concretamente para 18 nudos).
Añadir que el valor de la resistencia para la velocidad de 3 nudos es la obtenida
con el modelo ILES debido a problemas de convergencia, como ya se indicó en el
Capítulo V.
226
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
Velocidad Van Oortmersen Holtrop Calisal (CP) Calisal (PC) Amadeo (CB) Amadeo (SB) K-Epsilon HR
V RT RT RT RT RT RT RT
(nudos) (kg) (kg) (kg) (kg) (kg) (kg) (kg)
3 196,8 252,3 - - 239,3 236,0 394,9*
4 352,7 433,2 183,8 - 427,7 431,2 -
5 606,1 659,5 428,9 - 688,8 700,0 2018,3
6 969,5 931,8 724,4 - 1048,4 1056,5 -
7 1410,2 1258,5 1082,0 - 1543,5 1519,3 -
8 1896,4 1665,5 1519,0 1424,6 2212,9 2113,0 -
9 2430,1 2208,5 2062,3 2693,5 3081,8 2871,4 -
10 3134,8 2961,6 2752,3 4112,3 4155,3 3841,1 10937,8
11 4110,7 4123,9 3652,2 5743,2 5442,7 5089,2 -
12 6208,3 5816,1 4862,3 7682,5 7008,3 6715,1 -
13 8358,9 7572,7 6542,3 10081,1 9054,6 8872,4 -
14 10986,2 9658,1 8946,9 13176,0 12156,2 11812,9 -
15 16166,7 13132,3 12478,7 17337,5 18436,9 15978,0 29439,3
16 24229,6 19243,0 17765,1 23137,4 - - -
17 33714,8 27884,1 25767,3 31446,2 - - -
18 43005,5 38150,7 37932,1 43568,0 - - -
19 51229,7 48545,5 56398,4 61423,6 - - -
20 58264,0 57035,5 84273,2 87796,1 - - 91315,0
Tabla 6.1. Comparación de los valores de las resistencias obtenidos mediante los distintos métodos de predicción de resistencia.
227
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
100000
90000
80000
70000
60000
Resistencia (kg)
50000
40000
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 6.1. Comparación de las curvas de resistencia obtenidas mediante los distintos métodos de predicción de resistencia.
228
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
2. CONCLUSIONES
229
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
90000
80000
70000
60000
Resistencia (kg)
50000
40000
30000
20000
10000
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Velocidad (nudos)
Gráfica 6.2. Comparación de las curvas de resistencia obtenidas mediante los distintos métodos de predicción de resistencia.
230
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
A continuación se muestra el error que se obtiene tomando como referencia los
resultados obtenidos por el método de Calisal (plena carga). Las velocidades de 3 y 5
nudos no vienen representadas debido a que, para estas velocidades, no se obtenían unos
valores coherentes de los resultados en el método de Calisal.
231
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
3. TRABAJOS FUTUROS
Cabe por último comentar que, de cara a trabajos futuros, hay aspectos del
proyecto susceptibles de mejora o nuevos aspectos que puedan ser estudiados.
Nombrando algunos ejemplos:
- Comparación de buques con forma en ‘U’ y en ‘V’, viendo cuál de ellos tiene
mejor comportamiento.
- Aplicación de otras metodologías.
- Aplicación de distintas soluciones de la superficie libre.
- Comparación de distintos tipos de mallado.
- Validación con experimentación.
- Cálculo de trimado dinámico y ensayos.
- Estudio de remolque.
- Estudio y optimización de una carena mediante el acople de diversos bulbos.
- Estudio de casos en los que sea necesario la creación de una ‘cuchara’.
- Cálculo de los coeficientes hidrodinámicos.
- Estudio de diversos apéndices, timón, hélice, etc.
232
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
4. BIBLIOGRAFÍA
233
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
[19] Gobierno de España, 2 de abril de 1993, Protocolo de Torremolinos de 1993
relativo al convenio internacional de Torremolinos para la seguridad de los
buques pesqueros (1977), Torremolinos.
[20] Gutiérrez Romero, J. E., 17 de enero de 2011, Introduction to CFD in Naval
Hydrodinamic: A little review.
[21] Gutiérrez Romero, J. E., 2011, Introducción a los CFD, Presentación Prezi.
[22] Gutiérrez Romero, J. E., 2015, Práctica 4: Ensayo de remolque de una carena,
ETSINO.
[23] H. Shames, I., 1995, Mecánica de fluidos, 3ª edición, Editorial McGraw-Hill
Interamericana S.A.
[24] Holtrop, J. & Mennen, G.G.J., july 1982, An approximate power prediction
method, International Shipbuilding Progress, vol. 29.
[25] Holtrop, J. & Mennen, G.G.J., october 1978, An statistical power prediction
method, International Shipbuilding Progress, vol. 25.
[26] Holtrop, J., A statistical re–analysis of resistance and propulsion data,
International Shipbuilding Progress.
[27] J. Helmore, P., Update on van Oortmerssen’s Resistance Prediction,
University of New South Wales, Australia.
[28] Micó Barba, F., 1978, Introducción a la pesca marítima, publicación del
ETSIN, Madrid.
[29] Minguito Cardeña, E. y Martínez Rubio, J. L., noviembre de 1990, Las formas
de los buques de pesca: Presente y futuro, XXX Sesiones Técnicas de
Ingeniería Naval, Bilbao.
[30] Palomo López, O., Teoría del Buque, ETSINO.
[31] Pérez Rojas, L., Los CFD como herramienta en Hidrodinámica, Aulas del
Mar, UPCT.
[32] Pérez Rojas, L., Sánchez Sánchez, J. M., Zamora Rodríguez, R., Bermejo, R. y
Souto Iglesias, A., 1996, La calidad en el software. Validación de CFD
(Computational Fluid Dynamics), Madrid.
[33] Pérez Rojas, L., Souto Iglesias, A. y Roca Fernández-Vizarra, P., septiembre
de 2001, Los CFD (Computational Fluid Dynamics) en la formación
hidrodinámica del Ingeniero Naval, IV Simposium Panamericano de
Educación en Ingeniería Naval “SYMMEDUC”, Veracruz (Ver. México).
[34] Rodríguez Martínez, R. M., El buque pesquero de litoral La construcción
naval en Almería. Nuevas tecnologías, Almería.
[35] Sámano Tirado, D. A. & Sen, M., 15 de julio de 2009, Mecánica de fluidos.
[36] Sánchez Nieto, M. M., 2013, Apuntes Mecánica de Fluidos. ETSINO,
Cartagena.
[37] Sánchez Nieto, M. M., abril de 2009, Mecánica de fluidos general, ETSINO,
Cartagena.
[38] Sánchez Vázquez, J., Simulación numérica directa en paralelo de las
ecuaciones de Navier-Stokes en flujos con capa límite desprendida.
234
CAPÍTULO VI – COMPARACIÓN, CONCLUSIONES Y TRABAJOS FUTUROS
Aplicaciones en instalaciones deportivas con gradas cubiertas, Capítulo 2,
Ecuaciones de Navier Stokes. Fénomenos Turbulentos, pp.31-45.
[39] Sander M., Calisal & Dan McGreer, october 1993, A Resistance Study on a
Systematic Series of Low L/B Vessels, Marine Technology, vol. 30, n. 4, pp.
286-296.
[40] Souto Iglesias, A, 2001, Nuevas herramientas de diseño de formas de buques
basadas en códigos de flujo potencial, Departamento de Arquitectura y
Construcción Navales, ETSIN.
[41] V. Giles, R., B. Evett, J. & Lui, C., 2003, Mecánica de los fluidos e hidráulica,
serie Schaum, 3ª edición, Editorial McGraw-Hill Interamericana S.A.
[42] Van Oortmerssen, G., 1971, A Power Prediction Method and its Application to
Small Ships, International Shipbuilding Progress, vol 18, n. 217.
[43] Vázquez, Juan L., 2004, La ecuación de Navier-Stokes. Un reto físico-
matemático para el siglo XXI, Monografías de la Real Academia de Ciencias
de Zaragoza, 26: 31–56, Departamento de Matemáticas, Universidad
Autónoma de Madrid.
[44] Velázquez Sánchez, J. A., octubre de 2010, Arrastrero congelador factoría de
fondo y pelágico, ETSIN, Madrid.
[45] WS Atkins Consultants & members of the NSC, 2004, Best Practice
Guidelines for Marine Applications of Computational Fluid Dynamics.
- www.compassis.com
- www.cfd-online.com
235