[DZvents] alerte givre et calcul humidité absolue

Réécriture du script lua en DZvents avec gestion des systèmes de notification, spécification des plages horaires de notification et calcul de l’humidité absolue, (le fameux « poids d’eau ») afin de corriger le seuil d’alerte en fonction de celle-ci.

La définition du point de givrage d’après Wikipédia

Le point de givrage, connu également comme la température du point givrage ou le point de gelée, est une donnée météorologique calculée à partir de l’humidité, la pression et la température. Le point de givrage de l’air est la température à laquelle, tout en gardant inchangées les conditions barométriques courantes, l’air devient saturé de vapeur d’eau par rapport à la glace. Le point de givrage est donc l’équivalent du point de rosée pour la condensation de la vapeur d’eau directement en cristaux de glace et non en micro-gouttelettes. C’est le phénomène de déposition, qui survient lorsque le point de givre est atteint, qui crée la gelée blanche.

maintenant la formule :
point_givre


Il nous faut donc dans un premier temps calculer le point de rosée nécessaire au calcul du point de givrage.

la formule du point de rosée.

function dewPoint (T, RH)
  local b,c = 17.67, 243.5
  RH = math.max (RH or 0, 1e-3)
  local gamma = math.log (RH/100) + b * T / (c + T) 
  return c * gamma / (b - gamma)
end

et maintenant la formule point de givrage

function freezing_point(dp, t) 
T = t + 273.15
Td = dp + 273.15
return (Td + (2671.02 /((2954.61/T) + 2.193665 * math.log(T) - 13.3448))-T)-273.15
end

la formule pour calculer l’humidité absolue

Humidité absolue (grammes / m3 ) =  (6,122 * e^[(17,67 * T) / (T + 243,5)] * rh * 2,1674))/(273,15 + T)

Cette formule est précise à 0,1% près, dans la gamme de température de -30 ° C à + 35 ° C

Dans la formule ci-dessus, la température (T) est exprimée en degrés Celsius, l’humidité relative (hr) est exprimée en %, et e est la base des logarithmes naturels 2.71828 [élevée à la puissance du contenu des crochets].

le script complet pour Domoticz

--[[
~/home/pi/domoticz/scripts/dzVents/scripts/alerteGivre.lua
auteur : papoo
MAJ : 05/05/2019
création : 06/05/2016
Principe :
Calculer via les informations température et hygrométrie d'une sonde extérieure
 le point de rosée, le point de givre et l'humidité absolue
pour générer une alerte givre cohérente.
http://pon.fr/dzvents-alerte-givre-et-calcul-humidite-absolue
http://easydomoticz.com/forum/viewtopic.php?f=21&t=1085&start=10#p17545
https://github.com/papo-o/domoticz_scripts/blob/master/dzVents/scripts/alerteGivre.lua
]]--

--------------------------------------------
------------ Variables à éditer ------------
--------------------------------------------
local device_temp_ext       = 'Temperature exterieure'     -- nom (entre ' ') ou idx  de la sonde de température/humidité extérieure
local device_dew_point      = 'Point de rosée'             -- nom (entre ' ') ou idx de l'éventuel dummy température point de rosée si vous souhaitez le suivre sinon nil
local device_freeze_point   = 'Point de givrage'           -- nom (entre ' ') ou idx du dummy température point de givre si vous souhaitez le suivre sinon nil
local device_hum_abs        = 'Humidité absolue'           -- nom (entre ' ') ou idx de l'éventuel dummy humidité absolue si vous souhaitez le suivre sinon nil
local device_freeze_alert   = 'Risque de givre'            -- nom (entre ' ') ou idx du dummy alert point de givre
local thresholdHumidity     = 2.8                          -- seuil humidité absolue en dessous duquel il est peu probable qu'il givre
local subject               = "/!\\ Attention /!\\"        -- sujet des notifications
local scriptName = 'Alerte Givre'
local scriptVersion = '2.2'
local message = '.'


return {
    active =  true,  -- false,
    on = {
        timer   = { 'at 06:25', 'at 07:01', 'at 07:15'}, -- heures de déclenchement des notifications en cas de givre
        devices = { device_temp_ext }     -- nom de la sonde de température/humidité extérieure'
        },

    logging =   {   -- level    =   domoticz.LOG_DEBUG,
                    -- level    =   domoticz.LOG_INFO,             -- Seulement un niveau peut être actif; commenter les autres
                    -- level    =   domoticz.LOG_ERROR,            -- Only one level can be active; comment others
                    -- level    =   domoticz.LOG_MODULE_EXEC_INFO,

                    marker = scriptName..' v'..scriptVersion },

    execute = function(domoticz, item)

        local SubSystem = "TELEGRAM, HTTP"
                            --[[ Systèmes de notification disponibles :
                                GOOGLE_CLOUD_MESSAGING, HTTP, KODI, LOGITECH_MEDIASERVER, NMA, PROWL, PUSHALOT, PUSHBULLET,PUSHOVER, PUSHSAFER, TELEGRAM
                                Pour une notification sur plusieurs systèmes, séparez les systèmes par une virgule.
                                Exemple :"PUSHBULLET, HTTP"
                            --]]

--------------------------------------------
----------- Fin variables à éditer ---------
--------------------------------------------

        function split(s, delimiter)
        result = {};
        for match in (s..delimiter):gmatch("(.-)"..delimiter) do
            table.insert(result, match);
        end

        return result;
        end
        function notificationTable(str)
        --NSS_GOOGLE_CLOUD_MESSAGING, NSS_HTTP, NSS_KODI, NSS_LOGITECH_MEDIASERVER, NSS_NMA,NSS_PROWL, NSS_PUSHALOT, NSS_PUSHBULLET, NSS_PUSHOVER, NSS_PUSHSAFER, NSS_TELEGRAM

            if (str) then
            str = string.gsub (str,"GCM", domoticz.NSS_GOOGLE_CLOUD_MESSAGING)
            str = string.gsub (str,"GOOGLE_CLOUD_MESSAGING", domoticz.NSS_GOOGLE_CLOUD_MESSAGING)
            str = string.gsub (str,"HTTP", domoticz.NSS_HTTP)
            str = string.gsub (str,"LMS", domoticz.NSS_LOGITECH_MEDIASERVER)
            str = string.gsub (str,"LOGITECH_MEDIASERVER", domoticz.NSS_LOGITECH_MEDIASERVER)
            str = string.gsub (str,"NMA", domoticz.NSS_NMA)
            str = string.gsub (str,"PROWL", domoticz.NSS_PROWL)
            str = string.gsub (str,"PUSHALOT", domoticz.NSS_PUSHALOT)
            str = string.gsub (str,"PUSHOVER", domoticz.NSS_PUSHOVER)
            str = string.gsub (str,"PUSHSAFER", domoticz.NSS_PUSHSAFER)
            str = string.gsub (str,"PUSHBULLET", domoticz.NSS_PUSHBULLET)
            str = string.gsub (str,"TELEGRAM", domoticz.NSS_TELEGRAM)
            end
        return (split(str,','))
        end

        function freezing_point(dp, t)
            if not dp or not t or dp > t then
                return nil, " La température du point de rosée est supérieure à la température. Puisque la température du point de rosée ne peut être supérieure à la température de l'air , l\'humidité relative a été fixée à nil."
            end
            T = t + 273.15
            Td = dp + 273.15
            return domoticz.utils.round((Td + (2671.02 /((2954.61/T) + 2.193665 * math.log(T) - 13.3448))-T)-273.15, 2)
        end

        function hum_abs(t,hr)
            -- https://carnotcycle.wordpress.com/2012/08/04/how-to-convert-relative-humidity-to-absolute-humidity/
            -- Formule pour calculer l'humidité absolue
            -- Dans la formule ci-dessous, la température (T) est exprimée en degrés Celsius, l'humidité relative (hr) est exprimée en %, et e est la base des logarithmes naturels 2.71828 [élevée à la puissance du contenu des crochets]:
            -- Humidité absolue (grammes / m3 ) =  (6,122 * e^[(17,67 * T) / (T + 243,5)] * rh * 2,1674))/(273,15 + T)
            -- Cette formule est précise à 0,1% près, dans la gamme de température de -30 ° C à + 35 ° C
            ha = domoticz.utils.round((6.112 * math.exp((17.67 * t)/(t+243.5)) * hr * 2.1674)/ (273.15 + t),1)
            return ha
        end

        temperature = domoticz.devices(device_temp_ext).temperature
        humidity = domoticz.devices(device_temp_ext).humidity
        dewPoint = domoticz.utils.round(domoticz.devices(device_temp_ext).dewPoint, 2)
        freezingPoint = freezing_point(dewPoint, temperature)
        absHumidity = hum_abs(temperature,humidity)
        domoticz.log("--- --- --- Température Ext :     ".. domoticz.utils.round(temperature, 2) .." °C", domoticz.LOG_DEBUG)
        domoticz.log("--- --- --- Humidité :            ".. humidity .." %", domoticz.LOG_DEBUG)
        domoticz.log("--- --- --- point de rosée :      ".. dewPoint .." °C", domoticz.LOG_DEBUG)
        domoticz.log("--- --- --- point de givre :      ".. freezingPoint .." °C", domoticz.LOG_DEBUG)
        domoticz.log("--- --- --- Humidité Absolue :    ".. absHumidity .." g/m³", domoticz.LOG_DEBUG)

        if device_dew_point    ~= nil  then domoticz.devices(device_dew_point).updateTemperature(dewPoint) end -- Mise à jour point de rosée
        if device_freeze_point ~= nil  then domoticz.devices(device_freeze_point).updateTemperature(freezingPoint) end -- Mise à jour point de givrage
        if device_hum_abs      ~= nil  then domoticz.devices(device_hum_abs).updateCustomSensor(absHumidity) end -- Mise à jour humidité absolue

        if(tonumber(temperature) <= 1 and tonumber(freezingPoint) <= 0) and absHumidity > thresholdHumidity then
            domoticz.log("--- --- --- Givre --- --- ---", domoticz.LOG_DEBUG)
            message = "Présence de givre"
            if tonumber(domoticz.devices(device_freeze_alert).color) ~= 4 then domoticz.devices(device_freeze_alert).updateAlertSensor(domoticz.ALERTLEVEL_RED, freezingPoint.."°C") end

        elseif(tonumber(temperature) <= 1 and tonumber(freezingPoint) <= 0) and absHumidity < thresholdHumidity then
            domoticz.log("--- --- --- Givre peu probable --- --- ---", domoticz.LOG_DEBUG)
            message = "Givre peu probable malgré la température"
            if tonumber(domoticz.devices(device_freeze_alert).color) ~= 3 then domoticz.devices(device_freeze_alert).updateAlertSensor(domoticz.ALERTLEVEL_YELLOW, freezingPoint.."°C") end

        elseif(tonumber(temperature) <= 3 and tonumber(freezingPoint) <= 0 )then
            domoticz.log("--- --- --- Risque de Givre --- --- ---", domoticz.LOG_DEBUG)
            message = "Risque de givre"
            if tonumber(domoticz.devices(device_freeze_alert).color) ~= 2 then domoticz.devices(device_freeze_alert).updateAlertSensor(domoticz.ALERTLEVEL_ORANGE, freezingPoint.."°C") end
        else
            domoticz.log("--- --- --- Aucun risque de Givre --- --- ---", domoticz.LOG_DEBUG)
            message = nil
            if tonumber(domoticz.devices(device_freeze_alert).color) ~= 1 then domoticz.devices(device_freeze_alert).updateAlertSensor(domoticz.ALERTLEVEL_GREEN, 'Pas de givre') end
        end

        --if ((item.trigger == 'at 06:25') or (item.trigger == 'at 07:01') or (item.trigger == 'at 07:15')) and message ~= nil then
        if (item.isTimer) then domoticz.log("--- --- --- La règle qui a déclenché l'événement était :".. item.trigger .." --- --- ---", domoticz.LOG_DEBUG) end
        if (item.isTimer) and (message ~= nil) then
            domoticz.log("--- --- --- notification alerte Givre --- --- ---", domoticz.LOG_DEBUG)
            if SubSystem == nil then
                domoticz.notify(subject, message)
            else
            domoticz.notify(subject, message, domoticz.PRIORITY_NORMAL, domoticz.SOUND_INTERMISSION,"",  notificationTable(SubSystem))
            end
        end
end
}

Retrouver la derniere version de ce script sur



 Site Officiel :   github.com

2 thoughts on “[DZvents] alerte givre et calcul humidité absolue

  1. Bonsoir Vincent
    les devices température sont des dummy température simple
    le device humidité absolue est un dummy device Custom Sensor avec comme nom d’axe g/m3

  2. Salut Papoo,
    Peux tu préciser quel type de dummy device il faut utiliser pour les différentes variables stp ?
    Merci,
    Vincent

Laisser un commentaire