Ec FisicoMat C01 2021

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 50

METODOS NUMERICOS DE MODELIZACION

MATEMATICA
Introducción a las ecuaciones diferenciales parciales

UNALM

Noviembre 2021

A. Canchoa (Institute) 29 1 / 30
Muchos problemas en ciencias aplicadas, física e ingeniería se modelan
matemáticamente mediante ecuaciones diferenciales parciales.

A. Canchoa (Institute) 29 2 / 30
Muchos problemas en ciencias aplicadas, física e ingeniería se modelan
matemáticamente mediante ecuaciones diferenciales parciales.
Una ecuación diferencial donde aparecen dos o más variables
independientes se llama ecuación en derivadas parciales o una
ecuación diferencial parcial.

A. Canchoa (Institute) 29 2 / 30
Muchos problemas en ciencias aplicadas, física e ingeniería se modelan
matemáticamente mediante ecuaciones diferenciales parciales.
Una ecuación diferencial donde aparecen dos o más variables
independientes se llama ecuación en derivadas parciales o una
ecuación diferencial parcial.
Una ecuación en derivadas parciales de la forma
2 ∂2 Z ∂2 Z
A (x, y ) ∂∂xZ2 + B (x, y ) ∂x ∂Z ∂Z
∂y + C (x, y ) ∂y 2 + f x, y , Z , ∂x , ∂y =0
se clasi…can como una ecuación diferencial parcial elíptica, parabólica
o hiperbólica, como sigue

A. Canchoa (Institute) 29 2 / 30
Muchos problemas en ciencias aplicadas, física e ingeniería se modelan
matemáticamente mediante ecuaciones diferenciales parciales.
Una ecuación diferencial donde aparecen dos o más variables
independientes se llama ecuación en derivadas parciales o una
ecuación diferencial parcial.
Una ecuación en derivadas parciales de la forma
2 ∂2 Z ∂2 Z
A (x, y ) ∂∂xZ2 + B (x, y ) ∂x ∂Z ∂Z
∂y + C (x, y ) ∂y 2 + f x, y , Z , ∂x , ∂y =0
se clasi…can como una ecuación diferencial parcial elíptica, parabólica
o hiperbólica, como sigue
Si B 2 4AC < 0 la ecuación es elíptica.

A. Canchoa (Institute) 29 2 / 30
Muchos problemas en ciencias aplicadas, física e ingeniería se modelan
matemáticamente mediante ecuaciones diferenciales parciales.
Una ecuación diferencial donde aparecen dos o más variables
independientes se llama ecuación en derivadas parciales o una
ecuación diferencial parcial.
Una ecuación en derivadas parciales de la forma
2 ∂2 Z ∂2 Z
A (x, y ) ∂∂xZ2 + B (x, y ) ∂x ∂Z ∂Z
∂y + C (x, y ) ∂y 2 + f x, y , Z , ∂x , ∂y =0
se clasi…can como una ecuación diferencial parcial elíptica, parabólica
o hiperbólica, como sigue
Si B 2 4AC < 0 la ecuación es elíptica.
Si B2 4AC = 0 la ecuación es parabólica.

A. Canchoa (Institute) 29 2 / 30
Muchos problemas en ciencias aplicadas, física e ingeniería se modelan
matemáticamente mediante ecuaciones diferenciales parciales.
Una ecuación diferencial donde aparecen dos o más variables
independientes se llama ecuación en derivadas parciales o una
ecuación diferencial parcial.
Una ecuación en derivadas parciales de la forma
2 ∂2 Z ∂2 Z
A (x, y ) ∂∂xZ2 + B (x, y ) ∂x ∂Z ∂Z
∂y + C (x, y ) ∂y 2 + f x, y , Z , ∂x , ∂y =0
se clasi…can como una ecuación diferencial parcial elíptica, parabólica
o hiperbólica, como sigue
Si B 2 4AC < 0 la ecuación es elíptica.
Si B2 4AC = 0 la ecuación es parabólica.
Si B2 4AC > 0 la ecuación es hiperbólica.

A. Canchoa (Institute) 29 2 / 30
Ejemplo: Ecuación del calor:

∂2 U ∂2 U ∂2 U
+ + = Ut
∂x 2 ∂y 2 ∂z 2

A. Canchoa (Institute) 29 3 / 30
Ejemplo: Ecuación del calor:

∂2 U ∂2 U ∂2 U
+ + = Ut
∂x 2 ∂y 2 ∂z 2

Esta ecuación describe la temperatura U en el tiempo "t" en la


posición (x, y , z ) en un cuerpo tridimensional. El valor de U en el
punto (x, y , z, t ) lo denotaremos por U (x, y , z, t ).

A. Canchoa (Institute) 29 3 / 30
PREDICCION DE TIEMPO Y CLIMA COMO UN
SISTEMA FISICO MATEMATICO DE ECUACIONES

Para Comprender el comportamiento de los sistemas meteorológicos,


en las distintas regiones del mundo, es necesario utilizar un sistema de
ecuaciones físico matemáticos que represente a dichos sistemas.

A. Canchoa (Institute) 29 4 / 30
PREDICCION DE TIEMPO Y CLIMA COMO UN
SISTEMA FISICO MATEMATICO DE ECUACIONES

Para Comprender el comportamiento de los sistemas meteorológicos,


en las distintas regiones del mundo, es necesario utilizar un sistema de
ecuaciones físico matemáticos que represente a dichos sistemas.
Siendo la meteorología dinámica el que trata sobre el estudio de los
movimientos de la atmósfera asociados al tiempo y al clima.

A. Canchoa (Institute) 29 4 / 30
PREDICCION DE TIEMPO Y CLIMA COMO UN
SISTEMA FISICO MATEMATICO DE ECUACIONES

Para Comprender el comportamiento de los sistemas meteorológicos,


en las distintas regiones del mundo, es necesario utilizar un sistema de
ecuaciones físico matemáticos que represente a dichos sistemas.
Siendo la meteorología dinámica el que trata sobre el estudio de los
movimientos de la atmósfera asociados al tiempo y al clima.
Para entender el rol físico que desempeñan los movimientos
atmosféricos en la determinación del clima y el tiempo observados se
hace indispensables desarrollar modelos basados en una simpli…cación
sistemática de las ecuaciones fundamentales (las ecuaciones de
movimiento, ecuación de continuidad y la ecuación de la energía
termodinámica).

A. Canchoa (Institute) 29 4 / 30
1. ECUACIÓN DE MOVIMIENTO:

!
dU ! ! 1 !
= 2Ω U rρ + !
g + Fr (1)
dt ρ

La ecuación (1) constituye la aplicación de la segunda ley de Newton


para un sistema de coordenadas en rotación y se trabaja por unidad
de masa.

A. Canchoa (Institute) 29 5 / 30
1. ECUACIÓN DE MOVIMIENTO:

!
dU ! ! 1 !
= 2Ω U rρ + !
g + Fr (1)
dt ρ

La ecuación (1) constituye la aplicación de la segunda ley de Newton


para un sistema de coordenadas en rotación y se trabaja por unidad
de masa.
Esta ecuación establece que la aceleración correspondiente al
movimiento relativo (término izquierdo de (1)) al sistema en rotación
es igual a la suma de la fuerza de Coriolis ( primer término a la
derecha de (1)), la fuerza de gradiente de presión (segundo término a
la derecha de (1)), la fuerza gravitatoria efectiva (tercer término a la
derecha de (1)) y la fuerza de rozamiento (cuarto término a la
derecha de (1)).

A. Canchoa (Institute) 29 5 / 30
1. ECUACIÓN DE MOVIMIENTO:

!
dU ! ! 1 !
= 2Ω U rρ + !
g + Fr (1)
dt ρ

La ecuación (1) constituye la aplicación de la segunda ley de Newton


para un sistema de coordenadas en rotación y se trabaja por unidad
de masa.
Esta ecuación establece que la aceleración correspondiente al
movimiento relativo (término izquierdo de (1)) al sistema en rotación
es igual a la suma de la fuerza de Coriolis ( primer término a la
derecha de (1)), la fuerza de gradiente de presión (segundo término a
la derecha de (1)), la fuerza gravitatoria efectiva (tercer término a la
derecha de (1)) y la fuerza de rozamiento (cuarto término a la
derecha de (1)).
Esta forma de la ecuación del momento lineal es fundamental en la
mayoría de los aspectos que aborda la meteorología.
A. Canchoa (Institute) 29 5 / 30
El desarrollo de la ecuación (1) da lugar a las ecuaciones (2), (3) y
(4). Las expresiones desarrolladas de las variaciones totales de la
velocidad con respecto al tiempo se presentan en las ecuaciones (5),
(6) y (7).

A. Canchoa (Institute) 29 6 / 30
El desarrollo de la ecuación (1) da lugar a las ecuaciones (2), (3) y
(4). Las expresiones desarrolladas de las variaciones totales de la
velocidad con respecto al tiempo se presentan en las ecuaciones (5),
(6) y (7).

du uv uw 1 ∂p
tan φ + = + 2Ωv senφ 2Ωv cos φ + F(2)
x
dt a a ρ ∂x
dv u 2 uw 1 ∂p
+ tan φ + = 2Ωu senφ + Fy (3)
dt a a ρ ∂y
dw u + v2
2 1 ∂p
= g + 2Ωu cos φ + Fz (4)
dt a ρ ∂z

A. Canchoa (Institute) 29 6 / 30
El desarrollo de la ecuación (1) da lugar a las ecuaciones (2), (3) y
(4). Las expresiones desarrolladas de las variaciones totales de la
velocidad con respecto al tiempo se presentan en las ecuaciones (5),
(6) y (7).

du uv uw 1 ∂p
tan φ + = + 2Ωv senφ 2Ωv cos φ + F(2)
x
dt a a ρ ∂x
dv u 2 uw 1 ∂p
+ tan φ + = 2Ωu senφ + Fy (3)
dt a a ρ ∂y
dw u + v2
2 1 ∂p
= g + 2Ωu cos φ + Fz (4)
dt a ρ ∂z
Las ecuaciones (2) , (3) y (4) son las componentes del momento
lineal en las direcciones este ( eje x), norte (eje y) y vertical (eje z).

A. Canchoa (Institute) 29 6 / 30
El desarrollo de la ecuación (1) da lugar a las ecuaciones (2), (3) y
(4). Las expresiones desarrolladas de las variaciones totales de la
velocidad con respecto al tiempo se presentan en las ecuaciones (5),
(6) y (7).

du uv uw 1 ∂p
tan φ + = + 2Ωv senφ 2Ωv cos φ + F(2)
x
dt a a ρ ∂x
dv u 2 uw 1 ∂p
+ tan φ + = 2Ωu senφ + Fy (3)
dt a a ρ ∂y
dw u + v2
2 1 ∂p
= g + 2Ωu cos φ + Fz (4)
dt a ρ ∂z
Las ecuaciones (2) , (3) y (4) son las componentes del momento
lineal en las direcciones este ( eje x), norte (eje y) y vertical (eje z).
Los términos proporcionales a 1/a en los miembros de la izquierda de
las ecuaciones (2), (3) y (4) se denominan términos de curvatura,
estos términos carecen de importancia para movimientos de escala
sinóptica en latitudes media.
A. Canchoa (Institute) 29 6 / 30
donde

A. Canchoa (Institute) 29 7 / 30
donde
φ: latitud, 1/a: término de curvatura Ω: velocidad angular de la
tierra

A. Canchoa (Institute) 29 7 / 30
donde
φ: latitud, 1/a: término de curvatura Ω: velocidad angular de la
tierra

! ! ! !
U = u i +v j +w k
du ∂u ∂u ∂u ∂u
= +u +v +w (5)
dt ∂t ∂x ∂y ∂z
dv ∂v ∂v ∂v ∂v
= +u +v +w (6)
dt ∂t ∂x ∂y ∂z
dw ∂w ∂w ∂w ∂w
= +u +v +w (7)
dt ∂t ∂x ∂y ∂z

A. Canchoa (Institute) 29 7 / 30
donde
φ: latitud, 1/a: término de curvatura Ω: velocidad angular de la
tierra

! ! ! !
U = u i +v j +w k
du ∂u ∂u ∂u ∂u
= +u +v +w (5)
dt ∂t ∂x ∂y ∂z
dv ∂v ∂v ∂v ∂v
= +u +v +w (6)
dt ∂t ∂x ∂y ∂z
dw ∂w ∂w ∂w ∂w
= +u +v +w (7)
dt ∂t ∂x ∂y ∂z
Las ecuaciones (5), (6) y (7) que viene a ser el desarrollo de las
derivadas totales:

A. Canchoa (Institute) 29 7 / 30
donde
φ: latitud, 1/a: término de curvatura Ω: velocidad angular de la
tierra

! ! ! !
U = u i +v j +w k
du ∂u ∂u ∂u ∂u
= +u +v +w (5)
dt ∂t ∂x ∂y ∂z
dv ∂v ∂v ∂v ∂v
= +u +v +w (6)
dt ∂t ∂x ∂y ∂z
dw ∂w ∂w ∂w ∂w
= +u +v +w (7)
dt ∂t ∂x ∂y ∂z
Las ecuaciones (5), (6) y (7) que viene a ser el desarrollo de las
derivadas totales:
en sus componentes local: ( ∂u ∂v ∂w
∂t , ∂t , ∂t ) y

A. Canchoa (Institute) 29 7 / 30
donde
φ: latitud, 1/a: término de curvatura Ω: velocidad angular de la
tierra

! ! ! !
U = u i +v j +w k
du ∂u ∂u ∂u ∂u
= +u +v +w (5)
dt ∂t ∂x ∂y ∂z
dv ∂v ∂v ∂v ∂v
= +u +v +w (6)
dt ∂t ∂x ∂y ∂z
dw ∂w ∂w ∂w ∂w
= +u +v +w (7)
dt ∂t ∂x ∂y ∂z
Las ecuaciones (5), (6) y (7) que viene a ser el desarrollo de las
derivadas totales:
en sus componentes local: ( ∂u ∂v ∂w
∂t , ∂t , ∂t ) y
advectiva: (u ∂u ∂u ∂u ∂v ∂v ∂v ∂w ∂w ∂w
∂x + v ∂y + w ∂z ,u ∂x + v ∂y + w ∂z ,u ∂x + v ∂y + w ∂z ).

A. Canchoa (Institute) 29 7 / 30
donde
φ: latitud, 1/a: término de curvatura Ω: velocidad angular de la
tierra

! ! ! !
U = u i +v j +w k
du ∂u ∂u ∂u ∂u
= +u +v +w (5)
dt ∂t ∂x ∂y ∂z
dv ∂v ∂v ∂v ∂v
= +u +v +w (6)
dt ∂t ∂x ∂y ∂z
dw ∂w ∂w ∂w ∂w
= +u +v +w (7)
dt ∂t ∂x ∂y ∂z
Las ecuaciones (5), (6) y (7) que viene a ser el desarrollo de las
derivadas totales:
en sus componentes local: ( ∂u ∂v ∂w
∂t , ∂t , ∂t ) y
advectiva: (u ∂u ∂u ∂u ∂v ∂v ∂v ∂w ∂w ∂w
∂x + v ∂y + w ∂z ,u ∂x + v ∂y + w ∂z ,u ∂x + v ∂y + w ∂z ).
En general, los términos advectivos de la aceleración son en magnitud
comparable a la aceleración local.
A. Canchoa (Institute) 29 7 / 30
1.1 Análisis de escala de las ecuaciones horizontales del
momento lineal

En el análisis de escala de las ecuaciones del movimiento, con el


objeto de determinar si alguno de los términos de las ecuaciones es
despreciable, en lo que concierne a los fenómenos meteorológicos.

A. Canchoa (Institute) 29 8 / 30
1.1 Análisis de escala de las ecuaciones horizontales del
momento lineal

En el análisis de escala de las ecuaciones del movimiento, con el


objeto de determinar si alguno de los términos de las ecuaciones es
despreciable, en lo que concierne a los fenómenos meteorológicos.
Esta eliminación de términos en base a estas consideraciones no
solamente presenta la ventaja del tratamiento matemático, sino que
la eliminación de los términos de valores pequeños poseen la
importantísima propiedad en determinados casos de eliminar
completamente o …ltrar un tipo de movimiento indeseado.

A. Canchoa (Institute) 29 8 / 30
Con el objeto de simpli…car las ecuaciones de movimiento (2) y (3) de
escala sinóptica, se puede de…nir las siguientes escalas características de
las variables de campo en base a valores observados para sistemas
sinópticos en latitudes medias.
U 10ms 1 : escala de la velocidad horizontal
W 1cms 1 : escala de la velocidad vertical
L 6
10 m : escala de longitud
D 104 m : escala de altura
∆P 3 m2 s 2
ρ 10 : escala de la variación horizontal de la presión.
L 5
U 10 s : escala de tiempo ( es el tiempo requerido
para recorrer una distancia L a la velocidad U).
Es conveniente considerar las perturbaciones centradas en la latitud φ0 o
igual a 45o e introducir la expresión

4 1
f/0 = 2Ωv senφ0 = f/0 = 2Ωv cos φ0 10 s .

A. Canchoa (Institute) 29 9 / 30
En la tabla siguiente se indica la magnitud característica de cada término
de las ecuaciones (2) y (3) en base a las consideraciones de escala
mencionadas anteriormente.
No se incluyen los términos de rozamiento, porque para la escala de
tiempo sinóptica la disipación por fricción es de importancia secundaria
por encima del nivel de 100 metros.
A B C D E F
Ec. de la comp. du uw uv 1 ∂p
dt 2Ωv senφ + 2Ωv cos φ + a a tan φ =
x del mov. lineal ρ ∂x

Ec. de la comp. dv uw u2 1 ∂p
dt + 2Ωu senφ + a + a tan φ =
y del mov. lineal ρ ∂y

Escalas de térmi- U2 UW U2 ∆P
L f0 U f0 W a a ρL
nos individuales
Magnitudes de
términos 10 4 10 3 10 6 10 8 10 5 10 3

individuales
Tabla. Análisis de escala de las ecuaciones horizontales del momento lineal.

A. Canchoa (Institute) 29 10 / 30
1.2 APROXIMACION GEOSTROFICA Y EL VIENTO
GEOSTROFICO

Para sistemas atmosféricos de latitudes medias y escala sinóptica, teniendo


en cuenta la tabla anterior el término B es la fuerza de Coriolis y el
término F es la fuerza del gradiente de presión que son aproximadamente
del mismo orden de magnitud ( 10 3 ).
Considerando únicamente estos dos términos en las ecuaciones (2) y (3) y
du dv
dt = 0, dt = 0 se obtiene la relación geostró…ca :

1 ∂p
2Ωv senφ ' (8)
ρ ∂x
1 ∂p
2Ωu senφ ' (9)
ρ ∂y

A. Canchoa (Institute) 29 11 / 30
Reemplazando el parámetro de Coriolis f = 2Ωu senφ en las ecuaciones
(8) y (9), se obtienen
1 ∂p
f v' (10)
ρ ∂x
1 ∂p
f u' (11)
ρ ∂y

A. Canchoa (Institute) 29 12 / 30
El campo de la velocidad horizontal en una aproximación geostró…ca se
de…ne como:
! ! !
Vg = ug i + vg j ,
denominado viento geostró…co, que también satisface las ecuaciones (10)
y (11). Por consiguiente en forma vectorial
! ! 1
Vg = k rp. (12)
ρf
De la ecuación (12), se deduce que la distribución de la presión
atmosférica en cualquier instante determina el viento geostró…co.
Debe tenerse en cuenta que la ecuación (12) de…ne el viento geostró…co;
pero sólo para movimientos atmosféricos a gran escala, debe utilizarse el
viento geostró…co como aproximación al campo del viento horizontal real.

A. Canchoa (Institute) 29 13 / 30
2. ECUACIÓN DE ENERGIA

Utilizando el principio de la conservación de energía y aplicado a un


elemento de ‡uido en movimiento.
El primer principio de la termodinámica se postula por lo general
considerando un sistema en equilibrio termodinámico, es decir un sistema
que se encuentra inicialmente en reposo y a continuación tras intercambiar
calor con sus alrededores y realizar trabajo sobre ellos vuelve de nuevo al
reposo.
Para un sistema de este tipo el primer principio establece que la variación
de la energía interna del sistema es igual a la diferencia entre el calor
aportado al sistema y el trabajo realizado por el mismo.

∆U = ∆Q ∆W

A. Canchoa (Institute) 29 14 / 30
Ecuación de la energía termodinámica :

q
q:‡ujo de calor que ingresa a la parcela
q (J/s )
En un intervalo de tiempo 4t, el balance r
de calor es:q 4t = 4e + 4w U
(velocida

q 4t : Flujo de calor por unidad de tiempo

4e : cambio de energía interna


PARCELA MOVIL
:4w :trabajo debido a la expansión

A. Canchoa (Institute) 29 15 / 30
Fórmula para 4w : Suponga un gas encerrado en un cilindro con un
pistón movible.

F
A Fuerza sobre el pistón : F = pA
Trabajo en la expansión :
F 4x = pA4x
F 4x = p 4 α
4α :cambio de volumen
δx

Además: 4e = Cv 4T
Balance de calor : q 4t = Cv 4T + p 4α
Cuando 4t ! 0:
dT dα
q = Cv +p (13)
dt dt
A. Canchoa (Institute) 29 16 / 30
dT dα
q = Cv +p (13)
dt dt
La ecuación (13) es la forma usual de la ecuación de energía. De esta
forma la primera ley de la termodinámica es aplicable a un ‡uido en
movimiento.
El segundo término a la derecha de (13) representa el trabajo efectuado
por unidad de tiempo (y por unidad de masa) por el sistema ‡uido y que
representa un término de conversión de energía térmica a mecánica. Este
es el proceso de conversión que permite a la energía calorí…ca del sol
gobernar los movimientos de la atmósfera.
La información que se utiliza en la ecuación anteriormente descrita se
requieren:
-Datos de temperatura del aire
-Presión atmosférica
-Densidad del aire
-Cv (calor especi…co para el aire seco a volumen constante:
717 J kg 1 K 1 )
A. Canchoa (Institute) 29 17 / 30
3. ECUACION DE CONTINUIDAD

La ecuación de continuidad trata de expresar que si un conjunto de


partículas de un ‡uido ocupan en un instante determinado, un volumen
limitado por una super…cie ‡uida, esa super…cie se comporta como
impenetrable, aunque deformable; de manera que la masa de ‡uido que
hay en su interior no varía, aunque dicha super…cie sufra deformaciones
que hagan variar el volumen que delimita.

A. Canchoa (Institute) 29 18 / 30
3.1 Interpretación Física

La ecuación de continuidad se puede expresar de varias formas, por


ejemplo una forma más común es la siguiente:

1 ∂ρ !
+r U = 0 (14)
ρ ∂t

De la expresión (14) podemos mencionar que el primer termino nos indica


la variación relativa de la densidad por unidad de tiempo que acompaña al
movimiento de una burbuja de aire es compensado con la divergencia de la
velocidad.

A. Canchoa (Institute) 29 19 / 30
La ecuación de continuidad se puede expresar en coordenadas horizontales,
isobáricas, otros.
En el caso de coordenadas horizontales la ecuación de continuidad se
puede expresar como:

1 ∂α !
=r U (15)
α ∂t
donde :
α = 1/ρ : es el volumen especí…co
! ! !
U = u i +v j
El primer término de la ecuación (15) indica la variación local del volumen
especí…co, mientras que el segundo término nos indica la divergencia de la
velocidad. También podemos mencionar que la variación temporal del
volumen especí…co, esta compensado por la divergencia de la velocidad.
La ecuación (15) es válido para ‡uidos compresibles, como es el caso del
aire.

A. Canchoa (Institute) 29 20 / 30
Es muy importe mencionar que la divergencia de la velocidad mide la
dilatación relativa del aire por unidad de tiempo.
En coordenadas isobáricas la ecuación de continuidad se puede expresar de
la siguiente forma:

∂p ! ∂p
= rP U , p= =w (16)
∂t ∂t
En este caso la divergencia isobárica es el origen único de las variaciones
de la velocidad de ascenso y descenso a lo largo del eje z en términos de la
presión.

A. Canchoa (Institute) 29 21 / 30
4.REQUERIMIENTO DE DATOS, ECUACIONES,
RESOLUCIÓN Y TIPOS DE MODELOS NUMERICOS

Para la predicción numérica del tiempo y clima se utilizan las ecuaciones


físicas que gobiernan el comportamiento de la atmósfera las que fueron
descritas anteriormente (ecuación de movimiento, de continuidad y energía
termodinámica) y se formulan modelos matemáticos al cual ingresamos
datos sobre el estado pasado y actual del tiempo atmosférico y se
solucionan utilizando métodos físicos matemáticos para predecir un estado
futuro de la atmósfera.
Esto es sin duda un proceso muy complejo debido a la necesidad de
resolver ecuaciones diferenciales no lineales, no es posible aplicar los
métodos analíticos, en vez de ello es necesario el uso de aproximaciones
numéricas con un fuerte trabajo computacional.

A. Canchoa (Institute) 29 22 / 30
Además, la naturaleza no lineal de estas ecuaciones implica que pequeñas
diferencias en los datos de entrada al modelo producirá grandes diferencias
en los resultados, esto es lo que se conoce como un sistema caótico, una
fuerte sensibilidad a las condiciones iniciales, de allí la importancia de la
información que alimentará al modelo.
Esta información es recopilada de las estaciones meteorológicas alrededor
del mundo, de boyas, estaciones meteorológicas en barcos o aviones,
radiosondas, radares e información satelital, hacia un Sistema de
Telecomunicación Global (GTS), de manera que este disponible para
inicializar cualquier corrida de modelo con un previo proceso de
transformación a la grilla del modelo de pronóstico (esquema de
asimilación), distribuyendo los datos uniformemente en el tiempo y el
espacio.

A. Canchoa (Institute) 29 23 / 30
Otro aspecto importante en la resolución de modelos numéricos del tiempo
es la ocurrencia de fenómenos de pequeña escala (temporal o espacial) que
tienen in‡uencia en el comportamiento de la atmósfera no podrían ser
cubiertas por las grilla del modelo, procesos como la formación de nubes,
la precipitación, los procesos de radicación, la liberación de calor latente,
los procesos de capa limite, etc. son descritos mediante esquemas de
parametrización. Este es un aspecto de suma importancia en los
pronósticos numéricos de corto plazo, siendo poco considerados para
pronósticos de mediano y largo plazo.
Las ecuaciones básicas que rigen el movimiento atmosférico esta
gobernado por las leyes física que son descritas por la ecuaciones de la
hidrodinámica y termodinámica.

A. Canchoa (Institute) 29 24 / 30
La mayoría de los modelos numéricos atmosféricos consideran estas
ecuaciones para la simulación de la dinámica atmósferica. Las ecuaciones
fundamentales que gobiernan el movimiento de la atmósfera son:
- La ecuación del movimiento horizontal (ver ecuaciones (2) y (3))
- La ecuación hidrostática

∂p
= ρg
∂Z
- Ecuación de continuidad (ver ecuación (14))
- Ecuación de la energía termodinámica (ver ecuación (13))
- Ecuación de estado : pV = nRT donde R = Cp CV

A. Canchoa (Institute) 29 25 / 30
Estas ecuaciones no ingresan al modelo de esta manera, primero deben ser
transformadas a relaciones de series matemáticas de manera que lo
entienda el computador. Existen dos métodos muy usados al respecto :
- diferencias …nitas (series truncadas de Taylor)
- método espectral (mediante el uso de senos y cosenos de la
trasformada de Fourier y Legendre)
- La precisión de un modelo está directamente relacionado a la
resolución horizontal y vertical, mientras más pequeña sea la resolución el
modelo simulará mejor la atmósfera real. Pero a su vez el incremento de la
resolución trae consigo un mayor trabajo computacional, hasta inclusive
abortar los procesos si no se cuenta con una computadora adecuada para
el procesamiento, porque la reducción de la resolución horizontal provoca
inestabilidad en el cálculo de diferencias …nitas.

A. Canchoa (Institute) 29 26 / 30
Esto fue descubierto en 1928 por Courant, Friedrichs y Levi denominado
comúnmente CFL o criterio de estabilidad para las diferencias …nitas el
cual debe cumplir la condición c 4 t/ 4 x 1, esto signi…ca que para 4x
pequeños se requiere de 4t pequeños a …n de asegurar la estabilidad
computacional. Se consideran 4t pequeño para movimientos de ondas
rápidos, como las ondas de gravedad debido a que c(velocidad de la onda)
pequeña. El incremento de la resolución sin perder estabilidad se podría
considerar cuando resolvemos las ecuaciones para áreas pequeñas como
ciudades o países y para periodos cortos de tiempo. La manera más
utilizada y complicada es usando grillas …nas anidadas sobre pequeñas
áreas de interés particular y éstas dentro de las grillas iniciales.

A. Canchoa (Institute) 29 27 / 30
4.1 TIPOS DE MODELOS
En la actualidad en el mundo existen más de 20 modelos de pronóstico
numérico trabajando operativamente. Los modelos presentan una
secuencia similar, constan de tres pasos principales: El pre-proceso, el
proceso propiamente dicho y el post-proceso.
La diferencia entre los modelos radica en los métodos que emplea para
resolver las ecuaciones que gobiernan el movimiento del aire, la resolución,
el tiempo de pronóstico para corto plazo para uno, dos o siete días el
dominio o área de pronóstico regional, local o global. Con respecto a los
modelos de pronóstico numérico a corto plazo más usados en EEUU son
los siguientes:
El ETA, con 29m y 22 km de resolución, usa el sistema de
coordenada “eta”.
El AVN, Modelos de Aviación, con 100km de resolución, es un
modelo antiguo vigente hasta hoy, brinda pronósticos para más de 72
horas.
El NGM, el modelo de grilla anidada con 80 km de resolución,
similar al ETA, predice de 2 días en adelante.
A. Canchoa (Institute) 29 28 / 30
4.2 INFORMACIÓN REQUERIDA
Existen otros modelos para predicciones a mediano y largo plazo,
mencionaremos algunos de ellos:
El MRF (Medium Range Forecast) 150Km de resolución.
El Modelo Océano Global, produce pronósticos estacionales de
temperatura de agua de mar y corrientes oceánicas.
El MM5 (Mesoscale Model versión 5) es un modelo de investigación.
El RAMS (Regional Atmospheric Model System), de investigación.
Los modelos numéricos deben incluir la mayor cantidad posible de la
información de:
- Super…cie, información de las variables meteorológicas de las
estaciones automáticas y convencionales de super…cie (velocidad y
dirección de vientos, presión atmosférica, temperatura del aire,
temperatura del rocío y otros).
- Altura, información de las variables meteorológicas de las estaciones
automáticas y convencionales de altura y información de aeronaves-AIREP
(velocidad y dirección de vientos, altura geopotencial, temperatura del
aire, temperatura del rocío y otros).
A. Canchoa (Institute) 29 29 / 30
BIBLIOGRAFÍA :
- Sensores remotos, información meteorológica de la troposfera
(temperatura, temperatura de agua de mar, humedad, viento).
- Radar, información de fenómenos meteorológicos.
- Training Course on Weather Forecasting for Operational
Meteorologists. 10-30 April 2000. Meteorogical Training Division, Korea
Meteorogical Administration
- Proceedings International Seminar of Tropical Meteorology.
Campinas (SP)- Brazil September 25 – October 10 – 1969.
- Holton, James R. 1979 An introduction to Dynamic Meteorology.
ACADEMIC PRESS, INC. Orlando.
- Medina, M. 1984 Teoría de la predicción meteorológica. Instituto
Nacional de Meteorología – España.
- Gordon, A. H. 1965 Elementos de Meteorología Dinámica. UTEHA-
México.
- Pronóstico Numérico. Manual de PCGRIDS y GRADS.
Departamento de Física y Meteorología - UNALM 2001.
- Haltiner, G. J. y R.T. Williams, 1980: Numerical Prediction and
A. Canchoa (Institute) 29 30 / 30

También podría gustarte