51

Coeficientes de Influencia

52
Centrmonos ahora en la manera de obtener los coeficientes de influencia. En
nuestro caso esto se reduce al clculo de las integrales

\
|

k elemento
k
ik
d
r
z
_
1

,

\
|

k elemento
k
ik
d
r
z
x
_
0
1

y

\
|

k elemento
k
ik
d
r
z
y
_
0
1

sobre
un tringulo genrico definido por las coordenadas x,y,z de cada vrtice. Se
presentarn dos mtodos de resolucin: el de Hess y Smith y el de Newman.
Ambos mtodos conducen a soluciones analticas de las integrales y realmente
son equivalentes, salvo por la mayor generalidad del mtodo de Newman y su
mayor facilidad a la hora de programarlo y de extenderlo a integrales con
trminos lineales, bilineales y de cualquier grado. Comentaremos ambos
mtodos para el desarrollo de la primera integral, sin trminos lineales y
posteriormente ampliaremos los resultados de Newman para las integrales con
trminos lineales.

El mtodo de paneles implementado en el presente proyecto usa
principalmente el mtodo de Newman, aunque se utiliz el mtodo de Hess
para clculos sobre la estela donde el potencial se distribuye de forma
continua sobre los elementos.
Mtodo de Hess y Smith
2


La primera integral puede ponerse de la siguiente forma:

( )
( ) ( ) [ ]

+ +

k elemento k elemento
k
dy dx
z y y x x
z
d
r
z
_
0 0
2 / 3
2 2
0
2
0
_
1

(4.1)


2
A continuacin se traduce [7] para facilitar la comprensin del lector. Dicho artculo se encuentra aqu
de forma parcial, centrndonos nicamente en la obtencin de los coeficientes de influencia.
53
En esta ecuacin (x,y,z) son las coordenadas del punto donde hemos colocado
la singularidad (punto P). La distancia entre este punto y un punto genrico
del plano (x
0
,y
0
,0) viene dada por r.


Ilustracin 13: Elemento triangular, numeracin de los vrtices

Obsrvese que el orden de los vrtices es favorable al sentido de las agujas del
reloj a diferencia de nuestro sistema de organizacin de los elementos, es
necesario reordenarlos antes de proceder a la integracin.

Vamos a realizar la integral segn el procedimiento seguido en [7]. En este
caso lo que se obtiene es el coeficiente de influencia correspondiente a una
distribucin uniforme (y de intensidad unidad) de manantiales sobre el panel.
El coeficiente de influencia correspondiente a dobletes puede obtenerse por
derivacin.

( ) ( )

+ +
=
elemento elemento
dy dx
z y y x x
d
r
0 0
2 2
0
2
0
1 1
(4.2)
(x
1
,y
1
)
(x
2
,y
2
)
(x
3
,y
3
)
P(x,y,0)
54

Para la integracin se adopta un sistema de coordenadas cilndrico, cuyo eje z
es paralelo al sistema de referencia local del panel. El ngulo es creciente en
el sentido de las agujas del reloj y medido desde una referencia conveniente,
por ejemplo desde el eje x negativo y el origen se encuentra en el punto P.


Ilustracin 14: Definicin del sistema de coordenadas cilndrico

La distancia de un punto cualquiera del panel al punto se puede obtener
como:
2 2
z R r + = (4.3)

Donde R es la distancia media en el plano del panel desde un punto cualquiera
del panel hasta el punto P. Tras estos cambios la ecuacin se presentar del
siguiente modo:

+
R
z R
RdRd
0
2 2

(4.4)

(x
3
,y
3
)
(x
2
,y
2
)
(x
1
,y
1
)
P(x,y,0)

R
55
Donde la integracin se hace desde R=0 hasta un punto del permetro y la
integracin en se realiza recorriendo el permetro en el sentido de las agujas
del reloj.

La contribucin a la integral de cada lado del polgono representar el valor de
la integral para un tringulo definido por los puntos extremos del lado y el
punto P. A medida que el permetro es recorrido en sentido de las agujas del
reloj, el ngulo d es positivo si el punto P est a la derecha del lado del
polgono y negativo si est a la izquierda. Cuando los potenciales de los
tringulos correspondientes a los cuatro lados se sumen para obtener la
integral, las contribuciones de las partes de los tringulos fuera del polgono
sumarn cero y por tanto el resultado ser la integral total.



Ilustracin 15: Influencia de cada lado en la integracin

Ahora de la definicin de r tenemos que:
56

2 2
z R
RdR
dr
+
= (4.5)

Y nuestra integral se convierte en:
( )

=
(
(

d z r d dr
r
z
(4.6)

Como z no depende de la posicin a lo largo del permetro podemos escribir:
k z rd

(4.7)

Donde k valdr cero si se encuentra fuera del tringulo, 2 si est fuera o si
se encuentra justo en el lado.

El primer trmino de la ltima expresin se evala ahora calculando la
contribucin de un lado a la integral y sumando para los tres lados del
tringulo. (Obviamente estos resultados pueden ser generalizados a polgonos
que tengan cualquier nmero de lados). Para expresar la contribucin a la
integral del lado entre los puntos (x
1
,y
1
,0) y (x
2
,y
2
,0) introducimos las
siguientes definiciones geomtricas.

La longitud del lado es:
( ) ( )
2
1 2
2
1 2 12
y y x x d + = (4.8)

El coseno y el seno del ngulo que forma el lado con el eje se calculan de:
12
1 2
12
12
1 2
12
sin
cos
d
y y
d
x x

=
(4.9)

57
Trazamos ahora una perpendicular al lado que pasar por el punto (x,y,0)
(proyeccin de P en el plano), con lo que definimos la coordenada s
12
como la
distancia medida desde el punto de interseccin de la perpendicular con la
prolongacin del lado. El sentido positivo de s
12
va desde (x
1
,y
1
,0) a (x
2
,y
2
,0).
Por lo que para un punto (x,y) del lado del tringulo tenemos que la
coordenada s
12
se define como:
( ) ( )
12 1 12 1 12
sin cos y y x x s + = (4.10)

En particular, la coordenada s
12
asociada a los puntos extremos del lado
(x
1
,y
1
,0) y (x
2
,y
2
,0) es respectivamente:
( ) ( )
( ) ( )
12 2 12 2
) 2 (
12
12 1 12 1
) 1 (
12
sin cos
sin cos
y y x x s
y y x x s
+ =
+ =
(4.11)

La distancia (medida en el sentido perpendicular al lado) entre el punto (x,y,0)
y la prolongacin al lado ser:
( ) ( )
12 1 12 1 12
cos sin y y x x R = (4.12)

Que ser positiva si el punto (x,y,0) est a la derecha del lado con respecto a
la direccin formada por (x
1
,y
1
,0) a (x
2
,y
2
,0) y negativa si est a la izquierda. Las
distancias desde el punto P hasta los puntos extremos del lado son,
respectivamente:
( ) ( )
( ) ( )
2 2
2
2
2 2
2 2
1
2
1 1
z y y x x r
z y y x x r
+ + =
+ + =
(4.13)

Finalmente, la integral puede expresarse en trminos de las siguientes
cantidades
|
|

\
|
+
+ +
=
|
|

\
|
+
+
=
12 2 1
12 2 1
) 1 (
12 1
) 2 (
12 2
12
ln ln
d r r
d r r
s r
s r
Q (4.14)
58
( )
|
|

\
|
+

=
(
(

|
|

\
|

|
|

\
|
=

) 1 (
12
) 2 (
12
2 2
12 1 2
) 1 (
12 2
) 2 (
12 1 12 1
1
) 1 (
12
12
1
2
) 2 (
12
12
1
12 12
tan tan tan ) sgn(
s s z R r r
s r s r R z
r
s
R
z
r
s
R
z
R J (4.15)

Es preferible utilizar la segunda forma de la ecuacin superior, ya que la
primera no est definida a lo largo de la extensin del lado. En la primera
forma de (4.15), las arcotangentes estn definidas en la rama principal ([-/2,
/2]), y en la segunda forma en el rango ([-, ]) .
La contribucin del lado entre los puntos (x
1
,y
1
,0) y (x
2
,y
2
,0) a la integral de
lnea de r es :
12 12 12 12
J z Q R I + = (4.16)

Haciendo lo mismo para todos los lados del tringulo tenemos finalmente que
la integral vale:
+ + z I I I
31 23 12
(4.17)

Para el caso del doblete basta con derivar el anterior resultado con respecto a
z y el resultado ser:
( ) [ ]

=
elemento
J J J k z d x n r
31 23 12 0 0
) sgn( ) ( / 1
r r
r
(4.18)

Donde k valdr 0 si nos encontramos fuera del tringulo, si nos
encontramos justo en el lado y 2 si nos encontramos dentro de l.

Coeficientes de influencia segn Newman
3


J. N. Newman obtuvo una forma ms general de obtener estos coeficientes de
influencia (vase [8]). Se bas en la utilizacin del teorema de Gauss-Bonnet

3
A continuacin se traduce [8] para facilitar la comprensin del lector. Dicho artculo se encuentra aqu
de forma parcial, centrndonos nicamente en la obtencin de los coeficientes de influencia.

59
de geometra diferencial, de tal manera que se evita la integracin directa sobre
la superficie del panel.

Adems, con el desarrollo propuesto por Newman queda claro que la
integracin tambin es vlida para cualquier panel curvilneo cuyos contornos
estn compuestos por segmentos rectos.

Como en el caso de Hess partiremos de un panel cuyos vrtices se ordenan en
el sentido de las agujas del reloj (a diferencia de nuestra numeracin, que
habr que alterar antes de la integracin). El punto de evaluacin P ser un
punto cualquiera del espacio (x,y,z).

Recordemos la integral de partida, el potencial generado por una distribucin
de dobletes constante es el siguiente:

( ) ( ) ( ) ( ) [ ]


+ + =

=
S S S
dy dx z y y x x z d r
z
d n r
0 0
2 / 3
2 2
0
2
0 0
. . / 1 / 1
r
r
(4.19)
Por comodidad prescindiremos del signo negativo para llevar a cabo la
integral.

La ltima forma de la ecuacin puede asociarse con el flujo a travs del panel
debido a una fuente en punto P. Se sigue de esto que el valor de la integral
ser igual al ngulo slido del panel, visto desde el punto P, y estableciendo
que el signo de z ser el mismo que el del ngulo slido.

=
S
r
d n r
3

r r
(4.20)

El ngulo slido del panel puede ponerse en funcin de propiedades
independientes para cada vrtice. Para este propsito se definen tres sectores
con respecto a los correspondientes vrtices, como se muestra en la figura, de
60
forma que la diferencia entre los dominios de el primer y segundo sector, ms
el tercero, es el dominio del tringulo. Si el valor de la integral de superficie
sobre cada sector se define como I
n
, se sigue que:






3 2 1
I I I Total + =




De acuerdo a la definicin tenemos que cada una de las integrales es igual al
ngulo slido del correspondiente sector referido al punto P (donde estamos
midiendo el potencial). Del teorema de Gauss-Bonnet de geometra
diferencial se sigue que el ngulo slido del panel desde el punto P es igual a la
suma:
+ + = + =
3 2 1 3 2 1
I I I Total (4.21)

Donde es el ngulo interior de cada vrtice en el plano tangente a la esfera.
Para un panel rectangular, el ltimo trmino de la ecuacin ser 2, y para un
polgono con N lados ser (N-2) . Obtengamos el trmino
n

correspondiente a un vrtice cualquiera. Para simplificar las ecuaciones que
siguen se asumir que el vrtice est en el origen del sistema de referencia
generalizando el resultado a un vrtice cualquiera en (x
n
,y
n
,0) .

1
2
3
61
El objetivo principal es determinar el ngulo entre dos segmentos
consecutivos (estando definido cada segmento por dos vrtices consecutivos)
del polgono (ntese que los segmentos estn incluidos en el plano z=0 al
estarlo los vrtices). Se define el ngulo
n
como el ngulo polar del segmento
que une el vrtice n con el vrtice n+1.

Para determinar el ngulo
n
visto desde P lo que hacemos es girar el sistema
de referencia hasta hacerlo coincidir con un sistema de referencia (u,v,w). Para
ello se gira un ngulo (el ngulo viene definido por la proyeccin de P en
el plano y el eje x local) alrededor del eje z, y despus un ngulo alrededor
de eje y (ngulo entre P y el eje z), con lo que la transformacin sera:

|
|
|

\
|

=



cos sin sin cos sin
0 cos sin
sin sin cos cos cos
T (4.22)

cos
sin
R y
e R iz x
i
=
= +
(4.23)

De tal forma que el eje w coincide con el radio vector que va desde el origen
hasta el punto P.

Si aplicamos la transformacin a la ecuacin del segmento ( )
n
x y tan =
tenemos la relacin:

( )
( )

=
n
n
v
u
sin
cos cos
(4.24)

62
Lo que representa la proyeccin de un segmento con ngulo polar
n
en el
plano (u,v). Ahora podemos expresar el ngulo
n
(visto desde P) de la
siguiente forma:

( ) ( )
(

=
+

cos
tan
tan
cos
tan
tan
1 1 1 n n
n
(4.25)

Podemos obtener una forma alternativa definiendo un par de rotaciones
diferente, rotando primero un ngulo alrededor del eje y, y despus un
ngulo alrededor del eje z hasta llegar al sistema u,v,w donde el eje v
coincide con el radio vector que va desde el origen hasta el punto P. Ahora la
transformacin nos proporciona la siguiente ecuacin para el segmento:


sin cos
sin sin cos cos cos

=
w
u
(4.26)

Con lo que:
( ) ( )
(

+
=
(

=
+
+
Rz
xy z x
Rz
xy z x
n n
n n
n
2 2
1 1
2 2
1
1 1 1
tan
tan
tan
tan
sin
cos cos sin tan
tan
sin
cos cos sin tan
tan

(4.27)

Donde se ha tenido en cuenta que la transformacin puede expresarse como:

cos
sin
R y
e R iz x
i
=
= +
(4.28)

Debido a que la arcotangente no es nica, se debe determinar la rama correcta
a usar. El ngulo slido es una cantidad comprendida entre -2 y 2, con el
mismo signo que z.

63
Para un polgono con tres vrtices cncavos tenemos que <
n
0 . Para z>
0 los denominadores de la ecuacin son positivos, y con la rama principal de
la arcotangente ([-/2, /2]),
n
tendr el mismo signo que
1
tan tan
+

n n

independientemente de las coordenadas. Para un triangulo normal, dos de
estas diferencias son positivas y la otra negativa. Entonces para una de las tres
integrales hay que aadir a la ecuacin para obtener correctamente el ngulo
slido.

Al cancelarse este con el ltimo valor del teorema de Gauss-Bonnet
podemos reemplazar la ecuacin directamente por la suma de los cuatro
trminos de
n
.

Desde el punto de vista computacional puede ser interesante simplificar la
expresin dada por la ecuacin. Para ello notemos que:
1
1
1 2
1 2
1
1
1
tan
tan
+
+
+ +
+ +
+
+
+
=

=
=

=
n
n
n n
n n
n
n
n
n n
n n
n
x
y
x x
y y
x
y
x x
y y

(4.29)

Sustituyendo y agrupando tenemos que:
( )
( ) [ ] ( )( )

+
=
3
1
2 2
1
0
tan / 1
n n n
n n n n n
S
x z R
y y x x x z x x y
d n r

r
r

( ) [ ] ( )( )
(
(

(
+

+
+ + +
n n
n n n n n
x z R
y y x x x z x x y

1
1 1
2 2
1 1
tan (4.30)

Los pares de arcotangentes en la ecuacin pueden agruparse utilizando la
siguiente identidad trigonomtrica:
|
|

\
|
=
|
|

\
|

|
|

\
|

3
3 1
2
2 1
1
1 1
tan tan tan
c
s
c
s
c
s
(4.31)

64
Donde s
n
y c
n
son los correspondientes factores en el numerador y
denominador de cada trmino de la ecuacin para n=1,2 y s
3
y c
3
vienen dados
por:
2 1 2 1 3
1 2 2 1 3
s s c c c
c s c s s
+ =
=
(4.32)

Si utilizamos esta ltima expresin, la ltima arcotangente de la ecuacin debe
ser evaluada en la rama ([-, ]) sin considerar de forma separada los
argumentos de la arcotangente del miembro de la izquierda.

El primer miembro de (4.31) est mal condicionado numricamente cuando el
punto P est situado directamente encima o debajo de un vrtice, por lo que
es preferible la utilizacin del segundo miembro de (4.31).

En el anlisis anterior se ha supuesto un panel plano situado en el plano z=0.
Sin embargo, los resultados de las ecuaciones son vlidos para el caso ms
general de un panel curvo cuyo contorno son segmentos rectos.

Integrales de coeficientes lineales

Para extender el anlisis precedente e incluir una distribucin lineal de
dobletes en cada panel, consideraremos dos integrales definidas como:
|
|

\
|

S
r
dy dx
y
x
z
3
0 0
0
0
(4.33)

Lo que es equivalente a las integrales

\
|

k elemento
k
ik
d
r
z
x
_
0
1

y

\
|

k elemento
k
ik
d
r
z
y
_
0
1

colocadas en forma matricial.

65
La primera integral puede expresarse de la siguiente manera mediante el
teorema de Green (de nuevo obviamos el signo negativo):

( )
0 0 0 0 0
0
0
3
0 0
0
3
0 0
0
1 1
zxI dy
r
z xI dy dx
r x
z zxI
r
dy dx
x x z
r
dy dx
x z
C S S S
+ |

\
|
= + |

\
|

= + =

(4.34)

Donde I
0
representa el valor de la integral sin trminos lineales que se calculo
en anteriores apartados. Operando de la misma forma para la integral
dependiente de y
0
:
0 0
3
0 0
0
1
zyI dx
r
z
r
dy dx
y z
C S
+ |

\
|
=

(4.35)

Al estar definidos los vrtices y lados en el sentido de las agujas del reloj al
contrario que la definicin de la integral de lnea, debemos variar el signo para
llegar a la siguiente expresin, que presenta el valor de las integrales como la
suma del aporte de cada lado:


|

\
|
|
|

\
|

|
|

\
|
=
|
|

\
|
=
n
Lado S
n n
n
dl
r
z I
y
x
r
dy dx
y
x
z
1
cos
sin
3
1
0
3
0 0
0
0

(4.36)

La ltima integral se resuelve fcilmente siendo su valor:
|
|

\
|
+
+ +
= = |

\
|
+ +
+ +

n n n n
n n n n
n
Lado
d r r
d r r
Q dl
r
n
, 1 1
, 1 1
ln
1
(4.37)

Donde r
n
es la distancia entre el punto P y el vrtice n y d
n+1,n
la longitud del
segmento entre n y n+1.

El resultado final es:
n
n n
n
S S
Q z I
y
x
r
dy dx
y
x
z d n
r y
x

=
|
|

\
|
|
|

\
|
=
|
|

\
|
= |

\
|

|
|

\
|
3
1
0
3
0 0
0
0
0
0
0
cos
sin
1

m
r
r
(4.38)
66

De manera que ya podemos resolver todas las integrales necesarias para
afrontar el problema.