Tercer Congreso Internacional de Salicceas en Argentina Disertacin

CORTINAS FORESTALES:
NUEVOS ASPECTOS FLUIDODINMICOS

Jorge Lassig y Claudia Palese

Universidad Nacional del Comahue - Facultad de Ingeniera

INTRODUCCIN

El exceso de viento debe ser reducido a valores tolerables para el cultivo a proteger, de
all que se hace necesario la implementacin de las denominadas Cortinas Rompevientos.
Estas pueden ser naturales o artificiales, entre otras, la principal diferencia entre ambas es
la relacin costos/beneficios: los reparos artificiales requieren una inversin que no se
recupera y duran cierta cantidad de aos pues la exposicin al viento y al sol las deterioran; en
cambio las cortinas forestales tienen la gran ventaja de producir un stock forestal importante.

PRIMEROS ESTUDIOS

Hay numerosos trabajos donde se ha modelizado el efecto reparador al viento de las

cortinas forestales. Un resumen importante hasta la dcada de 1960 est en el Bulletin N 29
del Department of Forestry, Edinburgh University, Forestry Commission, titulado
SHELTERBELTS AND MICROCLIMATE, escrito por J. M. CABORN, y publicado en
Londres en 1957.
En l se presentan estudios realizados en tneles de viento y mediciones de campo de
barreras forestales de diversas especies.
En esos estudios se pueden encontrar conclusiones muy importantes, como la diferencia
entre el patrn de flujo que se desarrolla detrs de un cerco impermeable (figura 1) y el que se
produce en otro con porosidad (figura 2).

Figura 1: Flujo desarrollado por una cerca impermeable al ser embestida por el viento
(obtenida de la ref.1).
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Como se puede observar comparativamente entre ambas figuras, el cerco impermeable

genera dos grandes vrtices detrs de ella y acelera considerablemente el viento por sobre la
misma. La que tiene porosidad permite pasar parte del flujo a travs de ella debilitando los
vrtices y reduciendo la aceleracin superior.

Figura 2: Flujo desarrollado por una cerca porosa al viento (obtenida de la ref.1).

Consecuencia de ello se observ, que una barrera forestal menos densa (con mayor
porosidad) protega mayor superficie que otra muy densa (ver figura 3).

Figura 3: La barrera forestal ms densa se comporta como impermeable y provoca grandes

vrtices detrs de la misma, y la atenuacin del viento cubre menos distancia, comparada
con la menos densa.

Tambin llegaron a la conclusin de que es mejor tener una cortina forestal de pocas filas,
frente a otra de gran nmero de filas, como se muestra en la figura 4, donde la barrera forestal
de muchas filas se comporta como un cuerpo slido frente al viento y protege muy poca
distancia detrs de la misma.
Por esa poca, mediciones a campo en Suiza (figura 4) y en Rusia (figura 5), mostraron
atenuaciones del viento para cortinas forestales de distintas densidades, segn se indican en
las respectivas figuras.
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Figura 4: Mediciones realizadas en cortinas forestales en Suiza por Nageli (ref.2).

Figura 5: Mediciones realizadas en cortinas forestales en Rusia por Panfilov (ref.3).

Consecuencia de todos esos trabajos en el pasado, es muy tpico encontrar entre los
extensionistas el grfico orientador de la figura 6, donde se instruye del grado de atenuacin
del viento detrs de una cortina forestal de porosidad media.
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Figura 6: tpico esquema orientador del grado de atenuacin del viento detrs de una
barrera forestal.

ESTUDIOS MODERNOS

Lo que denominaremos Estudios Modernos se relacionan con un mayor conocimiento

del viento dentro de la capa lmite atmosfrica, y de las herramientas nuevas que se aplican
para su modelado, tanto fsico (tneles de viento meteorolgicos) como computacionales
(CFD: Computational Fluid Dynamics).

Capa Lmite Atmosfrica

La distribucin del viento con la altura, se realiza segn las leyes de la capa lmite
atmosfrica:

u* z
u= ln + ( z / L )
k z0

donde, u* es la denominada velocidad de friccin, k el parmetro de von Karman, z0 es la

rugosidad aerodinmica y depende de lo que hay sobre el terreno, y (z/L) depende de las
condiciones de estabilidad atmosfrica.

120

100

80
Altura (m)

60

40

20

0
0 2 4 6 8 10
Velocidad del viento (m/s)

Figura 7: distribucin esquemtica del viento con la altura, en las cercanas del suelo.

Tneles de viento meteorolgicos o de capa lmite

Los estudios de cortinas forestales mencionadas en la parte anterior, se realizaron en

tneles de viento aerodinmicos y no meteorolgicos.
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

La principal caracterstica de un tnel de viento aerodinmico es una distribucin de la

velocidad uniforme en la vertical y con baja o nula turbulencia, mientras que los tneles de
viento meteorolgicos reproducen la forma de la capa lmite atmosfrica (variacin del viento
con la altura y con un adecuado nivel de turbulencia).
Por lo tanto los resultados a obtener entre un tipo y otro de tneles de viento, diferirn un
poco. Cabe aclarar que los tneles de viento meteorolgicos nacen en la dcada de 1970 con
la denominada Ingeniera del Viento y recin se introducen en el campo agrcola-forestal a
principios de 1990 (ref. 4).

Figura 8: seccin de trabajo en un tnel de viento donde se ensayan barreras forestales

a escala.

Figura 9: seccin por donde ingresa el aire en un tnel de viento meteorolgico para
producir el perfil de velocidades tpico de la capa lmite atmosfrica.

La figura 9 muestra esquemticamente el interior de un tnel de viento meteorolgico en

donde los generadores de vrtices, las pantallas dentadas y los generadores de rugosidad
logran producir un perfil de velocidades al ingresar a la seccin de trabajo del tnel, es decir,
una distribucin del viento tpica semejante a la que se desarrolla en la capa lmite atmosfrica
y con la adecuada turbulencia.
Otra nueva herramienta introducida son los medidores de velocidad de mayor
sensibilidad, antes slo se registraban velocidades medias (y/o presiones medias) a travs de
manmetros diferenciales de columna hidrulica. A partir de la dcada de 1960 se
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

introdujeron los anemmetros de hilo caliente que permiten medir fluctuaciones de velocidad
a frecuencias de hasta 100 Hz. Posteriormente (dcada 1980) comenzaron a desarrollarse los
anemmetros lser que han permitido ver la estructura de los vrtices que se desarrollan en
las estelas. En la dcada de 1990 aparecen los anemmetros snicos.
Todas estas herramientas han producido un salto cualitativo y cuantitativo de la
determinacin del patrn de flujo que se desarrolla detrs de una barrera forestal.

Reparo al viento de una fila de rboles

En forma individual, detrs de un rbol se produce una atenuacin del viento que depende
fundamentalmente de su follaje y dimensiones.

Figura 10: velocidades medias detrs de un roble desarrollado, con y sin follaje, siendo
impactado con una intensidad del viento de 5 m/s.

La figura 10 representa mediciones de velocidad promedio detrs de un rbol de roble

adulto con y sin follaje.
Pero realmente detrs de un solo rbol el patrn de flujo es ms complejo, se establece
una estela vorticosa consecuencia de la diferencia de presiones entre adelante y detrs del
mismo. La figura 11 ejemplifica esto, dicha figura ha sido realizada por Bodo Ruck del
Laboratorium fr Gebude und Umweltaerodynamik, Institut fr Hydromechanik, Universitt
Karlsruhe, Kaiserstr, Alemania, donde han desarrollado importantes tcnicas en mediciones
de viento con lser (ref.5).
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Figura 11: patrn de flujo en la estela de un rbol tipo confera de gran densidad (ref.5).

Parte de los vrtices inferiores pueden ser barridos si la zona cercana al suelo est libre,
como ejemplo, la figura 12 lo ilustra y fue tomada de la referencia anteriormente mencionada
(ref.5).

Figura 12: obtenida de la referencia 5, por CFD, muestra el campo de velocidades medias
alrededor de un rbol semejante al roble.

Si se colocan 3 rboles similares distanciados a un dimetro de copa, entonces, el campo

de vientos ser como el indicado por la figura 13; de all a construir una barrera forestal es un
paso simple.
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Figura 13: obtenida de la referencia 5, por CFD, muestra el campo de velocidades medias
generadas por un conjunto de rboles similares y alineados entre s formando una pequea
barrera forestal.

Si la cortina forestal tiene una porosidad del 50%, entonces se puede obtener el grfico de
la figura 14, donde las curvas isotacas corresponden a una altura igual a la de la cerca (H).

Figura 14: curvas isotacas alrededor de una barrera forestal con un 50% de porosidad.

El viento de por s es turbulento, y el modelo ms simple para interpretarlo es expresar a

la velocidad instantnea u, como la suma de una velocidad media ms una velocidad
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

fluctuante u:

u = u + u
As se pueden reproducir en los siguientes grficos la distribucin de velocidades medias
en un plano perpendicular al suelo y axial al eje de la barrera forestal (figura 15a), y otro que
indique el grado de turbulencia en el mismo plano (figura 15b) donde se adimensionalizan las
isotacas con U / UH y U / UH, siendo UH la velocidad del viento sin perturbar.

Figura 15: superior (a) isotacas de velocidades medias adimensionalizadas en un plano

vertical y coaxial al eje de la barrera forestal; inferior (b) isotacas de la velocidad fluctuante
adimensionalizada en el mismo plano.

De esta forma se tiene una informacin ms completa, no slo de cmo se atena la

Tercer Congreso Internacional de Salicceas en Argentina Disertacin

velocidad del viento detrs de la cortina forestal, sino que adems se indica el grado de
turbulencia presente detrs de la barrera.
Existe una interaccin entre modelos CFD - Tnel de Viento - Mediciones de Campo, las
tres herramientas se deben emplear en estos estudios. Si bien los modelos computacionales
CFD son muy vistosos, la puesta a punto de ellos requiere una gran dedicacin y deben ser
calibrados con datos reales. Como las mediciones a campo son costosas y las variables en
juego no dependen del experimentador sino de la naturaleza, el uso del tnel de viento ayuda
a representar los datos obtenidos con bastante precisin respecto de la realidad.
En sntesis, datos de campo ayudan a verificar los ensayos en los tneles de viento, y los
datos de los tneles de viento ayudan a calibrar los modelos CFD.

ASPECTOS DE INTERS EN LA REGIN NORPATAGNICA

Con las modernas herramientas disponibles es posible realizar estudios con mucha mayor
precisin sobre temas de inters regional. A continuacin se describen dos temas de
importancia para la Norpatagonia.

Rameado en rboles frutales

La exportacin de frutos de pepita en el Alto Valle y ctricos en la Mesopotamia,

representan un importante ingreso para las economas regionales, no obstante, ao a ao se
descarta una cantidad importante de produccin por la mala calidad de la piel de los frutos
causada por el efecto denominado rameado. Esto se traduce en grandes prdidas
econmicas.
La utilizacin de barreras de proteccin elica vivas, que se componen de rboles altos
ubicados en filas a los efectos de disminuir la velocidad del viento dentro del cuadro de frutos,
ayuda a reducir el dao por rameado. Sin embargo, por distintos motivos estas barreras vivas,
bajo ciertas condiciones, parecen no cumplir totalmente con el propsito de reducir los daos
en las frutas.
Actualmente se estudian las barreras elicas vivas de distintas dimensiones y porosidad
con el objeto de evaluar la disminucin de la velocidad del viento y la turbulencia asociada,
en especial la energa que presenta la turbulencia en escalas espaciales significativas en
cuanto al proceso que da lugar al dao sobre los frutos.
Es necesario tener una descripcin ms detallada del tipo de vientos en el interior de la
parcela, adems de su valor medio y direccin, que permita determinar la turbulencia y las
escalas de la turbulencia, con la finalidad de identificar las variables significativas del
problema. Estas variables seran las que permitiran reconocer las caractersticas del viento
que resultan ms perjudiciales para los frutos.
El origen del rameado se puede atribuir a 3 efectos fluidodinmicos:
1) Vibracin de las ramas debido a desprendimientos de vrtices de von Karman.
2) Vibracin de ramas y conjunto de rbol debido a la turbulencia del viento (espectro de
turbulencia del viento incidente o del que logra pasar a travs de la barrera protectora).
3) Efecto aeroelstico (flutter) de flexin de las ramas debido al arrastre que el viento
produce sobre ellas y su restitucin debido a la resistencia a flexin propia de dichas ramas.

Efecto (1): Vibracin debido a desprendimientos de vrtices de von Karman

Al incidir el viento sobre un cilindro detrs del mismo se desprenden vrtices peridicos
conocidos como calles de von Karman. Estos inducen sobre el cilindro vibraciones y si el
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

cilindro es flexible comenzar a vibrar.

La frecuencia de ocurrencia de los vrtices se determina mediante el nmero de Strouhal
(St), que para nmeros de Reynolds prcticos es aproximadamente igual a 0,20.

Fc D
St =
V ,

donde Fc es la frecuencia de vibracin, V la velocidad del viento, y D el dimetro del

cilindro (o rama en nuestro caso).
Por otra parte el nmero de Reynolds se define como la relacin entre fuerzas inerciales y
fuerzas viscosas, expresado en la siguiente ecuacin:

V .D
Re = ,

donde V es la velocidad del viento, D el dimetro del cilindro (o rama), y la viscosidad

cinemtica del aire.
Por lo tanto, es posible determinar las frecuencias de vibracin de ramas de distintos
dimetros sometidas a distintas velocidades medias de viento, y determinar aquellas
velocidades en las cuales las frecuencias se acercan a las frecuencias naturales de vibracin de
las ramas, que dependen de la masa, dimensiones y del mdulo de elasticidad de la madera
viva del rbol.

Figura 16: representacin esquemtica de una rama vibrando y deformndose debido a los
vrtices que se desprenden a sotavento de la rama.
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Efecto (2): Vibracin de ramas debido a la turbulencia del viento

La densidad espectral de turbulencia Su permite visualizar la forma en la cual se

distribuye la energa contenida en la turbulencia entre las diferentes componentes que la
integran. A partir del valor del espectro en el origen (a frecuencia cero) es posible estimar la
escala de longitudes de los vrtices que contienen la mayor parte de la energa de la
turbulencia, es decir, la escala integral de la turbulencia de acuerdo a la siguiente ecuacin:

V .S (0)
L=
4. 2

donde

2 = 2. S v ( f ).df
0

y V es la velocidad media del intervalo.

La curva del espectro adimensionalizado en funcin de la frecuencia presenta un pico,
para una frecuencia fp, con la cual se calcula la escala integral de la turbulencia segn la
siguiente ecuacin:

V
L=
2. . f p

Esta escala suele presentar valores significativamente superiores a las dimensiones de un

rbol. Por tal motivo, el pasaje de un vrtice de tal dimensin solicitara en forma similar a
todas las partes del rbol, y generara respuestas de movimiento en cada una de ellas de
acuerdo a sus caractersticas elsticas. En particular, las partes de mayor porte del rbol se
movern con una frecuencia prxima a su frecuencia propia induciendo impacto en los frutos,
especialmente los de mayor tamao.

Figura 17: Espectro de potencia adimensionalizado de la componente

fluctuante de velocidad.
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Por lo tanto, se deben calcular las frecuencias de vibracin naturales de los rboles y
ramas, segn el mdulo de elasticidad, dimetro y masa.
Con esta informacin se pueden comparar las frecuencias de los vrtices detrs de la
barrera elica con las frecuencias naturales de las ramas de los rboles frutales, si estas
coinciden se amplificar el movimiento de las ramas, y lo que habra que hacer es darle a la
barrera viva la porosidad necesaria para que los vrtices generados sean distintos de las
frecuencias naturales de las ramas.
Esta metodologa se est implementando en el Instituto de Mecnica de los Fluidos e
Ingeniera Ambiental, de la Facultad de Ingeniera, Universidad de la Repblica, Uruguay,
para ctricos (ref.6).

Figura 18: el viento turbulento contiene vrtices de diversos tamaos, los de dimensiones
semejantes a las plantas son los que inducen movimientos globales sobre las mismas.

Efecto (3): Aeroelstico (flutter)

Por ltimo, el viento produce sobre las ramas una fuerza de arrastre FA que tiende a
deformarla. Como el viento es turbulento, las variaciones de velocidad alteran la intensidad de
dicha fuerza, produciendo la reaccin de las ramas para restituir su posicin original, por tal
motivo estas se mueven obedeciendo a una interaccin aeroelstica entre la fuerza
aerodinmica y la resistencia a la flexin de las maderas de las ramas.
Aqu interesa conocer las caractersticas medias y extremas de los vientos, como la
distribucin de Weibull, el desvo Standard, la relacin Vmedia vs. Vmx, etc., para
transformarlas en los valores de FA sobre las ramas.

1
FA = V 2 CA A
2 ,

donde V es la velocidad del viento, la densidad del aire, CA el coeficiente de arrastre de

un cilindro (rama o tronco en estos casos), y A la superficie frontal del cilindro expuesto al
viento.
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

Figura 19: el viento produce una carga aerodinmica sobre todo el rbol
deformndolo segn su capacidad de flexin.

La frecuencia a la cual tender a restituirse el rbol ser:

2 EI
Fc =
2. A

donde es el ngulo flexionado, E el mdulo de elasticidad de la madera viva del rbol, I

el momento de inercia de la rama (o parte del rbol, segn se considere en el estudio), A el
rea de la rama (o parte del rbol) transversal al viento, y la densidad de la madera viva del
rbol.

Cortinas forestales y heladas tardas en el Alto Valle

Por ltimo, comentaremos un aspecto fluidodinmico de las heladas tardas que ocurren
en el Alto Valle de Ro Negro y Neuqun, donde las barreras forestales juegan un rol muy
importante.
Entre los aos 1992 y 1997 realizamos experimentos de campo consistentes en medir la
velocidad del flujo de aire y su temperatura a varios niveles de altura dentro del Alto Valle y
sobre la meseta Sur y Norte (ref.7).
Utilizamos estaciones meteorolgicas automticas y sondeos con un globo cautivo que
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

subamos y bajbamos manualmente, registrando temperatura, humedad y velocidad del

viento.
As pudimos medir la configuracin fluidodinmica en las noches con ocurrencia de
heladas tardas.
Este experimento lo realizamos en varias localidades como: San Patricio del Chaar,
Colonia Valentina (Sur y Norte) en Neuqun y Cinco Saltos, Allen y Gral. Roca en Ro
Negro.

Figura 20: distribucin esquemtica del viento y la temperatura con la altura dentro
del valle cultivado.

La conclusin que obtuvimos de las mediciones realizadas es que en noches con

condiciones para que ocurran heladas tardas, existe en capas bajas de la atmsfera un chorro
de aire con temperaturas entre 5 y 7 C ms caliente que la existente a un metro y medio del
suelo. Este chorro de aire (o jet) de mayor temperatura se ubica entre 25 y 35 m de altura, y
tiene velocidades entre 4 y 7 m/s. Estas caractersticas se pueden observar en el esquema de la
figura 20.

Figura 21: ilustracin del modelo fluidodinmico de dos capas que se producen en el
Alto Valle en noches con condiciones para la ocurrencia de heladas tardas.

La configuracin fluidodinmico queda representada con un modelo de dos capas, como

se ilustra en la figura 21, en la que la capa de abajo es fra y sin velocidad, y la capa superior
(chorro) es ms caliente y con velocidad. Cuando se produce esta superposicin de las capas
de aire pueden generarse ondas de Kelvin-Heltmohltz. La expresin de la velocidad de
propagacin de estas ondas es:
Tercer Congreso Internacional de Salicceas en Argentina Disertacin

1 . u1 + 2 . u 2 g 1 - 2 1 . 2.( u1 - u 2 )
2

c= .
1 + 2 k 1 + 2 ( 1 + 2 )2

siendo c la velocidad de propagacin de la onda, g la aceleracin de la gravedad, k el

nmero de onda (que contiene la longitud de onda), u las velocidades de la capa 1 y 2, y las
densidades de cada capa.
Existen condiciones de velocidad, temperaturas (densidad) y longitudes de ondas en las
cuales la raz cuadrada se hace negativa significando la inestabilizacin de la onda.
Si este chorro o jet de aire caliente por algn motivo se inestabiliza la onda toca el suelo y
evita la helada, mezclando el aire fri estratificado en cercanas del suelo con la masa ms
clida proveniente de esta onda inestabilizada. Si no se inestabiliza avanza el enfriamiento en
el suelo y se produce la helada tarda.

Figura22: imagen de una onda de Kelvin-Heltmohltz cuya longitud de onda coincide con
la separacin entre cortinas rompevientos de cuadro con cultivares de frutales.

Lo interesante de este proceso es que el valor de esas longitudes de onda coinciden con la
separacin tpica de cortinas rompevientos instaladas en el Alto valle (125, 250 m, etc.).
Como este chorro de viento caliente est por encima de los lamos, estos inducen una
longitud de onda en valores que ayudan a inestabilizar el chorro haciendo que este llegue al
suelo y evite la helada.
Por lo tanto, la mejor recomendacin es que no se quiten las cortinas forestales en el Alto
Valle, pues de lo contrario, aumentarn los das con ocurrencia de heladas tardas, con las
sabidas consecuencias econmicas.

CONCLUSIONES

Existen hoy en da poderosas herramientas de medicin y de simulacin, tanto fsica

como computacional, para la determinacin del campo fluidodinmico que se genera en los
alrededores de una cortina forestal. La adecuada comprensin de los efectos del viento (medio
y turbulencia) sobre los rboles podra servir para alimentar a modelos agronmicos y
forestales, mejorando la calidad de sus resultados.

Agradecimientos: A los Organizadores de las Jornadas de Salicceas 2011.

Tercer Congreso Internacional de Salicceas en Argentina Disertacin

REFERENCIAS

1.- Caborn, J. M. (1957). Shelterbelts And Microclimate; Bulletin N 29, Department of

Forestry, Edinburgh University, Forestry Commission, London.

2.- Nageli, W. (1946). Weitere untersuchungen tiber die Windverhaltnisse im Bereich von
Windschutzstreifen, Mit!. sclllveiz. Anst.forstl. Versllchsw. 24(2), (659-737).

3.- Panfilov, Y. (1936). Experiments and Research of All Union Research Institute for
Agri-Forestry Melioration and Forest Management, N 6 (Sheltabelts), State Tech. Printing
Office, Moscow.

4.- Cermak, J. E. (1990). Physical-modeling investigations for wind engineering:

Applications beyond wind-load determinations; Journal of Wind Engineering and Industrial
Aerodynamics, Vol. 36, pp. 39-52.

5.- Ruck, B. (2005). ber Die Aerodynamik Der Bume; Fachtagung Lasermethoden in
der Strmungsmesstechnik; 6 8 September 2005, BTU Cottbus.

6.- Cataldo, J. y Duraona, V. (2007). Wind study, its characteristics and citric fruit
quality. Preliminary results, Twelfth International Conference on Wind Engineering, Cairns,
Australia, 1-6 julio 2007.

7.- Flores, A.; Lassig, J., Palese, C., y Cogliati, M. (2007). Proyecto Mecin: Alto Valle de
Ro Negro, Argentina; Revista Meteorolgica,Vol. 31, N 1 y 2, pp. 17-27. ISSN: 0325-187X.