Modelamiento matemático y simulación avanzada
Segunda sesión
William González Calderón1
1 UniversidadAutónoma de Bucaramanga (UNAB)
Facultad de ingenierías
Maestría en ingeniería en energía
Colombia, Bucaramanga, 2021
Colombia, Bucaramanga, 2021 1/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Contenidos
1 Sección anterior
2 Condiciones de frontera
3 Elementos finitos
4 Más métodos de solución
5 Métodos de elementos finitos
Colombia, Bucaramanga, 2021 2/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación del calor
Sección anterior
Modelo matemático de flujo calor sin fuentes en un alambre uniforme cuyos
extremos se mantienen a temperatura constante cero:
∂u ∂2u
(x, t) = β 2 (x, t), 0 < x < L, t > 0, (1)
∂t ∂x
u(0, t) = u(L, t) = 0, t ≥ 0, (2)
u(x, 0) = f (x), 0 < x < L. (3)
Solución,
∞
X nπ 2 nπx
u(x, t) = cn e −β( L
) t
sin ( ) (4)
L
n=1
Colombia, Bucaramanga, 2021 3/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Condición inicial de la Ecuación de calor
Una solución u de la forma (4) que satisface (3) debe cumplir
∞
X nπx
u(x, 0) = cn sin ( ) = f (x), 0<x <L (5)
L
n=1
Luego, el problema de flujo de calor en un alambre delgado se reduce al
problema determinar un desarrollo de f(x) de la forma:
∞
X nπx
f (x) = cn sin ( ), 0<x <L (6)
L
n=1
Tal desarrollo es una serie de Fourier.
Colombia, Bucaramanga, 2021 4/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Coeficientes de Fourier
Si se puede elegir los coeficientes cn de modo que satisfaga la condición
(6), entonces el desarrollo de u(x, t) en (4) es una solución formal del
problema de flujo de calor. Los coeficientes vienen dados para series de
Fourier de funciones senos
2 L
Z
cn = f (x) sin(nπx/L)dx,
L 0
Entonces, u(x, t) quedaría
∞
X nπ 2 nπx
u(x, t) = cn e −β( L
) t
sin ( )
L
n=1
Si este desarrollo converge a una función con segundas derivadas parciales
continuas, entonces la solución formal es la solución única del problema.
Colombia, Bucaramanga, 2021 5/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de Onda
Sección anterior
Modelo matemático de la deflexión u(x, t) de una cuerda en cualquier
punto x y en cualquier tiempo t.
∂2u 2
2∂ u
(x, t) = α (x, t), 0 < x < L, t > 0,
∂t 2 ∂x 2
u(0, t) = u(L, t) = 0, t ≥ 0,
u(x, 0) = f (x), 0 ≤ x ≤ L,
∂u
(x, 0) = g (x), 0 ≤ x ≤ L.
∂t
Solución,
∞ h
X nπα nπα i nπx
u(x, t) = an cos ( t) + bn sin ( t) sin ( )
L L L
n=1
Se determinaron los coeficientes an y bn
Colombia, Bucaramanga, 2021 6/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de Laplace1
Sección anterior
∂2u ∂2u
(x, y ) + (x, y ) = 0, 0 < x < a, 0 < y < b,
∂x 2 ∂y 2
1
[3]
Colombia, Bucaramanga, 2021 7/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
∂2u ∂2u
Ecuación de Laplace ∂x 2 + ∂y 2 =0
Aplicaciones
Surge en el caso estacionario de la ecuación del calor. Luego, sólo
importa especificar las condiciones de frontera.
Analizar el desplazamiento estático de un membrana estirada y
asegurada en el espacio a lo largo de la frontera de una región: u debe
satisfacer la ecuación de Laplace dentro de la región.
Analizar los potenciales electrostático y gravitacional en ciertos
campos de fuerzas: u debe satisfacer la ecuación de Laplace en
cualquier región libre de cargas eléctricas o de masa;
En Mecánica de fluidos de un fluido idealizado, la función de flujo
cuyas curvas de nivel u(x, y ) = kte, representan la trayectoria de
partículas en el fluido: de nuevo, u debe satisfacer la ecuación de
Laplace en la región de flujo.
Colombia, Bucaramanga, 2021 8/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de Laplace
Sección anterior
∂2u ∂2u
(x, y ) + (x, y ) = 0, 0 < x < a, 0 < y < b,
∂x 2 ∂y 2
u(0, y ) = u(a, y ) = u(x, b) = 0, 0 ≤ x ≤ a, 0 ≤ y ≤ b,
u(x, 0) = f (x), 0 ≤ x ≤ a,
Solución,
∞
X nπx nπ
u(x, y ) = cn sin ( ) sinh ( (b − y ))
a a
n=1
Se demuestra que
a
2
Z
nπx
cn = f (x) sin dx
a sinh (nπb/a) 0 a
Colombia, Bucaramanga, 2021 9/
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras aisladas
Se aborda el problema de flujo de calor cuando los extremos están aislados,
es decir, se mantiene a una temperatura constante no nula. Lo anterior,
implica condiciones de frontera no homogéneas en la frontera.
∂u ∂2u
(x, t) = β 2 (x, t), 0 < x < L, t > 0, (7)
∂t ∂x
∂u ∂u
(0, t) = (L, t) = 0, t > 0, (8)
∂x ∂x
u(x, 0) = f (x), 0 < x < L. (9)
Decir aislados significa que no fluye calor hacia afuera (o hacia adentro)
por los extremos del alambre. El principio de conducción del calor implica
que el gradiente de temperatura debe anularse en estos extremos (Ver
condición (8)).
Colombia, Bucaramanga, 2021 10 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras aisladas
Solución del problema de flujo de calor con extremos aislados,
∞
a0 X 2 nπx
u(x, t) = + an e −β(nπ/L) t cos .
2 L
n=1
Si se sustituye la condición inicial, se obtiene
∞
a0 X nπx
u(x, 0) = + an cos = f (x), 0 < x < L.
2 L
n=1
Esto significa que si eligen los coeficientes an en la serie de cosenos de
Fourier para f ,
L
2
Z
nπx
an = f (x) cos dx, n = 0, 1, 2, . . .
L 0 L
Colombia, Bucaramanga, 2021 11 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras no homogéneas
Cuando los extremos del alambre se mantienen a 0o C o cuando están
aislados, las condiciones de frontera son homogéneas. Pero, si los extremos
del alambre se mantienen a temperaturas constantes distintas de cero:
∂u ∂2u
(x, t) = β 2 (x, t), 0 < x < L, t > 0, (10)
∂t ∂x
u(0, t) = U1 , u(L, t) = U2 , t > 0, (11)
u(x, 0) = f (x), 0 < x < L. (12)
Entonces, las condiciones en la frontera son no homogéneas.
Colombia, Bucaramanga, 2021 12 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras no homogéneas
Similar a caso de vibraciones forzadas, se espera que una solución al
problema de flujo de calor con condiciones no homogéneas en la frontera,
conste de una solución de estado estacionario v (x) que satisfaga las
condiciones no homogéneas en la frontera (11), más una solución
transitoria w (x, t),
u(x, t) = v (x) + w (x, t),
donde w (x, t) y sus derivadas parciales tienden a cero cuando t → ∞.
Entonces, la función w (x, t) debe satisfacer las condiciones homogéneas en
la frontera.
Colombia, Bucaramanga, 2021 13 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras no homogéneas
Al sustituir u(x, t) = v (x) + w (x, t) en el problema de flujo de calor con
fronteras no homogéneas,
∂u ∂w ∂2w
(x, t) = (x, t) = βv 00 (x) + β 2 (x, t), 0 < x < L, t > 0, (13)
∂t ∂t ∂x
v (0) + w (0, t) = U1 , v (L) + w (L, t) = U2 , t > 0, (14)
v (x) + w (x, 0) = f (x), 0 < x < L. (15)
Al hacer t → ∞ en (13) y (14), luego la solución transitoria w (x, t) → 0,
luego para el estado estacionario, se obtiene
v 00 (x) = 0, 0 < x < L,
v (0) = U1 , v (L) = U2 .
Colombia, Bucaramanga, 2021 14 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras no homogéneas
se obtiene el siguiente problema
v 00 (x) = 0, 0 < x < L,
v (0) = U1 , v (L) = U2 .
cuya solución del estado estacionario está dada por
(U1 − U2 )x
v (x) = U1 +
L
Colombia, Bucaramanga, 2021 15 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras no homogéneas
(U1 −U2 )x
Al reemplazar v (x) = U1 + L , el problema anterior se reduce a:
∂w ∂2w
(x, t) = β 2 (x, t), 0 < x < L, t > 0, (16)
∂t ∂x
w (0, t) = w (L, t) = 0, t > 0, (17)
(U1 − U2 )x
w (x, 0) = f (x) − U1 − , 0 < x < L. (18)
L
Colombia, Bucaramanga, 2021 16 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fronteras no homogéneas
Donde la solución del estado transitorio está dada por,
∞
X 2 nπx
w (x, t) = an e −β(nπ/L) t sin .
L
n=1
y la solución general,
∞
(U1 − U2 )x X 2 nπx
u(x, t) = U1 + + an e −β(nπ/L) t sin .
L L
n=1
donde
L
2 (U1 − U2 )x
Z
nπx
an = f (x) − U1 − sin dx.
L 0 L L
Colombia, Bucaramanga, 2021 17 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fuentes de calor
Ejercicio de clase
Suponga el problema de flujo de calor cuando se tiene la influencia de una
fuente de calor independiente del tiempo,
∂u ∂2u
(x, t) = β 2 (x, t) + P(x), 0 < x < L, t > 0, (19)
∂t ∂x
u(0, t) = U1 , u(L, t) = U2 , t > 0, (20)
u(x, 0) = f (x), 0 < x < L. (21)
De manera similar, se considera una solución en estado estacionario v (x) y
una solución transitoria w (x, t), cuyas derivadas parciales tienden a cero
cuanto t → ∞, es decir,
u(x, t) = v (x) + w (x, t).
Colombia, Bucaramanga, 2021 18 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fuentes de calor
Al sustituir u(x, t) = v (x) + w (x, t) en el problema de flujo de calor,
∂w ∂2w
= βv 00 (x) + β 2 (x, t) + P(x), 0 < x < L, t > 0, (22)
∂t ∂x
v (0) + w (0, t) = U1 , v (L) + w (L, t) = U2 , t > 0, (23)
v (x) + w (x, 0) = f (x), 0 < x < L. (24)
Colombia, Bucaramanga, 2021 19 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fuentes de calor
Al hacer t → ∞ en (22) y (23), luego la solución transitoria w (x, t) → 0,
luego para el estado estacionario, se obtiene
1
v 00 (x) = − P(x), 0 < x < L,
β
v (0) = U1 , v (L) = U2 .
cuya solución del estado estacionario está dada por
Z L Z z Z x Z z
1 x 1
v (x) = U2 − U1 + P(s)ds dz +U1 − P(s)ds dz
0 0 β L 0 0 β
Colombia, Bucaramanga, 2021 20 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fuentes de calor
Al reemplazar v (x), el problema anterior quedaría:
∂w ∂2w
(x, t) = β 2 (x, t), 0 < x < L, t > 0, (25)
∂t ∂x
w (0, t) = w (L, t) = 0, t > 0, (26)
w (x, 0) = f (x) − v (x), 0 < x < L. (27)
Colombia, Bucaramanga, 2021 21 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de calor con fuentes de calor
Donde la solución del estado transitorio está dada por,
∞
X 2 nπx
w (x, t) = an e −β(nπ/L) t sin .
L
n=1
y la solución general,
u(x, t) = v (x) + w (x, t)
Colombia, Bucaramanga, 2021 22 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Condiciones de frontera
Condición de Dirichlet
Donde la solución u(x, y ) de la ecuación en un dominio D debe satisfacer
u(x, y ) = f (x, y ) en la frontera de D
Condición de Neumann
Donde además se le pide que la derivada direccional
∂u
(x, y ) = g (x, y ), en la frontera de D
∂n
donde g (x, y ) una función definida sobre la frontera en D.
Colombia, Bucaramanga, 2021 23 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Condiciones de frontera
Condición de Robin
∂u
(x, y ) + hu(x, y ) = g2 (x, y ), en la frontera de D
∂n
Condiciones mixtas
Decimos que las condiciones en la frontera están mezcladas si la solución
debe satisfacer u(x, y ) = f (x, y ) en una parte de la frontera y
∂u
∂n (x, y ) = g1 (x, y ) en la parte restante de la frontera.
Colombia, Bucaramanga, 2021 24 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Condiciones de frontera
Ejemplos
u(L, t) = u0 establece que la frontera x = L se conserva de alguna
manera a temperatura constante u0 en todo tiempo t > 0.
∂u
∂x |x=L = 0, significa que la frontera en x = L está aislada.
A partir de la ley empírica para la transferencia de calor, el flujo de
calor a través de la frontera (esto es, la cantidad de calor por unidad
de área por unidad de tiempo que se transfiere a través de la frontera)
es proporcional al valor de la derivada normal ∂u/∂n de la
temperatura u.
Colombia, Bucaramanga, 2021 25 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Condiciones de frontera
Ejemplo
∂u
La expresión ∂x |x=L = h(um − u(L, t)) donde h > 0 y um son constantes.
Se interpreta como el calor que se "pierde" en el extremo derecho de la
varilla por estar en contacto con un medio, como aire o agua, que se
encuentra a temperatura constante um .
A partir de la ley de Newton para el enfriamiento sabemos que el flujo de
calor hacia afuera de la varilla es proporcional a la diferencia entre la
temperatura u(L, t) en la frontera y la temperatura del medio um del medio
circundante.
Colombia, Bucaramanga, 2021 26 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Condiciones de frontera
Ejemplo
Observe que si existe pérdida de calor desde el extremo izquierdo de la
varilla, la condición de frontera es
∂u
|x=0 = −h(um − u(0, t))
∂x
donde h > 0 y um son constantes.
El cambio de signo es consistente con el supuesto de que la varilla se
encuentra a una temperatura más elevada que el medio circundante en los
extremos, por lo que u(0, t) > um y u(L, t) > um . En x = 0 y x = L, las
pendientes ux (0, t) y ux (L, t) deben ser positiva y negativa
respectivamente.
Colombia, Bucaramanga, 2021 27 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Variaciones en la ecuación diferencial
Casos
Si existe transferencia de calor desde la superficie lateral de una varilla
hacia un medio circundante a temperatura constante um , entonces la
ecuación de calor es
∂u ∂2u
= k 2 − h(um − u)
∂t ∂x
Colombia, Bucaramanga, 2021 28 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Variaciones en la ecuación diferencial
Casos
En la ecuación de onda, cuando las fuerzas externas de amortiguamiento y
restablecimiento elástico se tienen en cuenta, toma la forma
∂2u 2
2∂ u ∂u
2
= α 2
− f (x, t) − c − ku,
∂t ∂x ∂t
donde f es una fuerza externa, c constante de amortiguamiento y k
constante de resistencia o restablecimiento de la cuerda.
También, se puede considerar un término de forzamiento h(x, t)
dependiente del tiempo
∂2u 2
2∂ u
= α + h(x, t)
∂t 2 ∂x 2
Colombia, Bucaramanga, 2021 29 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación de Laplace con un problema mixto en la frontera2
∂2u ∂2u
+ = 0, 0 < x < a, 0 < y < b,
∂x 2 ∂y 2
∂u ∂u
(0, y ) = (a, y ) = 0 0 ≤ y ≤ b,
∂x ∂x
u(x, b) = 0, 0 ≤ x ≤ a,
u(x, 0) = f (x), 0 ≤ x ≤ a,
2
[3]
Colombia, Bucaramanga, 2021 30 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Laplaciano en coordenadas polares
En coordenadas rectangulares, el laplaciano tiene la forma
∂2u ∂2u
∆u = + ,
∂x 2 ∂y 2
En coordenadas polares (r , θ), hacemos
x = r cos θ y = r sin θ
El laplaciano en coordenadas polares está dada
∂ 2 u 1 ∂u 1 ∂2u
∆u = + +
∂r 2 r ∂r r 2 ∂θ2
Colombia, Bucaramanga, 2021 31 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Laplaciano en coordenadas polares3
Un disco metálico circular de radio a tiene sus caras aisladas. La orilla del
disco se mantiene a una temperatura dada, que varía de su posición θ. La
temperatura de estado estacionario dentro del disco satisface la ecuación
de Laplace.
∂ 2 u 1 ∂u 1 ∂2u
+ + = 0, 0 ≤ r ≤ a, −π ≤ θ ≤ π,
∂r 2 r ∂r r 2 ∂θ2
u(a, θ) = f (θ), −π ≤ θ ≤ π.
3
[3]
Colombia, Bucaramanga, 2021 32 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Método de solución
Si u(r , θ) = R(r )T (θ) compruebe que se establecen las siguientes
ecuaciones diferenciales ordinarios
r 2 R 00 (r ) + rR 0 (r ) − KR(r ) = 0,
T 00 (θ) + KT (θ) = 0,
donde K es una constante.
Suponga que se imponen las siguientes condiciones de frontera:
u(r , 0) = u(r , π) = 0,
u(r , θ) permanece acotada cuando r → 0+ .
Muestre que una solución no trivial de la forma u(r , θ) = R(r )T (θ) debe
satisfacer las condiciones de frontera
T (0) = T (π) = 0,
R(r ) permanece constante cuando r → 0+ .
Colombia, Bucaramanga, 2021 33 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Laplaciano en coordenadas cilíndricas
En coordenadas rectangulares, el laplaciano tiene la forma
∂2u ∂2u ∂2u
∆u = + + 2,
∂x 2 ∂y 2 ∂z
En coordenadas polares (r , θ), hacemos
x = r cos θ y = r sin θ, z = z,
La ecuación de Laplace se convierte en
∂ 2 u 1 ∂u 1 ∂2u ∂2u
∆u = + + + 2 = 0.
∂r 2 r ∂r r 2 ∂θ2 ∂z
Colombia, Bucaramanga, 2021 34 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Soluciones exactas vs. aproximadas
Los métodos de separación de variables pueden aplicarse tan sólo a un tipo
especial de problemas. Las soluciones que se obtienen mediante este
método, tienen apariencia de exactitud.
Sin embargo, estas exigen procesos de paso al límite que en la mayoría de
los casos no pueden llevarse a cabo, de modo que en cálculo efectivo se
presenta un error.
Colombia, Bucaramanga, 2021 35 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Solución mediante Series de Fourier
Aún cuando puedan encontrarse todos los coeficientes de Fourier de la
solución, sólo en casos muy especiales, es posible hallar la suma de la serie
resultante.
Entonces sólo se puede encontrar una aproximación de la serie de Fourier
de la solución, y por tanto, una solución aproximada. En general, la
solución es una suma parcial de la serie de Fourier.
Además, una función dada en un intervalo [0, L] puede ser representada por
más de una serie de Fourier.
Colombia, Bucaramanga, 2021 36 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Solución mediante Series de Fourier
Una solución en series de Fourier implica la búsqueda de los coeficientes de
Fourier de una función dada (por ejemplo, los valores iniciales). Esto
significa que deben calcularse una infinidad de integrales.
A menos que la función sea muy sencilla, tan sólo un número finito de tales
integraciones pueden efectuarse en un tiempo finito. Incluso, puede ocurrir
que ni siquiera es posible calcular estas integrales exactamente, sino
mediante métodos numéricos.
Colombia, Bucaramanga, 2021 37 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Métodos de solución
1 Métodos directos
2 Separación de variables
3 Transformadas de Integral: Fourier, Laplace
4 Solución de D’ Alembert para la ecuación de onda
5 Cambio de coordenadas
6 Transformación de la variable dependiente
7 Métodos de perturbación
8 Técnica de impulso-respuesta: "divide y reinarás"
Colombia, Bucaramanga, 2021 38 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Elementos finitos
Introducción4
4
[Link]
Colombia, Bucaramanga, 2021 39 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Elementos finitos
Introducción5
5
[Link]
Colombia, Bucaramanga, 2021 40 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Elementos finitos rectangulares6
Figura
6
[Link]/ojs/[Link]/Tecnura/article/view/11559/12387
Colombia, Bucaramanga, 2021 41 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Elementos finitos rectangulares
Un elemento rectangular finito con dimensiones a × b es definida en el
sistema de coordenadas xy . Asuma una función de la forma
φ = A + Bx + Cy + Dxy
y derive la función de "forma" de este elemento rectangular escrita como
φ(x, y ) = N1 (x, y )φ1 + N2 (x, y )φ2 + N3 (x, y )φ3 + N4 (x, y )φ4 .
Colombia, Bucaramanga, 2021 42 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Elementos finitos triangulares
Figura7
7
[Link]
Colombia, Bucaramanga, 2021 43 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Elementos finitos triangulares
Un elemento triangular se define un sistema de coordenadas xy con nodos
1, 2 y 3 localizados en (x1 , y1 ), (x2 , y2 ) y (x3 , y3 ). Derive una función de
forma φ.
The shape function is the function which interpolates the solution between
the discrete values obtained at the mesh nodes.
Asuma una función de interpolación para representar la variación de la
cantidad desconocida como
φ(x, y ) = A + Bx + Cy
Derive la función de "forma" de este elemento rectangular escrita como
φ(x, y ) = N1 (x, y )φ1 + N2 (x, y )φ2 + N3 (x, y )φ3 .
Colombia, Bucaramanga, 2021 44 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Elementos finitos triangulares
Ejemplo
Un elemento triangular tiene nodos en los puntos localizados en
(x1 = 1, y1 = 1), (x2 = 6, y2 = 1), y (x3 = 3, y3 = 4). Una función ha
sido computada para obtener los siguientes valores en los nodos: φ1 = 900,
φ2 = 600 y φ3 = 1200. Haga la interpolación para el elemento triangular y
computar el valor de φ en (x = 3, y = 2).
Colombia, Bucaramanga, 2021 45 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ejercicio8
Suponga que un problema de elemento finito contiene triángulos T1 y T2 .
Calcular φ1 y φ2 .
8
[Link]
Colombia, Bucaramanga, 2021 46 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Formas débiles
Introducción
Esta técnica cambia una ecuación diferencial parcial a una ecuación
integral, donde la incógnita se ubica dentro del integrando. La integral se
resuelve por varias técnicas.
Ejemplo. Formas variacionales
Una función u es solución única de la ecuación diferencial
d 2u
−a + cu = f , 0 < x < L, u(0) = u(L) = 0
dx 2
si, y sólo si, es solución única de
Z L Z L Z L
du dφ
a +c uφ = f φ, φ(0) = φ(L) = 0
0 dx dx 0 0
Colombia, Bucaramanga, 2021 47 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Formas débiles
Problema
Sea p ∈ C 1 , q, f ∈ C [0, 1], y p(x) ≥ δ > 0, q(x) ≥ 0, para 0 ≤ x ≤ 1.
La función y ∈ C 2 [0, 1] es solución única de la ecuación diferencial
d dy
− (p(x) ) + q(x)y = f (x), 0 ≤ x ≤ 1,
dx dx
si y sólo si, y es una función C 2 [0, 1] el cual es el mínimo de la integral
Z 1
I[u] = {p(x)[u 0 (x)]2 + q(x)[u(x)]2 − 2f (x)u(x)}dx
0
Colombia, Bucaramanga, 2021 48 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Método de elementos finitos
Es muy parecido al método de Rayleigh-Ritz. Originalmente se utilizó para
resolver problemas de ingeniería civil, pero sirve para aproximar soluciones
de las ecuaciones diferenciales parciales que surgen en todos los campos de
la ingeniería.
Muchos problemas tienen tipos de condiciones que incluyen derivadas y
fronteras irregulares en sus dominios. Resulta difícil aplicar esta clase de
condiciones de frontera con el método de diferencias finitas.
Dado que cada condición que contenga una derivada debe aproximarse
mediante un cociente de diferencias en los puntos de red, y la forma
irregular de la frontera dificulta la colocación de los puntos de red.
Colombia, Bucaramanga, 2021 49 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación diferencial parcial con condiciones de frontera
Ecuación diferencial parcial con condiciones de frontera
∂ ∂u ∂ ∂u
p(x, y ) + q(x, y ) + r (x, y )u(x, y ) = f (x, y ),
∂x ∂x ∂y ∂y
con (x, y ) ∈ D, donde D es una región plana con frontera s.
Las condiciones de frontera de la forma u(x, y ) = g (x, y ) se imponen en
una parte s1 de la frontera.
Colombia, Bucaramanga, 2021 50 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Ecuación diferencial parcial con condiciones de frontera9
En el resto de la frontera, s2 se requiere que la solución u(x, y ) satisfaga
∂u ∂u
p(x, y ) cos θ1 + q(x, y ) cos θ2 + g1 (x, y )u(x, y ) = g2 (x, y ),
∂x ∂y
donde θ1 y θ2 son los ángulos de dirección de la normal hacia afuera
respecto a la frontera en el punto (x, y ).
9
[4]
Colombia, Bucaramanga, 2021 51 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Método de elementos finitos
El método de elementos finitos incluye las condiciones de frontera como
integrales en una funcional que va reducirse al mínimo, de modo que el
procedimiento de construcción es independiente de las condiciones
particulares de frontera del problema.
Existen tres tipos de condiciones de frontera asociados a una ecuación
diferencial: condición de Dirichlet, condición de Neumann y condición de
Robin.
Colombia, Bucaramanga, 2021 52 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Método de elementos finitos10
El método de elementos finitos aproxima una solución al minimizar el
funcional I en una clase más pequeña de funciones.
El primer paso consiste en dividir la región en una cantidad finita de
elementos de tamaño regular, ya sea rectángulos o triángulos.
10
[4]
Colombia, Bucaramanga, 2021 53 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Método de elementos finitos
Se asume que se subdivide la región D en elementos triangulares definidos
como
(i) (i) (i)
φ(i) (x, y ) = N1 (x, y )φ1 + N2 (x, y )φ2 + N3 (x, y )φ3 ,
(i)
y a sus vértices se les da el nombre de nodos Nk . El método busca una
aproximación de la forma
m
(i)
X
φ(x, y ) = γj φj (x, y )
i=j
donde γ1 , γ2 , . . . , γm son constantes.
Colombia, Bucaramanga, 2021 54 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Método de elementos finitos
Observación
Algunas de estas constantes, se denotan γn+1 , γn+2 , . . . , γm sirven para
asegurarse de que la condición de frontera,
φ(x, y ) = g (x, y ),
se satisfaga en s1 , y las constantes restantes γ1 , γ2 , . . . , γn se emplean para
minimizar el funcional
Xm
I[φ] = I[ γi φ(i) ]
i=1
Colombia, Bucaramanga, 2021 55 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
Sistema de ecuaciones
Para hallar los γi , se establece un sistema de ecuaciones Ac = b, donde
c = (γ1 , . . . , γn )t donde A = (αi j) y b = (β1 , . . . , βn )t se definen por
medio de
Z Z
∂φi ∂φj ∂φi ∂φj
αij = [p(x, y ) +q(x, y ) − r (x, y )φi φj ]dxdy + · · ·
D ∂x ∂x ∂y ∂y
Z
··· + g1 (x, y )φi φj dS,
s2
Z Z Z m
X
βi = − f (x, y )φi dxdy + g2 (x, y )φi dS − αik γk ,
D s2 k=n+1
para toda i = 1, . . . , n y para toda j = 1, 2, . . . , m.
Colombia, Bucaramanga, 2021 56 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
For Further Reading I
Galvis, J., Versieux, H.
Introducción a aproximación numérica de ecuaciones diferenciales
parciales.
28o Coloquio Brasileiro de Matemática
Stanley J. Farlow
Partial Differential Equations for Scientists and Engineers, 1993
Dover Publications, Inc.
Nagle, Saff & Snider.
Ecuaciones diferenciales y problemas con valores en la frontera, 2005
Cuarta Edición, Editorial Pearson, Addison Wesley.
Burden, R. & Douglas-Faires, J.
Análisis numérico, 2002
Séptima Edición, Editorial Thomson Learning.
Colombia, Bucaramanga, 2021 57 /
William González Calderón (UNAB) Modelamiento matemático y simulación avanzada 57
