Fórmula de cálculo de la temperatura del termómetro húmedo

Iniciado por Gerard Taulé, Miércoles 15 Octubre 2008 23:00:35 PM

Tema anterior - Siguiente tema

ManoloZ

*****
Cb Incus
Mensajes: 3,765
  1. manolokosh@yahoo.es
Ubicación: Zaragoza
En línea
#12
Por justificar un poco las formulas, que me parece mal sacarlas de la manga
la idea es que la entalpía perdida en el enfriamiento es la que se gasta en evaporar, o sea
(cp+m*c'p)(t-t')=L'(M'-m)
donde Cp escalor especifico a P cte aire seco
C'p idem vapor agua
L': calor latente a t'
M': razón de mezcla saturada a t'
m: razón mezcla a t
si haceis las cuentas sale algo parecido (bueno, sale 1488 en vez de 1510 pero es una formula aproximada)

edit:
m=0.622 e/p
M=0.622 E/p
y os recuerdo que es aproximado
Zaragoza: http://www.flickr.com/photos/manolokosh/
Uno que madrugó un mensaje se encontró; pero más madrugó el que lo posteó.
Lexibook SM1660 "La Cutrestación" (TM)-"Otra-cutrestación" (TM) ALDI (Globaltronics GT-WS-04)

FRENTEFRIO

*****
Cb Incus
Mensajes: 4,980
Ojáiz-Peñacastillo (Cantabria)
  1. ojaizmet@hotmail.com
En línea
#13
Saludos.

Podéis ver muchas fórmulas en el enlace siguiente:

http://www.aprweather.com/pages/calc.htm

Luego sólo es pasar las que interesen a una hoja de cálculo y sustituir las siglas de las fórmulas por las celdas. Así tengo yo hechas varias de ellas.
Espero que sea de ayuda.

PD: La fórmula que busca Gerard esta al final de todo.

Estación Meteorológica de Ojáiz-Peñacastillo (Cantabria) Página web: OJAIZMETEO
DETECTOR DE RAYOS
                                                   Vigilando el tiempo

Gerard Taulé

*****
Cb Incus
Mensajes: 4,418
Ubicación: Salt
En línea
Cita de: FRENTEFRIO en Sábado 18 Octubre 2008 13:12:50 PM
Saludos.

Podéis ver muchas fórmulas en el enlace siguiente:

http://www.aprweather.com/pages/calc.htm

Luego sólo es pasar las que interesen a una hoja de cálculo y sustituir las siglas de las fórmulas por las celdas. Así tengo yo hechas varias de ellas.
Espero que sea de ayuda.

PD: La fórmula que busca Gerard esta al final de todo.



Muchas gracias, Frentefrío, esta fórmula es relativamente sencilla comparando con las fórmulas que he encontrado por ejemplo en el Calendario Meteorológico 1999.
Salt, 3 km al oeste de Girona, 81 m, isla de calor urbana

JanusSX

*
Sol
Mensajes: 12
En línea
#15
Hola:
La única manera para hallar con buena precisión la temperatura de bulbo húmedo es usando la ecuación de Claysius Clayperon, derivarla  y resolverla por métodos iterativos como el de Rapson Newton, para ello requieres los siguientes datos previos:
Datos q tienes:
-HR (%)
-ta (temp del aire °C)
Necesitamos "e" y "es":
Con t--->es
con td--->ed   (td: temp de rocío en °C)
donde:
es: es la presión de vapor de saturación
ed: presión de vapor actual
HR=e/es*100 (forma aproximada ok)
es=6.11*10^(at/(b+t))     ... Teetens(1930) con t en °C, y "es" en mb.
donde : a=7.5 y b=237.3 para t>=0
y:          a=9.5 y b=265.5 para t<=0
Entonces:
e=100*HR/es ( "e" en mb)
Ahora usaremos el método usado por la NOAA-AWIPS:
Calculadas las presiones "e" y "es" hacemos:
s=(es-ed)/(T-Td), luego:
1ra. aproximación de Tw=(T.f.p+Td.s)/(f.p+s)
Donde:
Tw, T, Td: temp de bulbo húmedo, temp del aire y temp rocío en °K
p=presión en mb
f=0.0006355 = C p /(L*epsilon) (1/K)
Dado nuestro primer supuesto y mediante método del balance de energía, la diferencia "de" deberá ser evaluada:
de=f.p(T-Tw)-(ew-ed)  ....Irribarbe & Godson 1981
Donde:
ew: presión de saturación de bulbo húmedo en mb.
ew=EXP(C15-C1.Tw-C2/Tw)
EXP: base log neperiano
Tw: temp de bulo húmedo en °K
C15=26.66082
C1=0.0091379024
C2=6106.396
Cuando la diferencia "de" es menor de 1/10,000, detenemos la iteración, quedándonos con dicho valor.
Luego, hallamos la derivada "der" de la diferencia "de" con respecto la temp de bulbo húmedo, la cual es:
der=ew.(C1-C2/(Tw^2)).f.p
Siendo finalmente nuestra próxima aproximación:
Tw=Tw-de/der, este nuevo valor es probado hasta por lo menos 4 iteraciones.

Resumiendo, con tu primer "Tw0", calculas tu primer "de0" como no es menor de 1/10,000 prosigues calculando "der0" con lo que obtienes tu próximo:
Tw1=Tw0-de0/der0
luego de manera repetitiva hallas el nuevo "de1" y al no ser menor de 1/10,000 prosigues calculando "der1" con lo obtienes tu próximo:
Tw2=Tw1-de1/der1
y así sucesivamente hasta que la dieferencia "de" sea menor de 1/10,000
Por experiencia con 4 a 5 iteraciones es suficiente.
Estos cálculos son sencillos de realizar elaborando un programa u hoja de cálculo.
Si deseas envíame un emal a: [email protected] te enviaré un programa u hoja de cálculo q tengo hecho al respecto.
Saludos...