Métodos Matemáticos I

Segundo departamental.
Escuela Superior de Física y Matemáticas
Enrique de Jesús Guzmán Calzadilla [email protected]

Tarea I. Escribir el Laplaciano en coordenadas rectangulares, cilíndricas y esféri-

cas.
∂2u ∂2u ∂2u
Coordenadas rectangulares: ∇2(x,y,z) u = ∂x2 + ∂y 2 + ∂z 2
∂2u 1 ∂u 1 ∂2u ∂2u
Coordenadas cilíndricas: ∇2(r,θ,z) u = ∂r 2 + r ∂r + r 2 ∂θ 2 + ∂z 2
∂2u 2 ∂u 1 ∂2u 1 ∂2u cotθ ∂u
Coordenadas esféricas: ∇2(r,ϕ,θ) u = ∂r 2 + r ∂r + r 2 sen2 θ ∂ϕ2 + r 2 ∂θ 2 + r 2 ∂θ

Tarea II. Escribir la ecuación de onda, difusión y Laplace en coordenadas esfer-

icas.
La ecuación de onda es:

∂2u
a2 ∇2(r,ϕ,θ) u = (1)
∂t2
escribiendo el Laplaciano en coordenadas cilíndricas, tenemos que la ecuación de onda en coordenadas cilindricas es:

2
1 ∂2u ∂2u ∂2u
a2 ( ∂∂ru2 + 1 ∂u
r ∂r + r 2 ∂θ 2 + ∂z 2 ) = ∂t2

La ecuación de difusión en coordenadas cilíndricas es:

∂u
k∇2(r,ϕ,θ) u = (2)
∂t
escribiendo el Laplaciano en coordenadas cilíndricas, tenemos que la ecuación de difusión en coordenadas cilindricas
es:
2
1 ∂2u ∂2u
k( ∂∂ru2 + 1 ∂u
r ∂r + r 2 ∂θ 2 + ∂z 2 ) = ∂u
∂t

Por último, la ecuación de Laplace en coordenadas cilíndricas es:

∇2(r,ϕ,θ) u = 0 (3)

escribiendo el Laplaciano en coordenadas cilíndricas, tenemos que la ecuación de Laplace en coordenadas cilindricas es:

∂2u 1 ∂u 1 ∂2u ∂2u

∂r 2 + r ∂r + r 2 ∂θ 2 + ∂z 2 =0
Xovee Xu 2

Tarea III. Escribir la ecuación diferencial de segundo orden de Cauchy-Euler, y

escribir solución.
La ecuación diferencial de segundo orden de Cauchy - Euler es:

d2 y dy
ax2 + bx + cy = 0 (4)
dx2 dx

Proponemos la solución de la ecuación (4) como:

y(x) = xk (5)

Luego, calculamos:

dy
dx = kxk−1
2
d y
dx2 = k(k − 1)xk−2

Sustituyendo en (4) tenemos:

ax2 (k(k − 1)xk−2 ) + bx(kxk−1 ) + cxk = 0

De donde se sigue:

xk (ak 2 + (b − a)k + c) = 0

O bien:

ak 2 + (b − a)k + c = 0 (6)

Así, (5) es solución de (4) si k es solución de la ecuación auxiliar. Tenmos tres casos a considerar:

• Raíces reales diferentes

Sean k1 y k2 las raíces reales de (6), tales que k1 ̸= k2 . Entonces

y1 (x) = c1 xk1 y2 (x) = c2 xk2

Por el principio de superposición, la solución de (4) es:

y(x) = c1 xk1 + c2 xk2 (7)

• Raíces reales iguales

Si las raíces de (6) son repetidas, es decir k1 = k2 , entonces se obtiene sólo una solución particular

y1 (x) = xk1 = xk2 = xk

Cuando las raíces de la ecuación auxiliar (6) son iguales, el discriminante necesariamente es cero, es así que de (6) se
deduce que las raíces deben ser

k = − (b−a)
2a
Xovee Xu 3

• Raíces complejas

Si las raíces de (6) son el par conjugado

k1 = α + βi k2 = α − βi

Donde α y β > 0 son reales, entonces una solución es

y(x) = c1 xα+βi + c2 xα−βi (8)

Tarea IV. Escribir la ecuación paramétrica de Bessel y su solución.

La ecuación paramétrica de Bessel de orden n es:

d2 y dy
x2 +x + (α2 x2 − n2 )y = 0 (9)
dx2 dx
donde n es un entero no negativo y α es un parámetro positivo.

La solución de la ecuación de Bessel es

y = C1 · Jn (αx) + C2 · Yn (αx)

Series de Fourier-Bessel. Para un número fijo n funciones de Bessel {Jn (αi x)} con i = 1, 2, · · · es ortogonal respecto
a la función peso p(x) = x en un intervalo [0,b] siempre que los αi estén definidos por la condición de frontera

A2 · Jn (αb) + B2 · α · Jn′ (αb) = 0

Los eigenvalores del correspondiente problema de Sturm-Liouville λi = αi2 . La serie ortogonal o serie generalizada de
Fourier del desarrollo de una función f definida en (0, b) en términos de este conjunto ortogonal es:
P∞
f (x) = i1 ci Jn (αi x)dx (2)

donde ci están definidos como:

Rb
xJn (αi x)f (x)dx
ci = 0
||Jn (αi x)||2 (3)

La norma cuadrada de la función Jn (αi x) está definida como:

Rb
||Jn (αi x)||2 = 0
xJn2 (αi x)dx

La serie descrita por la ecuación (2) con los coeficientes descritos por la ecuación (3), se llama Serie de Fourier-Bessel o
bien, simplemente serie de Bessel.

Series de Fourier-Bessel. La serie de Fourier-Bessel en un intervalo (0, b) está dado por:

P∞
f (x) = i=1 ci Jn (αi x)dx

1
Rb
ci = 2
b2 ·Jn+1 (αi b) 0
xJn (αi x)f (x) · dx

donde los αi son las raíces positivas de la ecuación Jn (αb) = 0.

Xovee Xu 4

P∞
f (x) = i=1 ci Jn (αi x)dx

2α2i Rb
ci = (α2i b2 −n2 +h2 )·Jn
2 (α b)
i 0
xJn (αi x)f (x) · dx

donde los αi son las raíces positivas de la ecuación hJn (αb) + αbJn′ (αb) = 0.

P∞
f (x) = c1 + i=2 ci J0 (αi x)dx

2
Rb 2
Rb
ci = b2 )·J02 (αi b) 0
xJ0 (αi x)f (x) · dx c1 = b2 0
xf (x) · dx

donde los αi son las raíces positivas de la ecuación J0′ (αb) = 0.

Relación de ortogonalidad.
Rb
0
xJn (αi x)Jn (αj x)dx = 0 i ̸= j
Rb
0
xJn (αi x)Jn (αj x)dx = 21 Jn+1
2
(αi ) i=j

Derivada de la función de Bessel:

d
dx Jn (αx) = Jn′ (αx) dx
d
αx = αJn′ (αx)

Tarea V. Resolver ecuación de onda por separación de variables.

La ecuación de onda es:

∂2u
a2 ∇2(r,θ,z) u = (10)
∂t2
escribiendola en coordenadas cilindricas es:
2
1 ∂2u ∂2u ∂2u
a2 ( ∂∂ru2 + 1 ∂u
r ∂r + r 2 ∂θ 2 + ∂z 2 ) = ∂t2

Sea u(r, θ, z, t) = R(r) · Θ(θ) · Z(z) · T (t) la solución de la ecuación (10), calculando:

∂2u
• ∂r 2 = R′′ (r) · Θ(θ) · Z(z) · T (t)

• ∂u
∂r = R′ (r) · Θ(θ) · Z(z) · T (t)
∂2u
• ∂θ 2 = R(r) · Θ′′ (θ) · Z(z) · T (t)
∂2u
• ∂z 2 = R(r) · Θ(θ) · Z ′′ (z) · T (t)
∂2u
• ∂t2 = R(r) · Θ(θ) · Z(z) · T ′′ (t)

Sustituyendo en (10), tenemos:

a2 (R′′ · Θ · Z · T + 1r R′ · Θ · Z · T + 1
r2 R · Θ′′ · Z · T + R · Θ · Z ′′ · T ) = R · Θ · Z · T ′′

Dividiendo la ecuación anterior entre a2 R · Θ · Z · T , tenemos:

Xovee Xu 5

R′′ 1 R′ 1 Θ′′ Z ′′ T ′′
R + r R + r2 Θ + Z = a2 T =A

Separando T (t), tenemos:

T ′′
= A −→ T ′′ − Aa2 T = 0 (11)
a2 T
Nos queda:

R′′ 1 R′ 1 Θ′′ Z ′′
R + r R + r2 Θ + Z =A

Separando Z(z), tenemos:

Z ′′ R′′ 1 R′ 1 Θ′′
=A− − − 2 = B −→ Z ′′ − BZ = 0 (12)
Z R r R r Θ
Nos queda:

1 Θ′′ R′′ 1 R′ Θ′′ ′′ ′

r2 Θ =A−B− R − r R −→ Θ = (A − B)r2 − r2 RR − r RR = C

De donde se sigue:

Θ′′
= C −→ Θ − CΘ = 0 (13)
Θ
Separando R(r), tenemos:

r2 R′′ − rR′
(A − B)r2 − = C −→ r2 R′′ + rR′ + [(B − A)r2 + C]R = 0 (14)
R
Seleccionando

• A = −λ2
• B = µ2
• C = −ν 2

De esta forma las ecuaciones (11), (12), (13) y (14) nos quedan como:

T ′′ + λ2 a2 T = 0 T (t) = C1 cos(λat) + C2 sen(λat) (15)

Z ′′ − µ2 Z = 0 Z(z) = C3 cosh(µz) + C4 senh(µz) (16)

Θ′′ + ν 2 Θ = 0 Θ(θ) = C5 cos(νθ) + C6 sen(νθ) (17)

p p
r2 R′′ + rR′ + [(µ2 + λ2 )r2 − ν 2 ]R = 0 R(r) = C7 Jν ( µ2 + λ2 r) + C8 Yν ( µ2 + λ2 r) (18)

La solución entonces es:

p p
u(r, θ, z, t) = [C7 Jν ( µ2 + λ2 r) + C8 Yν ( µ2 + λ2 r)] · [C5 cos(νθ) + C6 sen(νθ)] · [C3 cosh(µz) +
C4 senh(µz)] · [C1 cos(λat) + C2 sen(λat)]
Xovee Xu 6

Tarea VI. Resolver ecuación de difusión por separación de variables.

La ecuación de disfusión es:

∂u
k∇2(r,θ,z) u = (19)
∂t
escribiendola en coordenadas cilindricas es:
2
1 ∂2u ∂2u
k( ∂∂ru2 + 1 ∂u
r ∂r + r 2 ∂θ 2 + ∂z 2 ) = ∂u
∂t

Sea u(r, θ, z, t) = R(r) · Θ(θ) · Z(z) · T (t) la solución de la ecuación (19), calculando:
∂2u
• ∂r 2 = R′′ (r) · Θ(θ) · Z(z) · T (t)
• ∂u
∂r = R′ (r) · Θ(θ) · Z(z) · T (t)
∂2u
• ∂θ 2 = R(r) · Θ′′ (θ) · Z(z) · T (t)
∂2u
• ∂z 2 = R(r) · Θ(θ) · Z ′′ (z) · T (t)
• ∂u
∂t = R(r) · Θ(θ) · Z(z) · T ′ (t)

Sustituyendo en (19), tenemos:

k(R′′ · Θ · Z · T + 1r R′ · Θ · Z · T + 1
r2 R · Θ′′ · Z · T + R · Θ · Z ′′ · T ) = R · Θ · Z · T ′

Dividiendo la ecuación anterior entre k · R · Θ · Z · T , tenemos:

R′′ 1 R′ 1 Θ′′ Z ′′ T′
R + r R + r2 Θ + Z = kT =A
Separando T (t), tenemos:

T′
= A −→ T ′ − AkT = 0 (20)
kT
Nos queda:
R′′ 1 R′ 1 Θ′′ Z ′′
R + r R + r2 Θ + Z =A
Separando Z(z), tenemos:

Z ′′ R′′ 1 R′ 1 Θ′′
=A− − − 2 = B −→ Z ′′ − BZ = 0 (21)
Z R r R r Θ
Separando Θ(θ)
1 Θ′′ R′′ 1 R′ Θ′′ ′′ ′

r2 Θ =A−B− R − r R −→ Θ = (A − B)r2 − r2 RR − r RR = C

De donde se sigue:

Θ′′
= C −→ Θ − CΘ = 0 (22)
Θ
Separando R(r), tenemos:

r2 R′′ − rR′
(A − B)r2 − = C −→ r2 R′′ + rR′ + [(B − A)r2 + C]R = 0 (23)
R
Seleccionando
Xovee Xu 7

• A = −λ2
• B = µ2
• C = −ν 2

De esta forma las ecuaciones (20), (21), (22) y (23) nos quedan como:
2
T ′ + kλ2 T = 0 T (t) = C1 e−kλ t
(24)

Z ′′ − µ2 Z = 0 Z(z) = C2 cosh(µz) + C3 senh(µz) (25)

Θ′′ + ν 2 Θ = 0 Θ(θ) = C4 cos(νθ) + C5 sen(νθ) (26)

p p
r2 R′′ + rR′ + [(µ2 + λ2 )r2 − ν 2 ]R = 0 R(r) = C6 Jν ( µ2 + λ2 r) + C7 Yν ( µ2 + λ2 r) (27)

La solución entonces es:

u(r, θ, z, t) =
p p 2
[C6 Jν ( µ2 + λ2 r) + C7 Yν ( µ2 + λ2 r)] · [C4 cos(νθ) + C5 sen(νθ)] · [C2 cosh(µz) + C3 senh(µz)] · C1 · e−kλ t

Tarea VII. Resolver ecuación de Laplace por separación de variables.

La ecuación de Laplace es:

∇2(r,θ,z) u = 0 (28)

escribiendola en coordenadas cilindricas es:

∂2u 1 ∂u 1 ∂2u ∂2u

∂r 2 + r ∂r + r 2 ∂θ 2 + ∂z 2 =0

Sea u(r, θ, z) = R(r) · Θ(θ) · Z(z) la solución de la ecuación (28), calculando:

∂2u
• ∂r 2 = R′′ (r) · Θ(θ) · Z(z)

• ∂u
∂r = R′ (r) · Θ(θ) · Z(z)
∂2u
• ∂θ 2 = R(r) · Θ′′ (θ) · Z(z)
∂2u
• ∂z 2 = R(r) · Θ(θ) · Z ′′ (z)

Sustituyendo en (28), tenemos:

R′′ · Θ · Z + 1r R′ · Θ · Z + 1
r2 R · Θ′′ · Z + R · Θ · Z ′′ = 0

Dividiendo la ecuación anterior entre R · Θ · Z, tenemos:

R′′ 1 R′ 1 Θ′′ Z ′′
R + r R + r2 Θ + Z =0
Xovee Xu 8

Separando Z(z), tenemos:

Z ′′ R′′ 1 R′ 1 Θ′′
=− − − 2 = A −→ Z ′′ − AZ = 0 (29)
Z R r R r Θ
Separando Θ(θ)
1 Θ′′ ′′
1 R′ Θ′′ ′′ ′

r2 Θ = − RR − r R − A −→ Θ = −r2 RR − r RR − r2 A = B

De donde se sigue:

Θ′′
= B −→ Θ − BΘ = 0 (30)
Θ
Separando R(r), tenemos:

r2 R′′ + rR′
− − r2 A = B −→ r2 R′′ + rR′ + [Ar2 + B]R = 0 (31)
R
Seleccionando
• A = λ2
• B = −µ2
De esta forma las ecuaciones (29), (30) y (31) nos quedan como:

Z ′′ − λ2 Z = 0 Z(z) = C1 cosh(λz) + C2 senh(λz) (32)

Θ′′ + µ2 Θ = 0 Θ(θ) = C3 cos(µθ) + C4 sen(µθ) (33)

r2 R′′ + rR′ + [λ2 r2 − µ2 ]R = 0 R(r) = C5 Jµ (λr) + C6 Yµ (λr) (34)

La solución entonces es:

u(r, θ, z) = [C5 Jµ (λr) + C6 Yµ (λr)] · [C3 cos(µθ) + C4 sen(µθ)] · [C1 cosh(λz) + C2 senh(λz)]

Tarea VIII. Calcular la temperatura de estado estacionario de la siguiente figura.

Xovee Xu 9

Las condiciones de frontera de este problema son entonces:

u(c, θ, z) = 0 (35)

u(r, θ, 0) = 0 (36)

u(r, θ, a) = V0 (37)

La ecuación de disfusión es:

∂u
k∇2(r,θ,z) u = (38)
∂t
recordemos que cuando un fenómeno físico se encuentra en estado estacionario, no depende del tiempo, entonces la
ecuación (38) nos queda como:

∇2(r,θ,z) u = 0 (39)
que corresponde a la ecuación de Laplace. Escribiendola en coordenadas cilindricas es:

∂2u 1 ∂u 1 ∂2u ∂2u

∂r 2 + r ∂r + r 2 ∂θ 2 + ∂z 2 =0

Sea u(r, θ, z) = R(r) · Θ(θ) · Z(z) la solución de la ecuación (39), calculando:

∂2u
• ∂r 2 = R′′ (r) · Θ(θ) · Z(z)

• ∂u
∂r = R′ (r) · Θ(θ) · Z(z)
∂2u
• ∂θ 2 = R(r) · Θ′′ (θ) · Z(z)
∂2u
• ∂z 2 = R(r) · Θ(θ) · Z ′′ (z)

Sustituyendo en (19), tenemos:

R′′ · Θ · Z + 1r R′ · Θ · Z + 1
r2 R · Θ′′ · Z + R · Θ · Z ′′ = 0

Dividiendo la ecuación anterior entre R · Θ · Z, tenemos:

R′′ 1 R′ 1 Θ′′ Z ′′
R + r R + r2 Θ + Z =0

Separando Z(z), tenemos:

Z ′′ R′′ 1 R′ 1 Θ′′
=− − − 2 = A −→ Z ′′ − AZ = 0 (40)
Z R r R r Θ
Separando Θ(θ)

1 Θ′′ ′′
1 R′ Θ′′ ′′ ′

r2 Θ = − RR − r R − A −→ Θ = −r2 RR − r RR − r2 A = B
Xovee Xu 10

De donde se sigue:

Θ′′
= B −→ Θ − BΘ = 0 (41)
Θ
Separando R(r), tenemos:

r2 R′′ + rR′
− − r2 A = B −→ r2 R′′ + rR′ + [Ar2 + B]R = 0 (42)
R
Seleccionando

• A = λ2
• B = −µ2

De esta forma las ecuaciones (40), (41) y (42) nos quedan como:

Z ′′ − λ2 Z = 0 Z(z) = C1 cosh(λz) + C2 senh(λz) (43)

Θ′′ + µ2 Θ = 0 Θ(θ) = C3 cos(µθ) + C4 sen(µθ) (44)

r2 R′′ + rR′ + [λ2 r2 − µ2 ]R = 0 R(r) = C5 Jµ (λr) + C6 Yµ (λr) (45)

como la solución u debe estar acotada dentro de toda la circunferencia, en particular, también en el origen r = 0, la función
Yµ (λr) −→ ∞ cuando r −→ 0.

Para evitar esta catástrofe, elegimos C6 = 0, de esta manera la solución nos queda como:

R(r) = C5 Jµ (λr)

La solución entonces es:

u(r, θ, z) = C5 Jµ (λr) · [C3 cos(µθ) + C4 sen(µθ)] · [C1 cosh(λz) + C2 senh(λz)]

Aplicando la condición de frontera (35), tenemos:

u(c, θ, z) = C5 Jµ (λc) · [C3 cos(µθ) + C4 sen(µθ)] · [C1 cosh(λz) + C2 senh(λz)] = 0

O bien

u(c, θ, z) = R(c) · Θ(θ) · Z(z) = 0 −→ R(c) = 0 −→ C5 Jµ (λc) = 0

Aquí, C5 ̸= 0 pues de lo contrario tenemos la solución trivial u(r, θ, z) = 0, por tanto:

Jµ (λc) = 0

Si xµm = c · λµm son las raíces positivas de la ecuación anterior, entonces:

xµm
λµm = m = 1, 2, 3, · · · (eigenvalores) (46)
c
Además como
Xovee Xu 11

x2µm
A = λ2 Aµm = λ2µm = c2

En general, las raices de la ecuación 46 son diferentes para diferentes ordenes de µ.

De esta manera la solución nos queda:

u(r, θ, z) = C5 Jµ (λµm r) · [C3 cos(µθ) + C4 sen(µθ)] · [C1 cosh(λµm z) + C2 senh(λµm z)]

Hay una condición relacionada a Θ(θ) que no está escrita pero está inmersa, y es que u debe ser univaluada, esto quiere
decir que la función u debe tener el mismo valor para un punto (r, θ, z) y cualquier representación polar del mismo, esto es:

Θ(θ) = Θ(θ + 2π) (47)

De modo que Θ(θ) debe ser una función periódica y el orden µ debe ser un entero denotado con n con n = 0, 1, 2, · · ·
De este modo, la solución es:

u(r, θ, z) = C5 Jn (λnm r) · [C3 cos(nθ) + C4 sen(nθ)] · [C1 cosh(λnm z) + C2 senh(λnm z)]

Aplicando la condición de la ecuación (36), tenemos:

u(r, θ, 0) = C5 Jn (λnm r) · [C3 cos(nθ) + C4 sen(nθ)] · [C1 cosh(λnm · 0) + C2 senh(λnm · 0)] = 0

De donde se sigue que C1 = 0, de este modo:

u(r, θ, z) = C5 Jn (λnm r) · [C3 cos(nθ) + C4 sen(nθ)] · C2 senh(λnm z)

Haciendo Anm = C2 · C3 · C5 y Bnm = C2 · C4 · C5 y aplicando principio de superposición, tenemos:

∞ X
X ∞
u(r, θ, z) = [Anm cos(nθ) + Bnm sen(nθ)]Jn (λnm r)senh(λnm z) (48)
n=0 m=1

Aplicando la condición de frontera (37), tenemos:

P∞ P∞
u(r, θ, a) = n=0 m=1 [Anm cos(nθ) + Bnm sen(nθ)]Jn (λnm r)senh(λnm a) = V0

Agrupando de la siguiente manera, tenemos:

P∞ P∞
n=0 m=1 [Anm Jn (λnm r)senh(λnm a) · cos(nθ) + Bnm Jn (λnm r)senh(λnm a) · sen(nθ)] = V0

De donde se sigue:
P∞ P∞ P∞
n=0 [{ m=1 Anm Jn (λnm r)senh(λnm a)} · cos(nθ) + { m=1 Bnm Jn (λnm r)senh(λnm a)} · sen(nθ)] = V0

haciendo
∞
X
αn = Anm Jn (λnm r)senh(λnm a) (49)
m=1

∞
X
βn = Bnm Jn (λnm r)senh(λnm a) (50)
m=1

Haciendo las sustituciones correspondientes, nos queda:

Xovee Xu 12

P∞
n=0 [αn · cos(nθ) + βn · sen(nθ)] = V0

Para obtener los coeficientes de Fourier, hacemos:

Z 2π
1 2π
Z
2
αn = V0 · cos(nθ)dθ = V0 · cos(nθ)dθ (51)
2π 0 π 0
Z 2π Z 2π
2 1
βn = V0 · sen(nθ)dθ = V0 · sen(nθ)dθ (52)
2π 0 π 0

o bien

∞ Z 2π
X 1
Anm Jn (λnm r)senh(λnm a) = V0 · cos(nθ)dθ (53)
m=1
π 0

∞ Z 2π
X 1
Bnm Jn (λnm r)senh(λnm a) = V0 · sen(nθ)dθ (54)
m=1
π 0

Obteniendo el coeficiente n = 0, tenemos:

∞ Z 2π
X 1
A0m J0 (λ0m r)senh(λ0m a) = V0 · dθ (55)
m=1
π 0

Definimos:
Z 2π
1
f1 (θ) = V0 · cos(nθ)dθ (56)
π 0

Z 2π
1
f2 (θ) = V0 · sen(nθ)dθ (57)
π 0
y definimos los coeficientes

γnm = Anm senh(λnm a) (58)

ωnm = Bnm senh(λnm a) (59)

Entonces las escuaciones (53) y (54) quedan como:

∞
X
γnm Jn (λnm r) = f1 (θ) (60)
m=1

∞
X
ωnm Jn (λnm r) = f2 (θ) (61)
m=1

Para ortener los coeficientes de Fourier de las ecuaciones (60) y (61), se procede como:
Z c
2
γnm = r · Jn (λnm r) · f1 (θ)dr (62)
π[Jn+1 (λnm )]2 0
Z c
2
ωnm = r · Jn (λnm r) · f2 (θ)dr (63)
π[Jn+1 (λnm )]2 0
Xovee Xu 13

De donde se sigue:
Z c Z 2π
2 1
γnm = r · Jn (λnm r) · { V0 · cos(nθ) · dθ} · dr (64)
π[Jn+1 (λnm )]2 0 π 0

Z c Z 2π
2 1
ωnm = r · Jn (λnm r) · { V0 · sen(nθ) · dθ} · dr (65)
π[Jn+1 (λnm )]2 0 π 0

Luego por las ecuaciones (58) y (59), se sigue que:

Z c Z 2π
2 1
Anm = r · Jn (λnm r) · { V0 · cos(nθ) · dθ} · dr (66)
π[Jn+1 (λnm )]2 senh(λnm a) 0 π 0

Z c Z 2π
2 1
Bnm = r · Jn (λnm r) · { V0 · sen(nθ) · dθ} · dr (67)
π[Jn+1 (λnm )]2 senh(λnm a) 0 π 0

De esta manera sustituyendo las ecuaciones (66) y (67) en (48) tenemos la solución deseada.

Tarea IX. Resolver ecuación de onda en un tambor circular.

Las condiciones de frontera son:
u(a, θ, z, t) = 0 (68)

u(r, θ, 0, t) = 0 (69)

u(r, θ, b, t) = 0 (70)

u(r, θ, z, 0) = f (r, θ, z) (71)

ut (r, θ, z, 0) = 0 (72)

La ecuación de onda es:

∂2u
a2 ∇2(r,θ,z) u = (73)
∂t2
escribiendola en coordenadas cilindricas es:
2
1 ∂2u ∂2u ∂2u
a2 ( ∂∂ru2 + 1 ∂u
r ∂r + r 2 ∂θ 2 + ∂z 2 ) = ∂t2

Sea u(r, θ, z, t) = R(r) · Θ(θ) · Z(z) · T (t) la solución de la ecuación (73), calculando:

∂2u
• ∂r 2 = R′′ (r) · Θ(θ) · Z(z) · T (t)

• ∂u
∂r = R′ (r) · Θ(θ) · Z(z) · T (t)
∂2u
• ∂θ 2 = R(r) · Θ′′ (θ) · Z(z) · T (t)
∂2u
• ∂z 2 = R(r) · Θ(θ) · Z ′′ (z) · T (t)
Xovee Xu 14

∂2u
• ∂t2 = R(r) · Θ(θ) · Z(z) · T ′′ (t)

Sustituyendo en (73), tenemos:

a2 (R′′ · Θ · Z · T + 1r R′ · Θ · Z · T + 1
r2 R · Θ′′ · Z · T + R · Θ · Z ′′ · T ) = R · Θ · Z · T ′′

Dividiendo la ecuación anterior entre a2 R · Θ · Z · T , tenemos:

R′′ 1 R′ 1 Θ′′ Z ′′ T ′′
R + r R + r2 Θ + Z = a2 T =A

Separando T (t), tenemos:

T ′′
= A −→ T ′′ − Aa2 T = 0 (74)
a2 T
Nos queda:

R′′ 1 R′ 1 Θ′′ Z ′′
R + r R + r2 Θ + Z =A

Separando Z(z), tenemos:

Z ′′ R′′ 1 R′ 1 Θ′′
=A− − − 2 = B −→ Z ′′ − BZ = 0 (75)
Z R r R r Θ
Nos queda:

1 Θ′′ R′′ 1 R′ Θ′′ ′′ ′

r2 Θ =A−B− R − r R −→ Θ = (A − B)r2 − r2 RR − r RR = C

De donde se sigue:

Θ′′
= C −→ Θ − CΘ = 0 (76)
Θ
Separando R(r), tenemos:

r2 R′′ − rR′
(A − B)r2 − = C −→ r2 R′′ + rR′ + [(B − A)r2 + C]R = 0 (77)
R
Seleccionando

• A = −λ2
• B = −µ2

• C = −ν 2

De esta forma las ecuaciones (74), (75), (76) y (77) nos quedan como:

T ′′ + λ2 a2 T = 0 T (t) = C1 cos(λat) + C2 sen(λat) (78)

Z ′′ + µ2 Z = 0 Z(z) = C3 cos(µz) + C4 sen(µz) (79)

Θ′′ + ν 2 Θ = 0 Θ(θ) = C5 cos(νθ) + C6 sen(νθ) (80)

Xovee Xu 15

p p
r2 R′′ + rR′ + [(λ2 − µ2 )r2 − ν 2 ]R = 0 R(r) = C7 Jν ( λ2 − µ2 r) + C8 Yν ( λ2 − µ2 r) (81)

o bien
p
r2 R′′ + rR′ + [α2 r2 − ν 2 ]R = 0 α2 = λ2 − µ2 −→ α = λ2 − µ2 R(r) = C7 Jν (αr) + C8 Yν (αr) (82)

La solución entonces es:

u(r, θ, z, t) =
[C7 Jν (αr) + C8 Yν (αr)] · [C5 cos(νθ) + C6 sen(νθ)] · [C3 cos(µz) + C4 sen(µz)] · [C1 cos(λat) + C2 sen(λat)]

Ya que la solución u debe estar acotada en toda la circunferencia, en particular en r = 0, sabemos que Yν −→ ∞
cuando r −→ 0. Para evitar esta catástrofe, elegimos C8 = 0. De esta manera, la solución es:

u(r, θ, z, t) = C7 Jν (αr) · [C5 cos(νθ) + C6 sen(νθ)] · [C3 cos(µz) + C4 sen(µz)] · [C1 cos(λat) + C2 sen(λat)]

Comenzamos evaluando la primera condición en la frontera:

u(a, θ, z, t) = R(a) · Θ(θ) · Z(z) · T (t) = 0 −→ R(a) = 0 −→ C7 Jν (αa) = 0

Se observa que C7 ̸= 0, de lo contrario obtenemos la solución trivial y no podemos aplicar las condiciones de frontera
restantes. De este modo:
Jν (αa) = 0 (83)

Sean xmµν = a · αmµν las raíces positivas de la ecuación (83). De este modo:

xmµν
αmµν = m = 1, 2, 3, · · · (84)
a
x2mµν
Además α2 = λ2 − µ2 = a2 .

La solución entonces es:

u(r, θ, z, t) = C7 Jν (αmµν r) · [C5 cos(νθ) + C6 sen(νθ)] · [C3 cos(µz) + C4 sen(µz)] · [C1 cos(λat) + C2 sen(λat)]

Hay una condición de frontera relacionada a Θ(θ) que no está escrita pero está inmersa, y es que u debe ser univaluada,
esto quiere decir que la función u debe tener el mismo valor para un punto (r, θ, z) y cualquier representación polar del
mismo, esto es:
Θ(θ) = Θ(θ + 2π) (85)
De modo que Θ(θ) debe ser una función periodica y el orden ν debe ser un entero denotado con n con n = 0, 1, 2, · · · De
este modo:

u(r, θ, z, t) = C7 Jn (αmin r) · [C5 cos(nθ) + C6 sen(nθ)] · [C3 cos(µz) + C4 sen(µz)] · [C1 cos(λat) + C2 sen(λat)]

Aplicado la condición de frontera descrita por la ecuación (69), tenemos:

u(r, θ, 0, t) = C7 Jn (αmµn r) · [C5 cos(nθ) + C6 sen(nθ)] · Z(0) · [C1 cos(λat) + C2 sen(λat)] = 0

O sea,
Z(0) = C3 = 0 (86)

De modo que la solución se reduce a:

Xovee Xu 16

u(r, θ, z, t) = C7 · C4 Jn (αmµn r) · sen(µz) · [C5 cos(nθ) + C6 sen(nθ)] · [C1 cos(λat) + C2 sen(λat)]

Aplicando la coondición de frontera expresada por la ecuación (70), tenemos:

i·π
Z(b) = 0 −→ sen(µb) = 0 −→ µb = i · π −→ µ = i = 0, 1, 2, · · · (87)
b
De modo que la solución se reduce a:

u(r, θ, z, t) = C7 · C4 Jn (αmin r) · sen( i·π

b z) · [C5 cos(nθ) + C6 sen(nθ)] · [C1 cos(λat) + C2 sen(λat)]

Para determinar el valor de λ, recurrimos a la expresión:

q
i2 ·π 2 i2 ·π 2
α2 = λ2 − µ2 −→ λ2 = α2 + µ2 = αmin
2
+ b2 −→ λ = 2
αmin + b2

De esta manera, la solución nos queda como:

u(r, θ, z, t) = q q
i2 ·π 2
C7 · C4 Jn (αmin r) · sen( i·π
b z) · [C5 cos(nθ) + C6 sen(nθ)] · [C1 cos(
2
αmin + i·π
b at) + C2 sen(
2
αmin + b2 at)]

Derivando la solución y aplicando la última condición de frontera, tenemos:

q
i2 ·π 2
ut (r, θ, z, t) = C7 · C4 Jn (αmin r) · sen( i·π z) · [C5 cos(nθ) + C 6 sen(nθ)] · [−C1 sen( 2
αmin + b2 at) +
q b q
2 2 2 2
C2 cos( αmin 2 + i b·π
2 at)](
2
αmin + i b·π
2 a)

Luego:
q
ut (r, θ, z, 0) = C7 · C4 Jn (αmin r) · sen( i·π
b z) · [C 5 cos(nθ) + C 6 sen(nθ)] · ( 2
αmin + i·π
b a)C2 =0

De donde se sigue que C2 = 0, de este modo la solución nos queda como:

q
i·π
2
u(r, θ, z, t) = C7 · C4 · C1 Jn (αmin r) · cos( αmin + b at) · sen( i·π
b z) · [C5 cos(nθ) + C6 sen(nθ)]

Haciendo Anmi = C7 · C4 · C1 · C5 y Bnmi = C7 · C4 · C1 · C6 tenemos la solución:

q
i·π
2
u(r, θ, z, t) = [Anmi cos(nθ) + Bnmi sen(nθ)]Jn (αmin r) · cos( αmin + b at) · sen( i·π
b z)

Aplicando principio de superposición y la condicion de frontera temporal, tenemos:

P∞ P∞ P∞
u(r, θ, z, 0) = n=0 m=1 i=0 [Anmi cos(nθ) + Bnmi sen(nθ)]Jn (αmin r) · sen( i·π
b z) = f (r, θ, z)

Agrupando, tenemos:
P∞ P∞ P∞ P∞
u(r, θ, z, 0) = n=0 m=1 [{ i=0 Anmi Jn (αmin r) · sen( i·π
b z)}cos(nθ) + { i=0 Bnmi Jn (αmin r) ·
sen( i·π
b z)}sen(nθ)] = f (r, θ, z)

Definimos:
P∞
γnm = i=0 Anmi sen( i·π
b z)
Xovee Xu 17

P∞
δnm = i=0 Bnmi · sen( i·π
b z)

De este modo,sustituyendo, tenemos:

P∞ P∞
u(r, θ, z, 0) = n=0 m=1 [γnm Jn (αmin r) · cos(nθ) + δnm Jn (αmin r) · sen(nθ)] = f (r, θ, z)

Agrupando, tenemos:
P∞ P∞ P∞
u(r, θ, z, 0) = n=0 [{ m=1 γnm Jn (αmin r)} · cos(nθ) + { m=1 δnm Jn (αmin r)} · sen(nθ)] = f (r, θ, z)

Definimos:
P∞
αn = m=1 γnm Jn (αmin r)
P∞
βn = m=1 δnm Jn (αmin r)

Sustituyendo, tenemos:
P∞
u(r, θ, z, 0) = n=0 [αn · cos(nθ) + βn · sen(nθ)] = f (r, θ, z)

Obteniendo los coeficientes de Bessel-Fourier, tenemos:

2
R 2π 1
R 2π
αn = 2π 0
f (r, θ, z)cos(nθ)dθ = π 0
f (r, θ, z)cos(nθ)dθ

2
R 2π 1
R 2π
βn = 2π 0
f (r, θ, z)sen(nθ)dθ = π 0
f (r, θ, z)sen(nθ)dθ

De esta forma, se sigue que:

2
Ra 2
R a R 2π
γnm = 2
a2 Jn+1 (αmin a) 0
rαn Jn (αmin r)dr = 2
a2 πJn+1 (αmin a) 0 0
rJn (αmin r)f (r, θ, z)cos(nθ)drdθ

2
Ra 2
R a R 2π
δnm = 2
a2 Jn+1 (αmin a) 0
rβn Jn (αmin r)dr = 2
a2 πJn+1 (αmin a) 0 0
rJn (αmin r)f (r, θ, z)sen(nθ)drdθ

Por último para obtener los coeficientes Anmi y Bnmi , hacemos:

2
Rb
Anmi = b 0
γnm sen( i·π
b z)dz

2
Rb
Bnmi = b 0
δnm sen( i·π
b z)dz

De donde se sigue:

4
R b R a R 2π
Anmi = 2
a2 πbJn+1 (αmin a) 0 0 0
rJn (αmin r)f (r, θ, z)cos(nθ)sen( i·π
b z)dθdrdz

4
R b R a R 2π
Bnmi = 2
a2 πbJn+1 (αmin a) 0 0 0
rJn (αmin r)f (r, θ, z)sen(nθ)sen( i·π
b z)dθdrdz

Al sustituir estos coeficientes en la solución obtenemos la solución general con condiciones de frontera.
Xovee Xu 18

Tarea X. Resolver ecuación de onda en una membrana circular.

La ecuación de onda para una membrana circular es:

∂ 2 u 1 ∂u 1 ∂2u ∂2u
a2 ( + + ) = (88)
∂r2 r ∂r r2 ∂θ2 ∂t2
Con condiciones en la frontera:
u(a, θ, t) = 0 (89)

u(r, θ, 0) = f (r, θ) (90)

ut (r, θ, 0) = 0 (91)

Para resolver la ecuación (88) procedemos por el método de separación de variables. Sea u(r, θ, t) = R(r) · Θ(θ) · T (t)
la solución de la ecuación, calculando:
∂2u
• ∂r 2 = R′′ (r) · Θ(θ) · T (t)

• ∂u
∂r = R′ (r) · Θ(θ) · T (t)
∂2u
• ∂θ 2 = R(r) · Θ′′ (θ) · T (t)
∂2u
• ∂t2 = R(r) · Θ(θ) · T ′′ (t)

Sustituyendo en (88) y dividiendo entre a2 · R(r) · Θ(θ) · T (t) tenemos:

R′′ 1 R′ 1 Θ′′ T ′′
+ + 2 = 2 (92)
R r R r Θ a T
Separando T (t), tenemos:

R′′ 1 R′ 1 Θ′′ T ′′
+ + 2 = 2 = A −→ T ′′ − Aa2 T = 0 (93)
R r R r Θ a T
Nos queda:

R′′ 1 R′ 1 Θ′′
R + r R + r 2 Θ =A

Separando Θ(θ), tenemos:

1 Θ′′ R′′ 1 R′ Θ′′ R′′ R′

2
=− − + A −→ = −r2 −r + Ar2 = B −→ Θ′′ − BΘ = 0 (94)
r Θ R r R Θ R R
Nos queda:
′′ ′ 2
R′′ +rR′
−r2 RR − r RR + Ar2 = B −→ − r R + Ar2 = B

Separando R(r), tenemos:

r2 R′′ + rR′ = (Ar2 − B)R −→ r2 R′′ + rR′ + (−Ar2 + B)R = 0 (95)

Haciendo las elecciones

Xovee Xu 19

• A = −λ2
• B = −µ2

Las ecuaciones (93), (94) y (95), nos quedan como:

T ′′ + λ2 a2 T = 0 −→ T (t) = C1 cos(λat) + C2 sen(λat) (96)

Θ′′ + µ2 Θ = 0 −→ Θ(θ) = C3 cos(µθ) + C4 sen(µθ) (97)

r2 R′′ + rR′ + (λ2 r2 − µ2 )R = 0 −→ R(r) = C5 Jµ (λr) + C6 Yµ (λr) (98)

Dado que u debe estar acotada en toda la circunferencia, incluyendo r = 0, sabemos que Yµ (λr) −→ ∞ cuando
r −→ 0. Para evitar esta catástrofe, elegimos C6 = 0. De este modo la solución general del problema es:

u(r, θ, t) = C5 Jµ (λr)[C3 cos(µθ) + C4 sen(µθ)][C1 cos(λat) + C2 sen(λat)] (99)

Aplicando la primera condicón de frontera, tenemos que:

u(a, θ, t) = R(a) · Θ(θ) · T (t) = 0 −→ R(a) = 0 −→ C5 Jµ (λa) = 0

De aquí C5 ̸= 0, de lo contrario obtenemos la solución trivial u(r, θ, t) = 0 y ya no se satisfacen las demás condiciones de
frontera. De esta manera:

Jµ (λa) = 0 (100)

Sean xmµ = λmµ a los ceros positivos de la ecuación anterior, de este modo:

x2mµ
Amµ = −λ2mµ = m = 1, 2, 3, · · · (101)
a2
En general, las raíces de (100) son diferentes para diferentes valores del orden de µ.

De este modo la solución queda como:

u(r, θ, t) = C5 Jµ (λmµ r)[C3 cos(µθ) + C4 sen(µθ)][C1 cos(λmµ at) + C2 sen(λmµ at)] (102)

Hay una condición implícita en el problema y es que u debe ser continua y acotada, por ende univaluada. Esto quiere
decir que el valor de u debe ser el mismo para un punto (r, θ) y cualquier representación polar de este punto. esto es

Θ(θ) = Θ(θ + 2π) (103)

de este modo la solución de Θ(θ) debe ser periodica, esto se logra si:

µ=n n = 0, 1, 2, . . .

Por tanto, la solución queda como:

u(r, θ, t) = C5 Jn (λmn r)[C3 cos(nθ) + C4 sen(nθ)][C1 cos(λmn at) + C2 sen(λmn at)] (104)

Aplicando la última condición de frontera, tenemos:

ut (r, θ, t) = C5 Jn (λmn r)[C3 cos(nθ) + C4 sen(nθ)][−C1 sen(λmn at) + C2 cos(λmn at)](λmn a) = 0

Xovee Xu 20

De donde se sigue:

ut (r, θ, 0) = C5 Jn (λmn r)[C3 cos(nθ) + C4 sen(nθ)](λmn a)C2 = 0 −→ C2 = 0

la solución queda como:

u(r, θ, t) = C5 · C1 · Jn (λmn r)cos(λmn at)[C3 cos(nθ) + C4 sen(nθ)] (105)

Haciendo Amn = C1 · C3 · C5 y Bmn = C1 · C4 · C5 y aplicando principio de superposición, tenemos:

∞ X
X ∞
u(r, θ, t) = [Amn cos(nθ) + Bmn sen(nθ)]Jn (λmn r)cos(λmn at) (106)
n=0 m=1

Aplicando la condición de frontera

P∞ P∞
u(r, θ, 0) = n=0 m=1 [Amn cos(nθ) + Bmn sen(nθ)]Jn (λmn r) = f (r, θ)

Agrupando como se sigue, tenemos:

P∞ P∞ P∞
u(r, θ, 0) = n=0 [{ m=1 Amn Jn (λmn r)}cos(nθ) + { m=1 Bmn Jn (λmn r)}sen(nθ)] = f (r, θ)

Definamos:
∞
X
αn = Amn Jn (λmn r) (107)
m=1

∞
X
βn = Bmn Jn (λmn r) (108)
m=1

Sustituyendo, tenemos:
P∞
u(r, θ, 0) = n=0 [αn cos(nθ) + βn sen(nθ)] = f (r, θ)

Para obtener los coeficientes de Fourier, hacemos:

Z 2π Z 2π
2 1
αn = f (r, θ)cos(nθ) · dθ = f (r, θ)cos(nθ) · dθ (109)
2π 0 π 0

Z 2π Z 2π
2 1
βn = f (r, θ)sen(nθ) · dθ = f (r, θ)sen(nθ) · dθ (110)
2π 0 π 0

De esta manera, para obtener los coeficientes Amn y Bmn , hacemos

Z a Z a Z 2π
2 2
Amn = 2 rJn (λnm )αn dr = 2 rJn (λnm )f (r, θ)cos(nθ) · dθdr (111)
Jn+1 (λmn ) 0 πJn+1 (λmn ) 0 0

Z a Z a Z 2π
2 2
Bmn = 2 rJn (λnm )βn dr = 2 rJn (λnm )f (r, θ)sen(nθ) · dθdr (112)
Jn+1 (λmn ) 0 πJn+1 (λmn ) 0 0

Sustituyendo estos coeficientes en la solución, obtenemos la función deseada.

Xovee Xu 21

Tarea XI. Resolver ecuación de difusión en una placa circular.

Tenemos una placa circular metálica con las siguientes condiciones en la frontera:

u(a, θ, t) = 0 (113)

u(r, θ, 0) = U0 (114)

La ecuación de difusión para una placa metálica circular

∂ 2 u 1 ∂u 1 ∂2u ∂u
k( 2
+ + 2 2
)= (115)
∂r r ∂r r ∂θ ∂t

Sea u(r, θ, t) = R(r) · Θ(θ) · T (t) la solución de (115), calculando:

∂2u
∂r 2 = R′′ (r) · Θ(θ) · T (t)

∂u
∂r = R′ (r) · Θ(θ) · T (t)

∂2u
∂θ 2 = R(r) · Θ′′ (θ) · T (t)

∂u
∂t = R(r) · Θ(θ) · T ′ (t)

Sustituyendo y diviendo entre k · R(r) · Θ(θ) · T (t), tenemos:

R′′ (r) 1 R′ (r) 1 Θ′′ (θ) T ′ (t)

+ + 2 = (116)
R(r) r R(r) r Θ(θ) kT (t)
Separando T (t) tenemos:

T ′ (t)
= A −→ T ′ − AkT = 0 (117)
kT (t)
Nos queda:

R′′ (r) ′
1 R (r)
′′
1 Θ (θ)
R(r) + r R(r) + r 2 Θ(θ) =A

Separando Θ(θ), tenemos:

1 Θ′′ (θ) R′′ (r) 1 R′ (r) Θ′′ (θ) ′′

2 R (r) R′ (r)
= A − − −→ = Ar 2
− r − r = B −→ Θ′′ − BΘ = 0 (118)
r2 Θ(θ) R(r) r R(r) Θ(θ) R(r) R(r)

Nos queda:
′′ ′
r 2 R′′ +rR′
Ar2 − r2 RR(r)
(r)
− r RR(r)
(r)
= B −→ Ar2 − R =B

Separando R(r)

r2 R′′ + rR′
= Ar2 − B −→ r2 R′′ + rR′ = (Ar2 − B)R −→ r2 R′′ + rR′ + (−Ar2 + B)R = 0 (119)
R

Haciendo las elecciones

Xovee Xu 22

• A = −λ2
• B = −µ2

De esta manera las ecuaciones (117), (118) y (119) nos quedan:

2
T ′ + λ2 kT = 0 T (t) = C1 · e−λ kt
(120)

Θ′′ + µ2 Θ = 0 Θ(θ) = C2 cos(µθ) + C3 sen(µθ) (121)

r2 R′′ + rR′ + (λ2 r2 − µ2 )R = 0 R(r) = C4 Jµ (λr) + C5 Yµ (λr) (122)

Sabemos que u es acotada en toda la circunferencia, particularmente en r = 0. Luego, sabemos que Yµ −→ ∞ cuando
r −→ 0. Para evitar esta catástrofe, elegimos C5 = 0. De esta manera, la solución es:
2
u(r, θ, t) = C1 · C4 · e−λ kt
Jµ (λr)[C2 cos(µθ) + C3 sen(µθ)] (123)

Haciendo A = C1 · C2 · C4 y B = C1 · C3 · C4 , nos queda:

2
u(r, θ, t) = e−λ kt
Jµ (λr)[A · cos(µθ) + B · sen(µθ)] (124)

Aplicando la primera condición en la frontera, tenemos:

u(r, θ, t) = R(a) · Θ(θ) · T (t) = 0 −→ R(a) = 0 −→ Jµ (λa) = 0

λ λ λ λ
Sean a1µ , a2µ , a3µ · · · con m = 1, 2, · · · , en general mµ
a las raíces positivas de la ecuación anterior. Además sabemos
que las raícen son distintas si cambian el orden. De esta manera, la solución queda como:

λmµ 2 λmµ r
u(r, θ, t) = e−( a ) kt Jµ ( )[A · cos(µθ) + B · sen(µθ)] (125)
a
Hay una condición implícita en Θ(θ) y es que se sabe que u es continua, acotada y debe ser univaluada. Esto lo que nos
quiere decir es que la temperatura en un punto (r, θ) debe ser la misma para todas las representanciones polares del punto.
Es decir.
Θ(θ) = Θ(θ + 2π) (126)

De esta manera, las soluciones de Θ(θ) deben ser periodicas y el orden µ = n con n = 0, 1, 2, · · · .
De esta manera, la solución se reduce a:

λ2
mn kt ) r
u(r, θ, t) = e−( a2 Jn (λmµ )[A · cos(nθ) + B · sen(nθ)] (127)
a
De esta manera aplicando el principio de superposición, tenemos:

∞ X
∞ λ2
X mµ kt r
u(r, θ, t) = [Anm · cos(nθ) + Bnm · sen(nθ)]e−( a2
)
Jn (λmn ) (128)
n=0 m1
a

Aplicando la segunda condición en la frontera, tenemos:

∞ X∞
X r
u(r, θ, 0) = [Anm · cos(nθ) + Bnm · sen(nθ)]Jn (λmn ) = U0 (129)
n=0 m
a
1
Xovee Xu 23

Agrupando de la siguiente manera:

∞ ∞ ∞
X X r X r
u(r, θ, 0) = [{ Anm Jn (λmn )} · cos(nθ) + { Bnm Jn (λmn )} · sen(nθ)] = U0 (130)
n=0 m
a m
a
1 1

Hacemos los cambios:

∞
X r
αn = Anm Jn (λmn ) (131)
m1
a

∞
X r
βn = Bnm Jn (λmn ) (132)
m1
a

De esta manera, la solución nos queda como:

∞
X
u(r, θ, 0) = [αn · cos(nθ) + βn · sen(nθ)] = U0 (133)
n=0

Calculando los coeficientes de Fourier, tenemos:

Z 2π
1 2π
Z
2
αn = U0 · cos(nθ) · dθ = U0 · cos(nθ) · dθ (134)
2π 0 π 0

Z 2π Z 2π
2 1
βn = U0 · sen(nθ) · dθ = U0 · sen(nθ) · dθ (135)
2π 0 π 0

De esta manera, para determinar los factores Anm y Bnm , hacemos uso de las relaciones de ortogonalidad en las series
de Fourier-Bessel, como se muestra a continuación:
Z a
2 r
Anm = 2 2 r · αn · Jn (λmn )dr (136)
a Jn+1 (λmn a) 0 a
Z a
2 r
Bnm = 2 2 r · βn · Jn (λmn )dr (137)
a Jn+1 (λmn a) 0 a

Sustituyendo αn y βn , nos quedan entonces:

Z a Z 2π
2 r
Anm = 2 U0 · r · Jn (λmn ) · cos(nθ)drdθ (138)
a2 πJn+1 (λmn a) 0 0 a
Z a Z 2π
2 r
Bnm = 2 U0 · r · Jn (λmn ) · sen(nθ)drdθ (139)
a2 πJn+1 (λmn a) 0 0 a

Sustituyendo estos coeficientes en la solución general, obetenemos la solución al problema con condiciones de frontera.